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FOREWORD 


This  is  me  Final  Report  on  contract  AF  30(802)  4287,  ARPA  Order  837,  entitled,  "Optical 
Processor  for  Synthetic  Spectrum  Data”.  The  sponsoring  agency  for  this  work  was  the  Air 
Force  Systems  Command,  Research  and  Technology  Division,  Rome  Air  Development 
Center,  Grlfflss  Air  Force  Base,  New  York  13440,  and  Project  Defender,  Advanced  Re¬ 
search  Projects  Agency  (ARPA),  Department  of  Defense,  Washington,  D.C.  20301  Mr. 
Robert  MacMillan  (EMASS)  was  the  RADC  Project  Engineer  and  Mr.  Morris  Wltow'  was  the 
ARPA  monitor  for  this  work.  Work  was  performed  on  this  contract  between  May,  1966,  and 
May,  1967,  at  the  Westtnghouse  Defense  and  Space  Center,  Surface  Division,  P.  O  Box 
1897,  Baltimore,  Maryland  21203 

Personnel  who  contributed  to  the  analysis,  design  and  documentation  performed  on  this 
contract  were  Messrs  B  L.  Drummond,  H.  S.  Fltzhugh  II,  T.  C.  Grlskey,  T.  R.  Hughes, 
G  S.  Ley,  C  A.  McGrew,  L.  A.  Nix,  M.  W.  Partin,  B.  S.  Smltn,  J.  Szostak,  and  J.  E. 
Thompson.  In  addition,  the  processing  and  analysis  of  recorded  data  was  performed  by  Dr. 

R  B.  Crane  and  Mr.  R  Fedorowlcz  of  the  Institute  of  Science  and  Technology,  University 
of  Michigan,  Ann  Arbor,  Mlchglan 


ABSTRACT 


A  film  recorder  for  high  resolution  radar  echoes  has  been  designed  and  fabricated.  It 
employs  a  modulated  cathode  ray  tube  trace  which  Is  imaged  onto  a  moving  film  Individual 
pul6e  echoes  are  written  across  one  dimension  of  the  film,  corresponding  to  the  direction  of 
the  CRT  sweep.  Successive  echoes  are  accumulated  slde-by-slde  In  the  orthogonal  direction 
The  data  Is  recorded  In  a  format  suitable  for  subsequent  processing  In  a  coherent  optical 
system  The  purpose  is  to  obtain  doppler  resolution  by  Integration  over  a  sequence  of  pulse*. 
This  Is  to  be  carried  out  simultaneously  for  each  range  resolution  element,  giving  an 
output  with  two  dimensions  of  resolution. 

The  specifications,  design,  fabrication,  and  check-out  of  the  film  recording  equipment 
are  described  in  this  report.  Simulated  data  has  been  employed  to  verify  the  overall  per¬ 
formance  of  the  device  A  review  Is  given  of  the  data  frames  recorded  and  the  test  results 
obtained  from  them.  The  system  can  record  data  In  two  different  formats,  to  allow  some 
flexibility  In  implementing  the  coherent  imaging  processor. 
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SECTION  I 


REPORT  SUMMARY 


1.  INTRODUCTION 

This  is  the  final  report  on  contract  AF  30  (602)4287,  ARPA  Order  837.  The  purpose  of 
this  contract  was  to  determine  the  feasibility  of  optically  processing  synthetic  spectrum  ra¬ 
dar  data  to  yield  a  high  resolution  range/doppler  map.  On  this  contract,  computer  software 
was  written  to  utilize  the  Univac  1108  and  the  Scientific  Data  Systems,  Inc.  SDS-910  digital 
computers  to  format  the  radar  data  for  subsequent  film  recording.  Digital  to  analog  conver¬ 
sion  equipment  and  a  timing  oscillator  for  system  synchronization  were  bought  from  SDS 
for  use  with  the  SDS-910  computer.  A  film  recorder  was  designed  and  built  by  Westlnghouse 
which  takes  video  information  and  synchronization  signals  from  the  SDS  equipment  and  dis¬ 
plays  these  on  the  face  of  a  high  resolution  cathode  ray  tube.  A  special  camera  was  fabricat¬ 
ed  to  Westlnghouse  specifications  by  Photogrammetry,  Inc.  Rockville,  Maryland.  This 
camera  contains  a  precision  film  drive  motor  and  associated  control  electronic  purchased 
from  Sequential  Electronics  Systems,  Inc. ,  New  York. 

2.  EQUIPMENT  DESCRIPTION 

A  photograph  of  the  film  recorder  is  shown  in  Figure  1,  in  which  the  various  sub-systems 
are  identified.  The  digital -to-analog  conversion  equipment  and  the  timing  oscillators  for 
system  synchronization  are  loci, ted  In  the  SDS-910  computer  classls  and  are  not  shown  In 
Figure  1. 


The  equipment  utilized  In  performing  optical  processing  of  synthetic  spectrum  data  Is 
shown  in  Figure  2.  A  brief  description  of  the  functions  performed  In  each  sub-system  of  the 

equipment  will  be  given  here.  Detailed  descriptions  of  the  functions  will  be  found  In  the  body 
of  this  report.  1 

The  UNIVAC  1108  digital  computer  is  used  to  perform  several  operations  on  the  radar  data 
and  to  assemble  It  into  a  format  compatible  for  the  SDS-910  computer.  The  radar  data  is 
calibrated,  the  responses  are  aligned,  pre-compensation  for  film  recorder  non-llnearitles 
can  be  applied  and,  finally,  the  proper  forma*  is  selected.  Additional  operations  on  the  data 
required  for  processing  in  data  format  II  (described  in  section  1-3)  are  also  performed  In  the 
Univac  1108  computer. 


The  SDS -9!°  digital  computer  Is  used  to  read,  store  and  output  the  data  In  synchronism 
with  the  sweep  of  the  cathode  ray  tube  (CRT)  beam.  The  DIGITAL-TO- ANALOG  (D  A) 
CONVERTER,  housed  In  the  SDS-910  computer,  accepts  the  9-blt  digital  output  and  converts 

R  7  a,s  Wh‘Ch  lS  the  Vldeo  slRTiat  10  ,he  CRT-  The  tIWNG  UNIT,  also  hous¬ 

ed  In  the  SDS-910  Computer,  generates  the  trigger  signals  and  control  frequencies  used  to 

synchronize  system  operation.  A  pulse  of  data  is  read  Into  the  SDS-910  core,  stored  there 
for  a  short  time,  and  then  read  out  and  written  onto  film,  prior  to  the  next  pulse  being  used. 

TnI™  c0nsls,s  01  four  subsystems,  the  CRT,  the  CRT  CONTROL  ELEC- 

™°N  CS/  thf  CAi1E,RA  the  FILM  DRIVE  WIT  w“h  its  control  electronics.  Timing  and 
control  signals  and  video  Informations  signals  are  supplied  to  the  FILM  RECORDER  The 

m  is  moved  continuously  at  a  precisely  controlled  rate  and  the  CRT  beam  Is  swept  repet - 

em  a  SfmC  tr7C  0,1  ,he  CRT  f3Ce-  °"cc  the  ‘"^Ration  Is  recorded  on  film 

tre  film  (after  development)  is  ready  for  processing  In  an  optical  correlator  On  this  mn. 

WerMty^of  Mlchlga*'18  prOCeSsed  by  ,he  InS,t,u,e  ot  Science  and  Technology  (1ST)  of  the  Un- 
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3.  BRIEF  DESCRIPTION  OF  FORMATS 


The  data  film  transparency  employed  with  an  optical  processor  (using  coherent  light)  must 
form  a  diffraction  pattern.  Particular  fringes  of  this  pattern  are  isolated  by  spatirJ  filtering 
and  further  transformed  by  lens  elements,  to  obtain  a  processed  output.  The  recording  for¬ 
mat  is  determined  largely  by  the  equlrcment  of  forming  an  interference  fringe  which  is 
modulated  by  the  radar  echo  data. 

To  give  the  required  fringes,  a  "spatial  carrier"  is  introduced  into  the  recorded  data. 

This  appears  as  an  underlying  sinusoidal  oscillation  of  the  film  transmissivity.  It  acts  as  a 
diffraction  grating  when  the  film  in  placed  in  a  coherent  optical  correlator.  The  presence  of 
the  spatial  carrier  also  allows  a  radar  signal  having  both  amplitude  and  phase  dependence  to 
be  recorded.  The  amplitude  is  used  to  control  the  magnitude  of  the  spatial  carrier  oscillations 
(e.  g.  the  envelope  of  the  waveform).  The  phase  is  used  to  shift  ti  e  position  of  the  carrier 
oscillations  by  a  fraction  of  a  cycle  (with  a  full  cycle  representing  360'’  of  phase).  This  is  the 
natural  meaning  of  a  phase  difference,  namely  the  shift  of  an  oscillatory  waveform  with  re¬ 
spect  to  a  standard  oscillation  of  the  fame  frequency.  In  tne  case  of  the  radar  data  film,  the 
standard  (phase  reference)  signal  would  be  a  uniform  sinusoidal  diffraction  grating  with’a 
spacing  constant  equal  to  one  cycle  of  the  carrier  on  the  film. 

In  addition  to  the  carrier,  a  bias  light  level  is  used  when  exposing  the  film  in  the  record¬ 
er.  This  bias  is  set  to  an  operating  point  on  the  linear  region  of  the  film  response  charact¬ 
eristic.  It  allows  the  carrier  oscillations  to  appear  as  fluctuations  in  transmissivity  about 
the  bias  level.  There  is  no  need  for  negative  values  of  transmissivity,  which  are  not,  in  fact, 
achievable  in  practice. 


Figure  2.  Optical  Processor  for  Synthetic  Spectrum  Data 
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The  recorder  can  produce  data  in  either  of  two  formats.  In  both,  each  line  on  the  film 
(corresponding  to  a  CRT  sweep)  represents  the  echo  from  one  radar  pulse.  The  difference 
lies  in  the  nature  of  the  data  written  along  each  sweep.  In  Format  1,  the  outputs  of  range 
channels  in  the  radar  directly  form  the  modulation  on  the  sweep.  That  is,  distance  along  the 
sweep  Is  proportional  to  range  across  the  tracking  window  of  the  radar.  In  Format  II,  the 
fourier  Format  I  data)  is  recorded  instead.  Here,  distance  along  the  sweep  is  proportional 
a  frequency  location  across  the  radar  signal  tiand. 

Atypical  film  recordingand processed  output  for  two  scatterers  is  shown  in  Figure  3.  A 
sample  of  the  recorded  film  (Format  II)  in  normal  size  is  shown  in  3(a).  In  3(b)  a  3ox  en¬ 
largement  of  this  recorded  film  is  shown  in  which  the  recording  process  is  more  apparent 
The  diagonal  bands  are  due  to  the  changes  in  spectrum  due  to  rotation  of  the  target,  which 
was  a  dumhell  configuration.  In  3(c)  the  processed  output  of  the  optical  correlator  is  shown 
in  which  both  the  zeroth  order  fringe  (on  the  optical  axis  of  the  correlator)  and  the  first 
order  fringe  are  shown.  The  zeroth  order  fringe  is  normally  removed  by  a  spatial  filter  and 
is  included  only  for  illustration.  The  offset  of  the  target  image  frnm  the  zero  order  is  caus¬ 
ed  by  the  presence  of  the  spatial  carrier.  The  light  from  the  background  bias  level  on  the 
film  is  concentrated  in  the  bright  zero  order. 

Additional  processed  results  are  given  in  Section  V. 

4.  SUMMARY  OF  CAPABILITY'S 


The  equipment  has  been  used  to  record  simulated  data  in  both  Formats  I  and  II.  The  tests 
winch  have  been  made  on  the  equipment  and  the  film  which  has  been  processed,  have  deter¬ 
mined  the  recorder  capabilities  which  are  listed  in  Table  A.  In  two  instances,  the  values  of 
the  parameters  arc  significantly  different  from  tin  design  goals  and  improved  data  recording 
is  possible.  Fu  st,  'lie  film  drive  unit  has  operated  wr'h  a  control  frequency  of  less  than 
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(b)FILM  FORMAT  (  X30  ENLARGEMENT  i 

Figure  3.  Typical  Film  Recording  and  Processed  Output- Format  II 


TABLE  l.  SUMMARY  OF  EQUIPMENT  CAPABILITIES 


CRT  spot  size 

0.  0015  Inches 

Number  of  "spots"  per  line 
of  recorded  data 

1024-2048  (selectable) 

Number  of  lines  (pulses)  per 

50  nm  frame 

500-1250 

CRT  line  sweep  time 

27-81  msec 

Spatial  carrier 

10-20  'pr/mm 

Maximum  number  of  pulses  which 
can  be  recorded  continuously 
(magnetic  tape  limited) 

3550 

Input  control  frequency  for 
film  drive 

600-60,  000 

600  Hz  where  714  Hz  was  the  design  value.  By  virtue  of  this,  the  maximum  density  of  lines 
along  the  film  was  approximately  25  lines  per  millimeter  rather  than  the  design  r'  vl  of  20 
lines  per  mllllimeter  (which  permits  1250  lines  per  50  mn.  frame  of  recorded  dr  a. )  To 
achieve  this  required  an  external  signal  source  be  used  instead  of  the  signal  from  the  TIMING 
UNIT. 

In  all  categories,  except  one  the  equipment  has  performed  at  or  beyond  design  specificat¬ 
ions.  The  one  exception  is  in  the  behavior  of  the  spatial  carrier  in  Format  II  recording.  In 
tills  format,  excessive  range  stdelobes  in  the  processed  data  are  produced  due  to  a  variation 
in  the  recorded  spatial  carrier  with  distance  along  the  trace. 

Although  the  recording  in  Format  II  has  corrections  to  be  made  to  it  in  order  to  give  best 
processed  output,  the  recordings  which  have  been  made  have  demonstrated  conclusively  that 
optical  processing  of  synthetic  spectrum  data  is  feasible.  The  results  given  in  Section  V  and 
the  conclusions  drawn  In  Section  V  I  elaborate  on  this. 

5.  APPLICATIONS  TO  SPACE  OBJECT  IDENTIFICATION  (SOI) 

To  apply  this  device  to  th?  SOI  problem,  the  equipment  corrections  to  give  recording  In 
Format  II  should  be  made.  Additional  processing  of  simulated  data  using  more  complicated 
grometrlcai  configurations  typical  of  targets  of  interest  should  be  performed.  This,  and  the 
processing  of  radar  data,  should  be  done  both  to  have  data  on  more  complex  targets  and  to 
have  more  data  from  which  to  make  observations  on  the  quality  of  processed  outputs  (espec¬ 
ially  In  having  a  time  record  available  on  which  to  base  observations).  These  recommendat¬ 
ions  are  discussed  in  Section  VI. 
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SECTION  II 


THEORY  OF  OPERATION 


A  target,  when  illuminated  by  a  microwave  radar,  modifies  the  incident  signal  in  the 
backscattering  process.  The  radar  recei"er  may  be  designed  to  sense  certain  of  these 
changes  or  modulations  In  order  to  give  Information  about  the  target.  The  most  commonly 
used  parameters  are  echo  time  delay  and  doppier  shift  these  are  linearly  related  to  target 
range  Jnd  range  rate,  respectively.  The  received  signal  power,  proportional  to  target 
radar  cross-section,  is  often  processed  as  well.  In  observing  space  objects,  one  would 
like  to  extract  much  more  detailed  information  about  the  target  Increasing  the  resolution 
capability  of  the  radar  can  provide  more  target-induced  data,  especially  when  the  resolution 
cell  becomes  smaller  than  the  target  extent.  When  the  Inherent  resolution  cell  Is  smaller 
than  the  body,  the  possibility  of  forming  a  "radar  image"  of  the  configuration  arises.  This 
section  discusses  radar  range 'doppier  Imaging  and  the  applicability  of  optical  techniques  In 
the  signal  processor.  Suitable  film  recording  formats  and  a  summary  of  critical  error  typec 
are  also  given. 

1  CHARACTERISTICS  OF  HIGH  RESOLUTION  ECHOES 

At  short  radar  wavelengths,  rigid  bodies  (such  as  satellites)  exhibit  so-called  "scattering 
centers".  That  Is,  the  radar  echo  can  be  associated  with  one  or  more  localized  regions  of 
the  target.  When  illuminated  with  a  very  short  pulse,  these  scattering  centers  produce  In¬ 
dividual  echoes  with  different  time  delays  (due  to  their  differences  in  range).  When  illumin¬ 
ated  with  a  single  frequency,  the  individual  echoes  superimpose  vectorlally  since  they  are 
not  resolved  In  time.  If  the  target  orientation  changes,  the  relative  phases  of  the  echoes 
are  altered  and  the  vector  sum  exhibits  a  beat  pattern.  The  rapidly  lobing  cross-section  of 
a  rocket  body  as  it  tumbles  is  a  common  example  of  this. 

A  range/doppler  imaging  radar  makes  use  of  both  the  range  differences  between  scatter- 
ers  and  the  relative  phase  changes  occurring  as  the  target  changes  orientation  by  a  small 
amount.  A  wlde-band  transmission  Is  employed  to  obtain  a  range  resolution  ceil  small 
compared  to  total  target  extent.  The  phase  information  can  be  gathered  with  a  single  radar 
site  by  observing  over  an  interval  of  time,  during  which  the  target  changes  aspect.  The 
phase  information  could  also  be  gotten  with  a  multisite  system,  where  the  separate  receivers 
each  view  the  object  at  a  different  angle,  at  the  same  time.  The  first  (monostatic)  type  of 
data  has  been  considered  on  this  contract. 

While  range  resolution  is  inherent  in  the  single  pulse  transmissions,  an  equivalent  doppier 
resolution  fs  obtained  only  after  processing  a  sequence  of  echo  pulses.  The  spectrum 
(Fourier  transform)  of  the  amplitude  and  phase  behavior  of  a  scatterer  gives  its  doppier 
resolution  pattern.  When  a  small  fraction  of  ?  total  target  revolution  period  Is  processed, 
the  doppier  offset  of  a  scatterer  is  proportional  to  its  location  along  an  axis  perpendicular 
to  the  Une-of-slght  This  property  aids  in  utilizing  the  radar  Image  as  a  spatial  map  or 
projection  of  the  echo  centers.  An  analytical  model  for  this  high  resolution  data  will  now  be 
given. 

a.  Form  of  Echo  Pulses 

A  pulse  stretch/collapse  system  would  undoubtedly  be  used  in  a  high  resolution  radar 
designed  for  observing  space  vehicles.  The  pulse  modulation  will  be  denoted  by  p  (t)  and  the 
r-f  center  frequency  by  The  transmitted  signal  will  be  represented  by  the  phasor 

s(t)  =  p(t)  exp  (2  -  i  fc  t)  (1) 
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where  a  is  the  amplitude  of  the  scatterer 

D  is  the  phase  shift  in  scattering 
r  is  the  delay  of  the  scatterer  (2R/C). 

When  the  scattering  center  has  frequency  dependence,  the  shape  of  the  echo  pulse  is  not 
exactly  that  transmitted,  since  the  scatterer  introduces  some  weighting  of  the  spectrum. 
However,  this  distortion  is  slight  for  a  radar  pulse  of  less  than  20%  bandwidth,  and  will 
not  be  considered  here. 


An  extended  solid  body  usually  exhibits  several  scattering  centers,  whose  characteristics 
(amplitude,  phase,  polarization  and  frequency  dependence)  depend  on  the  direction  of  the 

line-of-sight  The  total  echo,  e(t),  for  a  given  orientation  is  a  summation  of  single  scatterer 
terms: 

N  10 

n 

c(t)  =  X  a  c  s(t -  r  ) 
n  =1  " 


2nifct  N  i(0  -2-f  r n) 

1  e  X  V  "  C  p(t-  r  )  (3) 

n  =  l  n  n 

where  Is  the  amplitude  of  the  n**1  scattering  center 

Phase  shift  in  scattering  for  the  n**1  scatterer 
rn  is  the  delay  of  the  nlil  scatterer 
N  is  the  total  number  of  scattering  centers 

S  is  the  transmitted  signal,  S(t)  =  p(t)e  2’rWct 
fc  Is  the  center  frequency  of  the  pulse 

p  is  the  (baseband)  modulation  of  the  radar  puluc. 

Each  term  In  Equation  (3)  can  usually  lx?  associated  with  a  localized  discontinuity  In  the 
target  surface  curvature  or  a  section  of  surface  normal  to  the  line-of  sight 


8 


In  the  radar  receiver,  this  echo  Is  processed  to  achieve  pulse  collapse.  For  a  "matched” 
receiver,  the  echo  is  correlated  with  the  conjugate  of  the  transmitted  signal.  Mathematically, 
the  radar  response  r{  v  )  Is  given  by 

r (r  )  =  J  e(t)s  (t-  r  )  dt  (4) 

This  reduces  to 

N  10  2elf  (  r  -  r  )  . 

r(r)»£aene  c  n  /  p(t-  r  )p(t-r)dt 

n-1  n  n 

The  integral  in  the  last  equation  is  the  matched  output  response,  r  (  r  -  r  ),  of  the  radar. 
There  are  N  such  terms,  one  for  each  scattering  center.  Therefore,  the  r3dar  response 
of  the  extended  target  is 

N  10  2rif  (r  -  r  ) 

r(  r  )  1  I  a  e  ne  c  n  r  ( r  -  r  )  (5) 

n  =  l  n  on 

where  rQ  is  the  matched  output  response,  defined  by 

ro  (  X  )  s/p(t)  p  (t  -  X  )  dt 

f  is  the  radar  center  frequency 


a  .  0  ,  ,  N  are  as  defined  above  (Equation  (3)) 

n  n  n 


is  the  running  delay  variable.  The  individual  response  peaks  occur 
when  r  *  r  ;n  =  l,  ...N. 


The  Synthetic  Spectrum  radar,  which  is  the  data  source  for  the  optical  recording  and 
processing  systems,  measures  the  target  response  in  a  bank  of  range  channels.  These  are 
spaced  at  equal  intervals  in  delay  (or  range).  A  sample  of  the  output  response  sign*., 
obtained  in  each  of  these  channels.  This  sample  Is  a  complex  number  In  polar  form;  that 
Is,  the  amplitude  and  phase  are  measured  and  recorded. 

lArgrm 

Output:  r(  rm)  =  |rjc  ;  rm  *  tq  *  m  rg,  m*  -43 . 42 

tq  =  gross  target  delay  established  by  the  real  time  tracker 
i  -  sampling  Increment 


The  single  pulse  information  (84  amplitudes,  84  phases,  and  much  auxiliary  data  like 
range,  range  rate,  azimuth,  elevation,  etc.  )  is  multiplexed  out  of  boxcar  stores  onto 
magnetic  tape  during  the  interpulse  period  This  data  is  in  digital,  not  analogue,  form. 
Amplitudes  are  recorded  on  a  logarithmic  scale  with  about  80  db  range  using  nine  bits. 
Phases  are  recorded  in  the  Interval  (0,  2ir)  using  eight  bits. 

In  the  past,  data  processing  operations  have  been  carried  out  on  a  large  digital  computer 
(IBM  7094  or  UNIVAC  1108).  The  amplitude  and  phase  information  Is  first  corrected  using 
calibration  data  (to  compensate  for  wlde-band  channel  differences  and  log-amp  and  phase 
detector  non-linearities).  The  relative  phases  between  radar  echoes  are  then  modified  to 
remove  as  much  of  ihe  phase  modulation  due  to  target  translation  as  possible.  Three 
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techniques  for  doing  this  are  discussed  in  References  1  and  2.  •  After  a  good  "pulse  align¬ 
ment  by  erne  of  these  methods,  the  pulse-to-pulse  modulation  of  the  high  resolution  echoes 
appears  to  arise-  from  target  rotation  about  a  fixed  point.  This  modulation  is  discussed  in 
the  next  section 

o  Pulse -To-  Pulse  Phase  Modulation 


C  onstder  one  of  our  target  scattering  centers  moving  with  respect  to  the  fixed  body 
phase  reference  Its  delay,  r  ,  changes  between  pulses  Over  a  small  interval  of  target 
rotation,  say  0.  1  radian  i3.  7  ),nthe  scatterer  delay  may  be  approximated  by  an  average 
value  and  a  term  linear  in  time. 

T  .  ,  (time  origin  at  middle  of 

n  3  n  '  n  small  rotation  Interval)  (6) 


Our  moving  scatterer  is  represented  by  one  term  in  Equation  (5),  namely 

10  2slf  (  r  -  T)  .  „  ,  -2  rif  „  t 

ae  e  c  'r  (  r  -T-  v  t)  e  c 


(7) 


where  the  subscript  n  labeling  the  scatterer  has  been  dropped.  Suppose  the  scatterer  moves 
by  one  tenth  resolution  cell,  i  e.  by  a  delay  ut  -  (10B)'1  where  B  is  the  signal  bandwidth 
Letting  the  center  frequency  f  the  phase  change  in  the  measured  echo  (contained  in  the 
last  factor  of  expression  (7)  Is  A 

'2'fcUt  *  ‘l'fef  1  *  2saT0"  ;  A  *  %  bandwidth. 


For  example,  if  A  is  10T  the  phase  change  Is  a  full  cycle.  Associated  with  this  remarkable 
Phase  change  is  a  small  shift  In  the  "short  pulse  envelope"  r  (r-T-vt)  which  amounts  to 
ai>out  one  tenth  its  3  db  width.  0 


i(0 


The  other  factors,  ae  '  c^  are  senslbiy  constant  over  the  small  rotation  inter¬ 

val  Therefore,  the  dominant  effect  of  the  relative  rotational  motion  of  the  individual  echo 
centers  Is  a  rapid  r-f  phase  change.  This  progressive  phase  shift  is  equivalent  to  a  doppler 
frequency  A  second  dimension  of  resolution  can  be  obtained  by  frequency  analysis  of  pulse- 
to-pulse  echo  data,  as  described  in  the  next  section. 


c.  Processing  for  Range/Doppler  Resolution 

The  echo  data  is  Inherently  two-dimensional.  In  Equation  (7)  the  two  variables  are: 


r  :  representing  single  pulse  delay  across  the  84  channel  "range  window" 

t  :  contained  in  each  puise-to-pulse  scatterer  delay  term,  r  *T  ♦  u  t 

n  n  n  ’ 

measuring  time  relative  to  the  pulse  in  the  middle  of  a  data  block 
(t  0  for  this  pulse). 

There  are  also  two  "variables"  or  parameters  associated  with  each  echo  center: 

T  :  the  delay  of  the  scatterer  relative  to  the  fixed  phase  reference  point, 
at  t  *  0. 


:  the  rate  of  change  of  delay  of  the  scatterer,  at  t  =  0, 

We  would  like  to  obtain  a  display  with  coordinates  delav  and  delay  rate  in  which  the  energy 

from  the  nlh  scatterer  is  concentrated  near  the  point  (T  ,  v  ).  The  coordinate  r  (relative 

— — -  n  n 

•  References  are  given  on  pag*  26. 


10 


delay)  already  has  this  property  In  that  the  "short  pulse  envelope",  r  (  r  -T  )  has  a  peak 

o  n 

when  Its  argument  Is  zero,  1.  e.  when  r  *  T  .  The  Fourier  transform  at  the  data  In  the 
t-dlrectlon(a  one-dimensional  transform)  concentrates  the  echo  energy  In  the  new  variable 
f  This  is  easily  seen  by  considering  the  expression: 


u(  r  ,  f) 


2r  If  (  r  -  T) 

V 


r -T- v t 


2-l(f-f  u 
c 


)  t. 


(8) 


with  tj  •*  €(At),  At  is  the  pulse  repetition  interval. 


A  sum  is  used  rather  than  an  integral  because  the  radar  transmits  discrete  pulses.  The 
slight  amplitude  modulation  of  the  term  in  brackets  will  be  neglected  (this  modulation  is 
caused  by  the  shift  of  the  envelope  function  r  ,  as  v  t  changes  values).  The  result  is 


1(5  2811  (  r  -T) 

u(  r  ,  f)  e»  ae  e  c  r  (  r  -  T) 


sin  [(2L  4  1)  8  (f  -  f  „  )At] 
sin  [8(f  -  fcu)At] 


(9) 


Suppose  that  the  radar  has  a  rectangular  power  spectrum.  An  explicit  form  then  results 
for  the  range  resolution  factor  r  (  r  -  T). 

B  2 

r  2sift  slngBt 

rQ(t)  =  /  e  df  -  rt  (10) 

-B/2 


Thus,  the  product  of  two  "resolution  patterns"  occurs  for  the  typical  echo  center: 

../  n_  ,r  l0r  2slfc(  r*T)  sinrB  (  r  -T)  sin  [(2L  ♦  1)  i  <f-fP„)At ] 

- v  '  —  v— - ' 

amplitude  &  phase  range  resolution  doppler  resolution 

The  peak  occurs  at  r  =  T  and  f  *  f^u  .  The  first  zero  in  the  range  pattern  is  at  (  r  -T)  =  g. 

The  first  zero  in  the  doppler  pattern  is  at  (f-fr»>  )  «  rer  }  =  f  Here  ^  ls  total 

processing  interval  since  it  is  the  number  of  pulses  times  tnetr  spacing  The  result  above 
is  perfectly  obvious,  namely  that  the  resolution  in  delay  time  is  roughly  the  reciprocal  of 
the  signal  bandwidth  and  that  the  resolution  In  doppler  frequency  is  roughly  the  reciprocal  of 
the  Integration  time. 

The  purpose  of  the  optical  processor  is  to  perform  the  doppler  frequency  analysis  (one- 
dlmensional  Fourier  transform)  simultaneously  for  ait  range  channels.  It  is  easily  recog¬ 
nized  that,  if  the  spectra  of  the  echo  pulses  are  used  in  lieu  of  the  pulses  themselves,  a 
processor  employing  a  two-dimensional  transform  gives  simultaneous  range  and  doppler 
resolution.  Both  these  configurations  are  well-known  in  optical  signal  processing  and  are 
briefly  described  in  the  next  section-  The  data  film  for  the  processor  requires  a  special 
format,  since  the  radar  data  initially  has  both  amplitude  and  phase  modulation.  To  record 
this  information  on  film  use  is  made  of  a  "spatial  carrier",  simitar  to  the  holographic  re¬ 
cording  technique  developed  at  University  of  Michigan.  The  format  and  functional  form  of 
the  writing  beam  intensity  are  discussed  in  a  following  section.  Finally,  the  effects  of 
errors  in  recording  are  summarized. 
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2  BASIC  OPTICAL  PROCESSOR 


a.  Introduction 


The  optical  processor  for  SOI  radar  data  should  transform  echo  Information,  recorded 
on  a  film  transparency  to  a  range  doppler  image.  The  radar  data  derives  from  the  echo 
pattern  for  each  radar  pulse.  The  transparency  is  made  by  a  CRT  film  recording  system. 

The  information  from  a  single  pulse  is  modulated  as  brightness  variations  along  a  sweep  of 
the  writing  beam  Successive  pulses  are  stacked  side  by  side  along  the  film.  The  developed 
transparency  is  processed  in  a  coherent  optical  system.  In  order  to  carry  out  the  opera¬ 
tions  outlined  in  the  iast  section,  the  Fourier  transforming  properties  of  lenses  are  used 
The  configuration  chosen  for  the  processor  has  a  bearing  on  the  film  format  and  is  d(scussed 
here  for  this  reason. 

b  Processor  Configuration 

The  processor  consists  of  a  source  and  collimating  assembly,  a  input  film  plane,  a 
two-dimensional  transforming  lens,  a  spatial  filter,  a  one-dlmensionai  transforming  lens, 
a  relay  lens  and  the  output  data  plane.  The  one -dimensional  transform  can  be  removed  by 
pre-transformlng  the  radar  data  before  the  input  transparency  is  written.  These  systems 
are  shown  schematically  in  Figure  4 


Figure  4.  Basic  Optical  Processor  Schematic 
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The  source  is  conveniently  a  laser,  which  provides  suitable  coherence  with  reason¬ 
able  intensity  The  collimating  assembly  includes  a  condensing  lens  followed  by  a  pinhole 
and  a  collimator.  The  pinhole  is  used  to  remove  "noise"  in  the  intensity  distribution  of  the 
laser  output  caused  by  reflections  from  the  mirrors  and  walls  of  the  discharge  tube  The 
collimator  produces  a  plane  wave,  covering  the  input  film  aperture,  from  the  light  leaving 
the  pinhole. 

The  input  film  plane  contains  the  radar  dara  transparency  immersed  in  an  index- 
matching  liquid  cell  or  "liquid  gate".  This  is  done  to  compensate  for  optical  path  differences 
caused  by  film  thickness  variations.  Successive  radar  pulses  are  recorded  on  separate 
traces,  side-by-side  across  the  film.  They  accumulate  in  the  direction  marked  "time"  in 
Figure  4  The  coordinate  in  the  other  direction  can  be  relative  range  (as  in  an  A-scope 
display)  or  the  Fourier  transform  variable,  wave  number.  These  two  possibilities  are 
Formats  1  and  11,  respectively. 

The  spherical  lens  following  the  input  data  plane  produces  the  Fourier  transform  of 
the  information  on  the  transparency.  This  two-dimensional  transform  appears  in  the  jack 
focal  plane  of  the  lens  The  data  must  be  recorded  in  a  way  which  permits  one  to  select  the 
desired  information  from  the  background  light  by  means  of  a  spatial  filter.  To  accomplish 
this  a  spatial  carrier  is  introduced  into  the  brightness  variations  along  the  data  sweep.  The 
spherical  lens  then  produces  fringes  corresponding  to  this  sinusoidal  diffraction  grating. 

The  spatial  filter  in  the  bark  focal  plane  of  the  transforming  lens  blocks  the  zero 
order  light  which  represents  the  background  transmitted  by  the  transparency.  An  opening 
in  this  "filter”  allows  a  region  about  the  first  order  fringe  of  the  spatial  carrier  to  pass 
through.  The  Fourier  transform  of  the  puise-to-pulse  phase  modulation  into  its  doppler 
components  is  achieved  in  this  plane  If  the  film  is  in  Format  1,  an  additional  one-dimen¬ 
sional  tiansform  to  reconstitute  the  range  resolution  is  done  by  the  cylindrical  lens.  Format 
fl  dees  not  require  this  lens.  In  Format  11,  the  spherical  lens  forms  the  range  pattern  as 
the  Fourier  transform  of  the  wavenumber  information  on  the  input  film. 

.  he  relay  lens  may  be  included  to  enlarge  the  output  image  and  focus  it  on  the  final 
observation  plane.  This  might  be  a  ground  glass  screen  or  a  film  In  a  camera,  if  a  perma¬ 
nent  record  is  desired. 

3  DATA  FILM  FORMAT  REQUIREMENTS 

a.  General  Considerations 

The  transmissivity  of  the  film  is  controlled  by  the  CRT  trace  exposure.  This  allows 
only  intensity  variations  of  the  light  through  the  transparency.  The  phase  variations  due  to 
the  film  itself  cannot  be  modulated  and  are,  in  fact,  smoothed  out  by  the  liquid  gate  The 
radar  data  contains  both  amplitude  and  phase  modulation  and  both  a.*e  necessary  in  pro¬ 
cessing  for  a  range, doppler  image.  Thus,  the  'irst  consideration  is  to  convert  the  com¬ 
bined  amplitude  phase  modulation  of  the  data  to  amplitude  or  intensity  fluctuations  only. 

It  is  also  necessary  to  record  using  a  section  of  the  film  transfer  curve  that  does 
not  severely  saturate  Saturation,  being  a  non  linearity,  will  produce  intermodulation 
between  the  Individual  echo  centers  and  give  rise  to  ghost  scatterers 

Finally,  an  underlying  grating  structure  (i.  e  a  periodicity)  should  be  introduced 
This  allows  the  background  light  to  be  spatially  filtered  leaving  a  diffraction  fringe  which 
carries  the  radar  data  All  three  of  these  considerations  can  be  achieved  by  recording  the 
data  as  modulated  on  a  spatial  carrier  with  a  bias  level  added  to  prevent  negative  ampli¬ 
tudes  A  combination  of  setting  the  bias  and  sealing  the  maximum  amplitude  is  used  to 
obtain  a  desirable  operating  range  on  the  film  transfer  curve 


13 


h — 0 — i 


(0) 

ECHO  PULSE  MODULATION  SPECTRUM 


lc(f)| 


(b) 

SPECTRUM  WITH  SPATIAL  CARRIER  AND  BIAS 


Figure  5.  Frequency  Specta  For  Data  Traces 
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b  Effect  of  Spatial  Carrier  and  Bias  Level 


As  pointed  out  above,  the  film  can  reproduce  only  an  amplitude  modulated  non- nega¬ 
tive  signal.  Simple  modulation  theory  allows  one  to  analyze  a  method  for  converting  Me 
radar  6tgnal  (having  both  A-M  and  P-M)  to  a  real  unipolar  waveform.  The  radar  echo  to  be 
recorded  will  be  denoted  by  r(  r  ).  This  has  a  spectrum  centered  on  zero  frequency  of  band¬ 
width  B,  as  represented  In  Figure  5  .  The  echo  amplitude  lr(  r  )l  and  phase  Arg  r  (  r  ) 
are  used  to  modulate  a  carrier.  The  carrier,  denoted  by  f  ,  Is  chosen  greater  than  B  '2 
giving  rise  to  a  waveform  0 

g(r)  *  |  r  ( r )  |  £cos  2h(qt  +  Arg  r  (r  )J  (11) 

This  real  function  has  positive  and  negative  frequency  bands  of  width  B  centered  on  t  f  as 
In  Figure  5.  It  can  change  sign  because  of  the  cosine  function.  To  make  a  positive  Signal, 
one  may  add  a  d.  c.  bias  R^  greater  than  the  maximum  value  of  the  signal  envelope. 

f(r)  *  Ro  ♦  |r  (  r  )|  cos[2rfQ  r  ♦  Argr  r  )]  ;  R  >max|r  (r)|(12) 

In  the  frequency  domain,  Rq  appears  as  a  line  at  zero  frequency. 


Suppose  a  film  Is  written  with  a  trace  whose  amplitude  transmission  function  varies 
as  Equation  12  Let  It  be  Illuminated  at  normal  Incidence  with  coherent  light  According  to 
simple  physical  optics,  the  Fraunhoffer  pattern  of  the  film  Is  formed  by  the  emerging  light. 
This  pattern  has  three  parts  which  correspond  to  the  three  sections  of  the  spectrum  of 
f  (t).  They  emerge  at  different  angles  as  In  Figure  6  and  can  be  separated  by  a  spatial 
filter  In  the  focal  plane  of  an  objective  lens  behind  the  film.  The  first  order  spectrum 
represents  the  function 

Y  |  R  (y)|cxp  [2nlfQy  ♦  I  Arg  R  (y)l  *  y  R  (y)  e  2"V 

V  =  v  R(y)  >  r  (  i-  ) 

2 

where  the  delay  variable  r  has  been  replaced  by  a  proportional  distance  variable  y.  This 
function  has  the  desired  factor,  r  (y),  which  carries  the  amplitude  and  phase  of  the  radar 
data  In  the  correct  form  This  part  of  the  focal  plane  distribution  enters  the  remaining 
section  of  the  optical  processor. 

The  film  s,  of  course,  two-dimensional  and  the  carrier  can  be  Introduced  In  either 
the  time  (x -dimension)  or  the  single  pulse  variable  (y -dimension)  or  even  In  a  skewed 
direction  The  single  pulse  radar  echoes  are  defined  by  84  or  less  samples  (complex 
numbers)  depending  on  the  number  of  radar  range  channels  used.  However,  more  samples 
are  needed  ln  writing  the  film  to  represent  the  spatial  carrier.  The  desirable  number  of 
pulses  to  be  used  In  one  Integration  frame  is  from  about  a  hundred  up  to  one  thousand  A 
50x50  mm  square  on  70  mm  film  contains  one  of  these  frames  of  radar  data.  A  spatial 
carrier  frequency  of  10  line-pair  per  mm  (cycles  mm)  or  greater  Is  desirable  to  get  the 
data-carrylng  fringe  out  of  the  "processor  noise  region  around  the  direct  beam. 


The  coherent  light  used  In  processing  the  radar  data  must  Interact  In  a  linear 
(amplitude)  manner  rather  than  as  addition  of  Intensities  (power).  For  proper  operation 
the  amplitude  transmission  coefficient,  T  of  the  film  (after  developing)  should  be  propor¬ 
tional  to  the  function  In  equation  (12), 


♦  R  (X  ,  y)  cos  [ 


*  Arg  R  (x.y)  ] 


15 


PROCESSOR 

entrance 

PUPIL 

tSPATlAL 

FILTER) 


Figure  G.  Diffraction  Pattern  Introduced  by  the  Spatial  Carrier 

"'here  :  light  intensity  transmitted  by  the  recording 

Io;  intensity  of  the  coherent  read-out  beam  in  the  optical  processor 
R«  :  bias  level 

|R(x,y)| :  envelope  of  the  radar  data  translated  to  film  coordinates 
Arg  R(x,  y) :  phase  of  the  radar  data  translated  to  film  coordinates 
f0  :  spatial  carrier 

y  :  distance  coordinate  on  the  film  proportional  to  fine  range  (or  delay) 
on  a  single  radar  echo 

X  :  dirtance  coordinate  on  the  film,  orthogonal  to  y,  which  takes  on  sucessive 

values  for  successive  radar  pulses  (I.  e.  proportional  to  time). 

If  the  spectra  of  ttr  radarpulses  are  recorded,  the  variable  y  Is  proportional  to  frequency 
measured  across  the  radar  signal  hand,  rather  then  to  delay  with  In  the  echo  response. 

The  transfer  characteristic  between  the  grid  drive  of  the  CRT  and  the  film  transmission 
coefficient  determines  the  operating  points  This  is  discussed,  with  reference  to  experi¬ 
mental  data,  in  Section  IV. 
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c.  Sampling  and  Spot  Size  Effects 

The  grid  drive  of  the  CRT  is  modulated  by  digital  information  which  has  been  D/A 
converted  from  a  magnetic  tape  The  individual  data  samples  are  nine-bit  quant'tles.  The 
actual  echo  data  shows  symmetric  deviations  about  a  central  bias  level  Thus,  the  max¬ 
imum  amplitude  is  half  t*ie  quantizing  range,  and  equal  to  255  times  the  quantizing  step.  This 
"contrast"  is  equivalent  to  48db.  A  signal  30dbbeiow  the  maximum  would  be  18  db  above  the 
quantizing  unit.  Its  oscillation  would  be  represented  by  the  first  8  quantizing  levels.  The 
rather  fine  subdivision  afforded  by  nine  bits  gives  some  freedom  in  operating  with  a  mark¬ 
edly  non-linear  transfer  characteristics,  when  pre-compensation  of  the  tape  data  is  used 

A  given  CRT-trace  (e  g.  radar  puise'  contains  many  of  the  individual  nine-bit  samples 
the  sample  density  must  be  greater  than  that  <ji  the  original  radar  data  because  the  mformat- 
ion  is  now  in  sidebands  of  the  spatial  carrier.  The  "sampling  theorem"  gives  a  maximum 
spacing  of  Ay*  j_  where  W  is  the  highest  spatiai  frequency  in  the  signal  However,  this 
2W 

rate  is  not  suitable  when  a  finite  section  of  a  signal  is  *o  be  represented.  A  higher  density 
is  needed  to  reduce  the  "truncation  errors"  as  developed  in  Reference  3. 

A  representation  of  this  effect  in  the  frequency  domain  is  shown  in  Figure  7.  The 
spectrum  of  the  band-limited  waveform  to  be  recorded  is  represented  in  Figure  7.  When 
this  v  iveform  is  sampled  with  a  spacing  A  y,  the  associated  spectrum  becomes  an  infinite 

set  of  replicas  separated  by^^  as  (n  Figure  7.  The  space  between  successive  replicas 

is  referred  to  as  the  guard  band.  Because  the  waveform  is  land  limited,  its  time  repre¬ 
sentation  cannot  be  identically  zero  outside  some  finite  interval.  Recording  only  a  segment 
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Figure  7.  Effects  of  Sampling 
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of  this  signal  thus  changes  the  spectnu:'  again.  If  N  samples  are  recorded,  the  signal  is 
effectively  gated  by  a  rectangular  wave  NAy  in  width.  The  frequency  domain  effect  Is  that 
the  set  of  replicas  (In  M3b)  Is  convolved  with  the  transform  of  the  gate  The  latter  scanning 
function  Is  sin  *  f  NAy  which  has  first  zero  crossings  at  f=  J_  .  It  is  Important  that  the 
i  IN  A  y  NAU 

guard  band  be  several  times  wider  than  the  main  lobe  on  this  sin  X  function.  For  example, 

X 

the  third  sidelobe  of  the  scanning  function  is  some  20  db  down  while  the  fifth  Is  suppressed 
about  25  db.  Sampling  at  the  rate  3f0  (f0  being  the  spatial  carrier)  gives  a  wide  guard  band. 
The  convolution  effect  Is  shown  In  Figure  7. 

Actually,  point  samples  .'.re  not  produced  by  the  D  A  converter  which  drives  the  CRT 
Intensity  control  Instead,  the  voltage  adjusts  to  the  digital  input  rapidly  (5  usee  or  less) 
and  remains  at  the  steady  level  until  the  next  sample.  The  effect  Is  of  a  "box  car'  where 
each  sample  Is  convolved  with  a  rectangle  of  width  Ay  .  the  space  between  samples.  This 
means  that  the  spectrum  is  multiplied  by  an  envelope  sin  *  fAy  as  In  Figure  7.  The  spot 

s  fAy 

size  of  the  CRT.  the  camera  lens,  and  the  film  itsell  Introduce  further  frequency  transfer 
functions.  These  should  Ik>  maintained  as  broad  as  possible  (l  e  resolution  not  degraded) 
to  minimize  spectral  distortion.  The  spot  size  appearing  on  the  film  and  the  associated 
spatial  frequency  transfer  function  are  reviewed  In  Section  IV 

To  recapitulate,  the  spatial  carrier  f0  must  l>e  chosen  larger  than  the  Inherent  data 
l»andwldth  so  that  the  sldelxinds  (fringes)  will  l>e  separated  The  digital  sampling  rate  must 
be  significantly  greater  than  2f0  Is  more  than  adequate 

The  maximum  Inherent  Ixindwidth  of  the  single  pulse  data  (when  on  film)  Is  easily 
found  The  actual  radar  response  has  a  bandwidth  not  greater  than  the  reciprocal  of  the 
resolution  cell  length.  When  this  cell  length  Is  expressed  tn  mm  on  the  film.  Its  reciprocal 
Is  the  spatial  Iwtndwldth  of  the  envelope  and  phase  modulation  of  the  data  trare  There  will 
be.  at  most,  84  cells  tn  50  nm  Thus  the  data  bandwidth  Is 

dj 

I  68  cvcles  nin  (line  pair  nm) 

50  nun 

because  an  actual  target  will  shift  alxjut  somewhat  In  the  real  time  range  gate,  some  rhan 
nets  on  the  ends  will  not  contain  useful  Information  Thus.  84  range  cells  Is  an  upper  limit, 
with  CO  75  cells  living  more  likely  The  data  bandwidth  reduces  tn  only  I  2  cycles  nm  with 
60  range  cells  If  the  spectrum  form  of  the  data  Is  used,  one  would  have  the  same  Inherent 
bandwidth  In  this  rase,  variations  having  n  42.  n  41.  .  n*40.  n*41  rycles  tn  50  nm 

(n  lielng  the  number  of  spatial  carrier  cycles)  are  required  to  represent  channels  42.  -41, 

.  .  40.  41  Again  the  maximum  difference  frequency  Is  84  eyries  tn  50  nm  or  I  68  cycles/ 
mm  This  inherent  bandwidth  Is  for  the  single  pulse  (y)  dimension  or  direction  or  direction 
on  the  film 

In  the  x-directlon.  the  bandwidth  of  the  doppler  variations  of  the  radar  signals  appears 
as  spatial  frequencies  This  Is  the  direction  of  film  motion  in  the  camera  The  data  band¬ 
width  in  this  dimension  depends  on  the  taryet  as  well  as  on  the  radar  waveform  and  prf 
It  Is  possible  (indeed,  a  common  otturance)  that  the  radar  prf  Is  too  low  for  unamblgous 
data  samples  Thus,  the  doppler  (\  direction)  bandwidth  will  be  taken  as  the  maximum  allow 
ed  by  the  prf  A  phase  change  lying  between  -180  and  *180  Is  allowable  from  one  pulse  to 
the  next  The  relative  difference  Is  180-  (-180)  or  one  cycle  Thus,  If  1000  pulse  traces  are 
w  ritten  In  the  50ntm  aperture,  one  could  have  a  spatial  doppler  frequencies  spanning  20 
cycles  nm  Lower  frequency  excursions  on  the  film  can  be  obtained  by  increasing  the  dis¬ 
tance  between  I  races  or  by  repeating  the  same  pulse  data  In  a  sequence  of  redundant  lines 

The  spatial  carrier  has  been  Introduced  in  the  single  pulse  trace  since  the  data  band 
width  there  Is  raiher  low-  In  synthetic  aperture  radar  processors  the  carrier  usually  appears 
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In  tne  time  ordoppler  dimension  Tltese  systems  have  many  more  range  cells  and  hence 
more  cross  film  bandwidth  than  in  the  present  case  In  tlie  present  case,  the  spatial  carrier. 
If  inserted  in  the  x -direction,  would  be  unrealistically  high,  considering  the  resulting  sam¬ 
ple  density  anti  the  tube  spot  size  We  thus  restrict  further  fc-mat  consideration  to  a 
carrier  In  the  y-direction  with  either  the  time  domain  or  frequency  domain  representation 
of  the  echoes. 


4  SUMMARY  OF  ERROR  SENSITIVITIES 
a  Introduction 

Positional  errors  in  the  training  beam,  film  drive  deviations  from  constant  speed, 
phospnor  and  film  noise,  and  ti*e  transfer  characteristics  of  the  CRT,  lens,  and  film  will 
all  causf  errors  In  the  transparency’  When  the  film  is  used  in  the  coherent  optical  processor, 
it  acts  as  an  extended,  two-dimensional  source.  That  is,  the  plane  wave  of  coherent  light 
impinging  on  it  is  modulated  by  the  local  transmissivity  of  each  region  of  the  film.  The  emer¬ 
ging  light  superimposes  coherently  (at  a  large  distance  from  the  film)  to  form  an  Interference 
pattern.  This  pattern  is  the  two-dimensional  Fourier  transform  of  the  amplitude  transmiss¬ 
ivity  of  the  film. 

The  presence  of  the  spatial  carrier  allows  one  to  separate  a  fringe,  carrying  the  ra¬ 
dar  data  modulation,  from  background  light  in  the  direct  beam.  However,  the  spatial 
carrier  has  the  effect  of  converting  small  errors  in  the  location  of  the  writing  beam  to 
relatively  large  phase  fluctuations  in  the  echo  data.  For  example,  if  the  carrier  is  10  line 
pa,r  per  mm,  one  cycle  (360  )  occupies  i  mm,  and  a  3.  60  phase  error  corresponds  to 

10 

jnly  J_  mm  (one  micron).  The  phase  error  due  to  a  positional  error  is  proportional  to  the 
1000 

spatial  carrier.  This  carrier  frequency  cannot  be  arbitrarily  reduced,  however,  becasue 
the  data-bcarlng  fringe  must  be  deviated  away  from  the  direct  beam  and  the  processor  noise 
region  near  zero  spatial  frequency. 


The  film  acts  somewhat  like  a  two-dimonstonal  antenna  composed  of  individual  elements. 
Both  have  radiation  patterns  given  by  the  Fourier  transform  of  their  aperture  distribution 
Analysis  of  error  effects  can  be  based  on  this  analogy.  The  film  Is  regarded  as  consisting 
of  cells  (elements)  each  with  amplitude,  phase,  and  positional  errors.  There  are  two  ways 
of  picking  these  elementary  cells.  The  first  is  to  regard  each  sample  point  (9-bit  number) 
transfered  from  tape  tc  film  as  a  separate  element.  There  are  between  1000  and  2000  of  these 
points  In  a  trace  and  between  100  and  1000  traces  or  lines  in  a  data  frame.  In  Appendix  II,  an 
error  analysis  is  carried  out  from  this  standpoint. 


The  second  viewpoint  is  to  think  of  each  trace  as  composed  of  radar  resolution  cells  or  in¬ 
dependent  spectral  components,  depending  on  the  format.  There  are,  at  most  84  such  cells 
in  a  trace,  with  60  being  a  more  likely  number.  Again,  there  are  100  to  1000  traces.  In  this 
second  approach,  the  true  fine  structure  within  the  cell  Is  by-passed  and  a  resultant  ampli¬ 
tude  and  phase  error  applied  to  the  cell  as  a  whole.  If  the  spatial  carrier  is  10  line  per/mm, 
and  60  cells  are  recorded  in  a  trace  there  will  be 


50  mm  10  cycles 

60cells  mm 


8.  3  cycles/cell 


There  is  little  difference  in  the  two  view  points  for  Format  I  data,  since  there  is  no 
integration  across  sub-elements  of  one  range  cell.  The  transformation  is  in  the  orthogonal 
direction.  The  difference  arises  in  Format  II  where  both  direction  are  transformed.  In  this 
format,  one  should  use  the  number  of  independent  samples  per  range  cell  in  defining  the  re¬ 
sultant  amplitude  and  phase  error  for  the  cell.  By  independent  samples,  we  mean  those 
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chosen  far  enough  apart  so  that  the  random  sweep  errors  are  not  correlated.  There  Is  corr¬ 
elation  between  successive  positions  of  the  actual  9-bit  data  words  so  that  using  thcce  data 
points  as  independent  (in  the  sense  of  positional  errors)  gives  an  optimistic  estimate  of  the 
performance  in  Format  II.  Using  the  range  cells  as  the  basic  Independent  elements  Is  more 
pessimistic,  since  there  are  fewer  samples  superimposing  in  the  transform  pattern. 

b.  Format  I  Frrors 


Random  errors  will  be  considered  first,  then  systematic  errors  which  are  of  lesser 
importance.  In  this  format,  a  range  cell  swath  along  the  film  (x-dlrectlon)  is  being  trans¬ 
formed  It  acts  as  a  linear  array  of  elements,  one  for  each  radar  pulse  trace.  Amplitude 
errors  arise  from  phosphor  and  film  irregularities  and  sweep  errors  which  stretch  out  c.r 
compress  the  9 -bit  data  samples.  Phase  errors  (which  are  more  critical)  arise  from  sweep 
positional  errors  and  from  film  drive  Jitter.  The  amplitude  errors  and  the  first  type  of  phase 
errors  are  similar  to  ••clement  current  errors”  In  an  antenna  array  while  the  second  type  of 
phase  error  in  similar  to  "element  placement  error"  in  the  array.  The  former  affect  all 
beams  (or  positions  along  the  output  doppler  axis)  equally.  The  latter  have  no  effect  for  the 
broadside  beam  (zero  doppler)  and  maximum  effect  on  the  most  deviated  beam  (maximum 
doppler).  The  worst  case  wili  be  considered.  Because  maximum  (positive)  doppler  corres¬ 
ponds  to  data  having  a  phase  shift  of  s  radians  between  traces,  a  fractional  error  In  the 
trace  spacing  Ax  (positional  error) /(trace  spacing)  results  in  a  worst  phase  error  of 
xo 


Actual  errors  for  a  beam  corresponding  to  a  smaller  doppler  are  reduced  by  the 


ratio  of  this  doppler  to  the  maximum  doppler.  Jitter  in  the  trace  sweep  produces  positional 
errors  in  the  y  -direction  If  this  error  is  Ay.  the  associated  phase  error  is  2rfoAy  where 
fo  is  the  spatial  carrier. 

The  one -dimensional  transform  used  to  process  Format  1  data  has  the  form: 

/'(«■„,  f)  *  £  T(m-f)e  *2',f ’{  m  fixed  (13) 

(  -I. 


where  (2L»I)  :  number  of  radar  pulses  integrated 

i.  :  output  light  distribution  (first  fringe) 

f  :  transform  variable  in  output  piane,  proportional  to  doppler  frequency 
m  :  delay  coordinate  of  the  mth  ranee  swath,  measured  across  the  film. 

* (  '■  positional  coordinate  of  the<  sample  point  in  the  range  swath,  measured 
along  the  length  of  the  film. 

complex  transmission  coefficient  for  the  film,  corresponding  to  the  upper 
sideband  (fhst  fringe)  produced  by  the  transmissivity  expressed  in  Fquat- 
ion  12. 


For  a  given  echo  center,  therm-  coordinate  (delay)  is  nearly  constant.  The  integratlor 
or  summation  is  done  for  a  fixed  value  of  delay  so  the  dependence  'n  this  coordinate  will  not 
be  carried  cxpllclty.  The  notation  can  be  simplified  to: 


M1?)  = 
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where  tj  =  -  2sf%0  :  the  cross- range  variable 

the  nominal  spacing  of  the  traces  (pulses)  on  the  film 
c  -a  con>plex  constant,  exp  (2s  i  f0  m  rQ.  representing  the  nominal  phase  factor 
for  the  m'e  range  ceil  ’ 

To  .-spacing  of  the  rfmge  cells  on  the  film. 
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q.  : amplitude,  a  random  quantity  with  mean  equal  to  the  nominal  (correct)  amplt- 

*  tude  and  standard  deviation  (a  fi.  ).  Thus,  the  standard  deviation  o(  the 

fractional  amplitude  error  ls<?. 

a 

9^:  phase,  a  random  quantity  with  mean  zero  and  standard  deviation  It  is 
composed  of  the  sweep  and  trace  spacing  error  terms  discussed  above. 

Error  analysis  which  apply  to  this  problem  have  been  done  many  times-  J.  Ruse  "Phy¬ 
sical  Limitations  on  Antennas"  TR248  Res.  Lab.  of  Elect  Mil  Oct.  30,  1952  contains  an 
early  example.  We  outline  the  error  analysis  but  do  not  carry  out  ail  the  steps.  Onedeflnes 
the  real  and  imaginary  parts  of  W  (r)  ; 


L 

cos(-tx*9| ) 

(15a) 
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sin(f  r 

(15b) 
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the  envelope  Iwl  and  phase  $: 
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IWL  (/■  *y2) 
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$>«  tan  (  Y/  X  ) 


(15c) 


For  large  N  (say  greater  than  10)  and  mild  assumptions  on  the  distribution  of  aj  and  6 
can  invoke  the  central  limit  theorem  and  treat  X  and  Y  as  Gaussian  variables.  The  pr 
ility  distributions  of  X  and  Y  are 
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When  the  amplitudes  have  an  even  distribution  (&i  “  a  |)  *«d  when  the  phase  errors  are 

symmetrical  about  zero  so  that  sin  01  ■  0.  one  has  the  mean  of  the  imaginary  component 
equal  to  zero. 


Y  -  O 


Under  these  conditions,  one  can  also  show  that  the  correlation  coefficient  of  X  and  Y  van¬ 
ishes.  Since  they  are  also  nearly  Gaussian,  they  are  almost  Independent,  and  their  Joint 
probability  p  (X.  Y)  will  be  taken  as  the  product  of  p(X)  and  p(Y).  In  order  to  f.nd  the  stat¬ 
istics  of  the  envelope  iWland  phase  <E>one  converts  to  polar  coordinates. 


The  distribution  for  the  envelope  becomes  the  modified  Rayleigh  distribution 


P  (1*0 
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(16) 


2  2  2 

The  approx i mat i on <ra><ry«ot  has  been  used.  It  turns  out  that  this  is  valid  in  the  sidelobc  re¬ 
gion  of  the  integrated  pattern  W  (x)  .  Wheno»x  (noise  power  much  greater  than  sidclobe 
power),  the  distribution  is  nearly  Rayleigh: 


P (  1*1) 
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When  <r«X  (noise  power  much  less  than  nominal  sldelobe  power),  the  distribution  ha*  s 
Gauss»an  shape  about  the  mean  sldelobe  pattern. 


P ( I w | ) .  f  |W| 
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(  (wi/x2 ) 

;<r«  X 


(18) 


The  limiting  form  for  noise  larger  than  nominal  stdelobes  will  be  used  for  estimating 
allowable  tolerances.  In  this  case,  the  mean  Is  approximately^^-  a- and  the  root  mean 
squar^/2  <r  .  There  is  a  99%  probability  of  all  values  being  less  than  the  3  <r  value. 

Tht  parameters  of  the  envelope  distribution,  X  and^,  can  be  simply  expressed  in 
terms  of  the  individual  phasor  errors.  For  simplicity,  the  phase  errois  are  considered  to 
be  normally  distributed  with  mean  zero  and  common  variance  cr2  The  percent  or  fractional 
amplitude  eirors  are  taken  to  be  normally  distributed  with  zero  mean  and  common  variance 


We  have 


.  x  o-1*2)  cos  { 


Y  -  0 


The  variances  of  X  and  Y  may  be  calculated  using  the  expressions 
2'-  -  2  ,  — „ — 2 
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The  average  o(<rx2  and  a  y2  is  independent  of  r  . 
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The  mean  square  sldelobe  power  level  (due  to  errors)  Is  2<r2.  The  ratio  of  this  mean  square 
power  level  to  the  expected  main  lobe  peak  power  is: 


--2 
X  (o) 


R.  (<To2 


~4 


4"e  1 


(20) 


The  (l-<re  )  term  in  the  denominator  arises  from  the  factor  COS  0r  ln  the  expected  peak. 
The  summation  in  the  denominator  is  the  nominal  peak  of  the  pattern.  The  factor  R  the 
ratio  of  the  two  summations  in  the  last  equation,  is  equal  to_J_  for  a  flat  amplitude  dls- 

2L*1 

tributlon  of  (2L*1)  components.  When  taper  is  employed  this  factor  Is  somewhat  greater 
than  (2 L ♦  1 )  **  ;  for  example,  it  is  0.  025  for  51  lines  with  33  db  sldelobe  suppression  A 
curve  of  R  Is  given  in  Appendix  II. 


The  side  power  has  37%,  probability  of  exceeding  the  mean  square  level  However 

Lh7!u‘8  w  ,y  [%  probab“uy  that  a  sldelobe  exceeds  the  mean  square  level  by  more  thaii 
6.  s  db.  If  we  set  the  mean  square  level  at  30  db,  the  *9%  confidence  level  is  23.  5  do  and 
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the  parameter  2cr  2/  x  2  (0)  Is  0.  001.  The  number  of  pulses  Integrated  Is  not  entirely  our 
choice  since  It  depends  on  the  rotation  rate  erf  the  object.  We  will  take  the  conservative 
estimate  of  100  and  make  R  =  0.  01.  Taking  the  major  error  source  to  be  In  the  phase  (trace 
and  carrier  positioning),  one  has 

2 

a  0.001  or  ^  Ko.  3  radian  =  1  cycle. 

'  o.  0i  9  20 

0 

When  the  phase  error  is  dominated  by  the  sweep  velocity  jitter,  the  tolerance  erf  1/20 
cycle  results  In  the  following  positioning  tolerances. 


Spatial  Carrier 

Film  Position  Error 

|  Fractional  Error  (50  ran  trace) 

15  line  pr/nm 

(1/300)  mm*  3.  3  micron 

6.7  x  10  *5 

10  line  pr/  nm 

(l/200)nm  =  5  micron 

10  *4 

(100  pulses  Integrated) 


If  the  phase  error  Is  attributed  solely  to  film  drive  non-llnearltles,  the  1/20  cycle 
tolerance  results  In  a  fractional  trace  position  error  of  10°’ >.  The  film  drive  Is  much  better 
than  this,  sc  this  error  source  may  be  neglected  In  practice.  To  be  comparable  to  the  phase 
error  of  »./20  cycle,  the  amplitude  error  would  be  about  30%.  Thus,  the  positional  jitter  due 
to  sweep  velocity  will  dominate  the  error  sour'-cs. 

In  Figure  8,  the  expected  sldelobe  ievel  Is  plotted  versus  the  design  level  (which  Is 
achieved  by  a  certain  amplitude  taper  across  the  aperture)  for  various  values  of  the  param-t- 
er  «r./  X  (0).  These  curves  show  a  saturation  level  when  the  error  power  dominates  the  de¬ 
sign  sldelobes. 

In  the  region  beyond  the  knee  In  the  curves  (where  the  error  power  dominates  the  de¬ 
sign  level),  the  rms  sldelobe  is  1  05  db  larger  than  the  expected  sldelobe  and  the  99%  con¬ 
fidence  level  is  7.  58  db  larger.  For  example,  If  the  expected  peak  to  sldelobe  ratio  Is  30  db, 
the  peak  to  rms  sldelobe  Is  29  db,  and  the  peak  to  99%  level  Is  22  db 

Format  1  Involves  a  one-dtmenslonal  integration  In  the  direction  of  film  motion  (x- 
dlrectlon)  with  the  spatial  carrier  in  the  y-dlrectlon.  The  envelope  data  In  the  y-directlon 
is  simply  imaged  onto  the  output  film.  We  have  seen  that  random  positional  errors  along  the 
sweep  should  not  exceed  one  part  on  104  Systematic  errors,  which  repeat  Identically  on  each 
trace  can  be  much  larger  In  order  to  estimate  peak  locations  to  1  10’th  of  a  range  cell,  the 
maximum  systematic  deviation  In  the  range  direction  should  not  exceed  one  part  In  8(K> 
Systematic  variations  In  the  film  velocity  have  the  same  effect  in  both  formats, 
and  the  analysis  will  l>c  outlined  Ir.  the  following  discussion  of  >:ormat  II 


c  Format  II  Errors 


Format  II  Involves  a  two-dimensional  Integration  for  each  output  cell.  Ruse’s  paper 
shows  that  In  this  case  the  factor  R  becomes 


M  2 

l 

m*-M/2 


M/2  J  \  2 

X  m>*  I  1 

m  =  -M/'2  I  M(2(L+1)) 


(21) 


As  was  pointed  out  in  the  Introduction  to  this  section,  the  number  M  should  be  the  numb®r 
of  statistically  independent  samples  along  the  data  trace.  A  pessimistic  estimate  Is  that  an 
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Figure  8.  Sidolobe  Level  Vs  Design  Level 


entire  range  cell  (roughly  10  cycles  of  the  spatial  carrier)  Is  correlated.  The  number  of 
Independent  samples  In  a  trace  will  be  taken  to  be  100  (which  Is  an  order  of  magnitude  smal- 
icr  than  the  number  of  samples  per  sweep. 

4  Uslnf>  100  pulses  as  before,  with  100 independent  sweep  samples,  one  finds  R  to  be 
0  Instead  of  10-  .  Thus,  the  random  error  tolerance,  o* ,  can  Increase  by  a  factor  of 
10  compared  to  Format  l  For  given  tolerances,  the  added  dimension  of  Integration  will 
raise  the  performance  limit  by  20  db.  This  can  be  seen  in  Figure  8  .  It  is  therefore  ex¬ 
pected  that  random  errors  wi.l  degrade  the  output  when  Format  II  is  used  as  severely  as 
in  Format  1.  ! 


The  systematic  errors  in  Format  II  can  create  a  problem.  The  repeatable  non-llnear- 
ities  in  the  sweep  are  not  simply  imaged  on  the  output  film.  Instead,  the  processor  does 

a  spectral  analysis  o i  ‘he  non-linearity  and  produces  mam  lobe  broadening  and  localized 
sidelobes  accordingly. 

One  divides  systematic  non-linearities  into  two  types  -  those  without  a  repetitive 

structure  and  those  with  repetition  or  oscillation.  The  first  are  analyzed  by  representation 

with  a  few  terms  of  a  power  series;  the  second  are  represented  by  a  Fourier  series  The 

power  series  type  produces  distortion  of  the  main  lobe  and  perhaps  close-in  sidelobes  (for 

*m*11  *rjS}  Th°  °Scl,,atory  ,yfH'  Klvcs  r*8<>  to  localized  sidelobes  having  the  appearance 
of  diminished  main  lobes  or  '  echoes”. 

Calculations  for  quadratic  phase  errors  appear  in  "Theory  and  Design  of  Chirp  Radars* 
by  Klauder  Price.  Darlington  and  Albersheim.  Reference  4.  A  combination  of  quadratic 

Ta  j! ‘  h.1^50  °ar0rS  ‘S  ,rt>a,<,d  in  R<‘,0rCnCe  2*  Whp"  spectrum  taper  is  employed 
appreciable  degradation  occurs  for  a  maximum  phase  difference  across  th-  aperture  of 

1  „  For  weighted  spectra  (say  25-30  di, nominal  sidelobes)  total  phase  deviation  of  » 

*2“ 

tov  is  toIerabie  Actually,  systematic  error  of  deviation*  will  produce  large  main  lobe 
broadening  (a  factor  of  1.  4  to  1.  5)  and  peak  reduction  of  15%.  A  deviation  of  Is  much  less 

2 

severe,  5%  reduction  of  peak  and  broadening  by  a  factor  of  1.  1. 

At  a  spatial  carrier  of  10  line  pr  ntn,  a  deviation  of  n  represents  1  '20  mm  This  r  - 
presents  (1  1000)  «h  of  the  sweep  length  and  is  a  more  stringent  than  in  Format  I.  In  fact, 
Format  II  has  comparable  tolerances  on  random  and  systematic  errors,  about  1  part  in  1000. 


Errors  of  the  oscillatory  type  are  best  represented  by  Fourier  components 
ample,  an  amplitude  ripple,  would  modlf\  the  pattern  as  follows* 

M  2 


For  ex- 


W  (T  ) 


I 
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-M  2 


ni 


[l  *A  cos  (niT*  a)J  e 


im  r 


(22) 


la 


=  wo  (  r  )  ♦  2~  Ae  WQ  ( r  )  .  i-  A  e 


-  t  a 


W0(r  -T  ) 


here  W0  Is  the  nominal  pattern  and  the  summation  represents  integration  in  the  dimention 
having  the  oscll  atory  error.  Thus,  in  addition  to  the  nominal  pattern, there  are  two  "paired 
echoes  of  relative  amplitude(l  2)A,  located  *T  away  from  the  central  peak.  If  we  require 
no  pa  red  echo  (due  to  amplitude  ripple)  be  greater  than  -26  db,  A  cannot  exceed  0  1  or 
10%  ripple  amplitude. 
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A  similar  analysis  applies  to  repetitive  phase  errors.  A  single  phase  oscillation  of  the 
form  B  cos  (mT  +/7)  will  be  considered. 

M/2  LBcos(mT+/3)  im  r 

W(  r)  =  z  Q  m  e  e 

m  =  -  M/2 


M/2 

I 

m-  -  M/2 


0 m[ J0  ♦2l Jj  (B)cos(mT ♦  )-2J2(B)  cos(2mT+2£)  4  ...  e  lmT 

(23) 


For  small  phase  ripple,  B  <  0.  4  radians,  one  uses  the  approximations 


VI.  J 1  (0)  =  1/2  B,  J2  (B)  =  0  etc.  Tne  expression  above  reduces  to 

i/9  1  a 

W  (r)  =  W0  (r)  4  1/2  iBe  W0  (r  +T )  ♦  1/2  i  Be  ~  W0  (r  -  T)  (24) 

There  are  two  "paired  echoes"  arising  from  t.*e  phase  ripple,  In  quadrature  with  those 
associated  with  a  similar  amplitude  ripple.  If  the  criterion  of  -26  db  suppression  of  these 
ecnoes  is  .applied,  one  must  keep  tne  ripple  coefficient  B  less  than  0.  1  radian  (5.7°). 
Suppression  of  -20  db  would  require  the  phase  ripple  less  than  0.  2  radian  (0.  032  cycles). 


The  tolerance  on  the  oscillatory  part  of  the  systematic  phase  error  is  quite  severe. 

At  a  carrier  of  10  line  pr/mm  the  0.  1  radian  tolerance  (-26  db  suppression)  represents  only 
(1/630)  mm  on  the  film,  when  in  the  sweep  length.  It  should  be  emphasized  that  these  effects 
are  associated  with  a  sinusoidal  sweep  or  film  velocity  variation  giving  range  lobes  and 
cross-range  lobes  respectively,  ft  does  not  represent  the  total  systematic  error  which  un- 
doubtly  has  a  non  oscillatory  part.  Degradation  caused  by  systematic  errors  in  »he  sweep 
velocity  have  been  observed  in  processed  outputs  with  Format  II  recordings.  Slmusoidal 
drive  speed  variations  would  produce  paired  echoes  in  cross-range  for  both  formats.  This 
has  not  been  observed  in  the  outputs. 


5.  REFERENCES  FOR  SECTION  11 
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SECTION  ni 


EQUIPMENT  DESIGN  PARAMETERS 


The  film  recorder  was  designed  to  provide  a  stable  and  accurate  trace  from  which  to  re¬ 
cord  data  on  film.  The  recording  on  film  is  accomplished  by  repeating  a  one-line  scan  across 
the  face  of  a  t  RT  and  moving  the  film  past  the  image  plane  at  a  constant  rate.  The  rate  of 
film  motion  is  set  such  that  lines  are  written  closely  spaced  with  some  overlapping, 

1.  CATHODE  RAY  TUBE  ASSEMBLY 

The  resolution  requirement  was  to  provide  a  dynamically  focused  spot  of  less  than  1.5 
mils  size  (measured  at  the  70%  points).  Sweep  accuracy  was  to  limit  spot  jitter  to  0.  5  mils 
from  sweep  to  sweep.  Intensity  compensation  was  provided  to  prevent  long  term  threshold 
Instability  of  the  CRT  trace. 


The  CRT  Assembly  contains  all  circuitry  necessary  to  display  video  information  by  inten¬ 
sity  modulating  a  single  line  trace  on  a  5CEPI1  cathode  ray  tube.  External  signals  necessary 
for  proper  operation  of  the  film  recorder  include  a  start  pulse,  stop  pulse,  film  drive  fre¬ 
quency  control,  and  video  input. 

The  start  and  stop  pulses  provide  the  proper  timing  and  voltage  amplitudes  required  to 
sweep  a  four  Inch  trace  on  the  face  of  the  CRT.  The  start  pulse  Is  reshaped  by  a  logic  differ¬ 
ential  receiver  and  expanded,  via  a  monostable  multivibrator.  Reshaping  and  expandli*  of 
the  stop  pulse  is  also  provided.  Both  pulses  provide  Inputs  to  two  nand  gates  connected  in  a 
Rip  flop  configuration.  The  output  of  this  nip  nop  determines  the  time  that  an  electronic 
switch  allows  the  integrating  capacitor  to  charge.  This  output  is  then  amplified  and  offset 
about  zero  volts  since  the  final  deflection  amplifier  Is  double  ended  (push-pull).  To  compen¬ 
sate  for  the  fact  that  the  5CEP  is  a  flat  faced  tube,  the  sweep  voltage  is  predistorted  prior  to 
final  amplification.  This  predistortio  i  is  accomplished  thru  a  diode -resistor  linearity  matrix 
which  changes  the  shape  of  the  sweep  lamp  such  that  the  trace  velocity  will  be  uniform  across 
the  face  of  the  tube. 

Variation  in  line  width  or  spot  growth  from  the  center  of  the  tube  to  the  edge  can  be  as 
great  as  5:1  over  the  center  specified  spot  size  (1.  5  mils  in  5CEP).  To  prevent  this  growth, 
dynamic  focusing  must  be  accomplished.  The  ideal  dynamic  voltage  waveshape  is  a  120  volt 
parabola  which  modulates  the  focus  power  supply.  To  generate  this  waveshape,  a  piece-wise 
linear  approximation  to  the  parabola  from  an  operational  amplifier  modulates  the  focus  power 
supply.  A  brightness  control  is  provided  so  that  proper  operating  bias  may  be  obtained  from 
tube  to  tube. 

To  insure  that  the  design  goals  are  set  forth  for  the  film  recorder  were  met,  the  com¬ 
mercial  power  supplies  were  specified  to  have  0.  05%  regulation.  The  anode  and  focus  power 
supplies  were  built  in-house  by  Westinghouse  and  were  0.01%  regulation.  This  tight  regula¬ 
tion  was  to  prevent  random  spot  growth.  An  attempt  was  made  In  the  design  of  the  analog 
circuits  to  molecularlze  as  much  circuitry  as  possible. 

2.  CAMERA  SYSTEM 

a.  Lens 

The  system  camera  was  provided  by  Photogrammetry,  Inc. ,  of  Rockville,  Maryland. 

It  was  designed  specifically  for  this  application  and  is  designed  to  record  intensity  modulated 
traces  from  a  5CEP11  CRT  onto  continuous  70  mm  film  at  a  2:1  reduction.  The  film  Is  con¬ 
tinuously  advanced  by  a  frequency  synchronized  capstan  during  exposure.  There  Is  no  shutter 
and  exposure  Is  controlled  by  sweep  intensity  and  lens  opening. 


27 


The  lens  Is  a  Boyer -Beryl  110  mm  f.  L.  Beryl  lens  having  a  maximum  aperture  if 
f/6.  8.  The  lens  r.ounting  allows  continuous  focus  adjustment  over  the  desired  range.  A  re¬ 
flex  mirror  is  so  mounted  that  it  may  be  rotated  to  reflect  light  collected  by  the  lens  on'o  a 
high  resolution,  Polacoat  Lenscreen.  which  Is  set  closely  in  the  focal  plane.  The  CRT  sweep 
was  focused  cn  this  screen  and  observed  with  a  120X  microscope  assuring  optimum  practical 
focus. 

b.  Film  Drive 

The  camera  drive  system,  manufactured  by  Sequential  Electronics  Sy  ms,  Inc. , 
consists  of  a  film  drive  servo  controlled  capstan  and  a  constant  torque  film  ten? 'on  take-up 
system.  This  may  be  s  en  in  Figure  9  .  The  capstan  servo  system,  completely  eliminates 
the  inherent  limiting  characteristics  of  sampled -data  phase  locking  systems  that  derive  feed¬ 
back  information  from  a  pulse  rate  tachometer. 

The  infinite  resolution  Band-Scan  Riadout  has  the  unit^ie  capability  of  ele'  tro-optlcally 
monitoring  continuous  speed/phase  coordinates  converting  this  information  into  the  required 
digitized  format  at  a  constant  conversion  rate  which  is  independent  of  shaft  speed.  It  com¬ 
putes  speed  with  zero  error  and  position  or  phase  with  a  maximum  peak  inaccuracy  of  0.33 
arc  seconds,  at  a  minimum  settable  conversion  rate  of  20  KHz.  This  ac'lon  produces  a  rate/ 
(Kisltion  control  system  with  a  constant  closed-loop  control  bandwidth  in  the  order  of  3  KHz. 

The  reference  input  to  the  Band -Scan  Reference  Generator  is  a  olgnal  whose  frequency 
directly  corresponds  to  commanded  motor  speed  and  whose  phase  relates  to  drive  shaft  posi¬ 
tion.  The  standard  speed -frequency  conversion  ratio  is  such  that  Motor  RPM  >  K  *  f,  where 
K  is  the  conversion  ratio  selected  for  a  given  application  and  f  is  the  reference  frequency. 

The  Reference  Generator  samples  and  stores  speed/ phase  command  Information.  Updating 
of  the  coordinates  in  storage  is  once  per  cycle  of  the  reference.  The  interrogation  rate  of 
the  Generator  for  information  fed  into  the  Band-Scan  FPL  is  automatically  made  compatible 
with  the  conversion  rate  of  the  Band -Scan  Readout. 

The  Band-Scan  Computer  compares  the  commanded  speedup,  ise  Information  produced 
by  tiie  Reference  Generator,  with  the  motor  shaft  speed  phase  coordinates  measured  by  the 
Band-Scan  Readout,  and  produces  a  digitalized  error  signal.  The  Digital  to  Analog  Converter 
converts  this  error  signal  into  an  analog  voltage,  with  a  conversion  rate  that  is  automatically 
made  compatible  with  the  conversion  rate  of  the  Band-Scan  Readout.  This  conversion  is  per¬ 
formed  with  no  significant  time  constant,  and  with  virtually  no  ripple,  which  permits  design 
of  an  extremely  wideband,  high  gain  control  loop.  The  Stabilization  Section  contains  stabili¬ 
zation  circuitry  which  applies  phase,  velocity,  and  acceleration  control  ir.  the  proper  ratio, 
and  with  the  proper  frequency  characteristics,  to  optimize  the  desired  closed  loop  performance 
of  the  control  loop.  The  Driver  unit  supplies  the  correct  excitations  to  the  motor  windings 
to  maintain  constant  shaft  velocity.  The  driver  and  motor  windings  are  operated  in  a  secon¬ 
dary,  closed -loop  manner,  to  linearize  the  torque  output  characteristics,  and  to  virtually 
eliminate  the  electrical  time  constant  of  the  motor  winding. 

The  position  sensitivity  of  the  Readout  is  such  that  if  the  electronic  gain  were  unity, 
less  than  100  arc  seconds  of  position  error  would  correspond  to  a  full  applied  torque  command. 
Achievable  electronic  gains  with  the  proper  stabilizing  element1-  are  on  the  order  of  40  db. 

This  results  in  a  total  open-loop  gain  of  1,  000,  000  (120  db)  i adlaus/sec/radlans;  it  produces 
an  absolule  closed  loop  position  control  accuracy  on  the  order  of  seconds  of  arc. 

3.  DATA  RATE 

The  input  to  the  SDS  910  is  a  tape  written  at  a  density  of  556  bits  per  inch.  On  this  tape 
are  records  of  date  (a  record  Is  one  pulse  of  radar  data  which  is  also  a  single  sweep  across 
the  face  of  the  CRT)  each  of  which  is  1024  or  2048  words  in  length.  There  are  nine  bits  per 
word.  The  highest  possible  input  rale  to  the  SDS  910  Is  5  records  per  second,  the  lowest  Is 
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2.  5  rerords  per  second-  This  corresponds  to  200  and  400  milliseconds  between  line  start 
pulses,  respectively.  The  highest  data  rate  at  the  output  of  the  SDS  910,  therefore,  is 
92,  160  bits  per  second  and  it  occurs  when  using  records  of  2048  words  and  200  milliseconds 
between  lines  star,  pulses  (i.e.  92,  160  bits/sec  =  5  X  9  X  2048).  The  lowest  usable  data  rate 
is  23,  040  bits  per  second  =  2.  5  X  9  X  1024;  this  occurs  when  using  records  of  1024  words  and 
400  milliseconds  between  line  start  pulses.  Normally,  the  film  recorder  is  operated  using 
the  highest  data  rate. 

It  is  necessary  to  know  how  many  pulses  of  radar  data  can  be  written  on  a  magnetic  tape. 
Since  one  word  requires  12  bits,  at  a  tape  density  of  556  bits  per  inch,  there  are  556/2  words 
per  inch.  If  the  record  length  is  2048  words,  then  have  (2/556)  2048  *  7.  37  inches/pulse. 
Between  each  pulse  there  is  a  .  75  inch  record  gap.  Therefore,  the  number  of  inches  per 
pulse  is  7.  37  ♦  .  75  =  8.  12.  The  magnetic  tapes  wc  use  are  2400  feet  long.  Therefore,  the 
maximum  number  of  pulses  per  tape  Is  (2400/8.  12)  12  =  3550.  At  a  record  length  of  1024 
words,  there  are  2(3550)  or  7100  pulses-  tape. 

4.  D/A  CONVERTER 

The  D/A  converter  (SDS  DX15,  SX11)  is  used  to  convert  the  digital  data  from  the  SDS  910 
to  an  analog  voltage  that  is  applied  to  the  grid  of  the  CRT  of  the  film  recorder.  The  block 
diagram  below  (Figure  10  )  shows  how  the  D/A  converter  is  used  in  conjunction  with  the 
SDS  910  and  the  film  recorder  (all  numtx>rs  are  SDS  numbers). 
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Figure  10  .  Block  Diagram  of  SDS  and  D/A  Converter 
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The  output  of  the  SDS  910  are  records  of  either  1024  or  2048  words  each.  Each  word  has 
an  amplitude  in  the  range  -255  to  256.  The  D/A  converts  these  digital  amplitudes  to  analog 
voltages  fanging  from  -10  to  +10  volts.  The  actual  signal  supplied  to  theCRT  is  shown  in  the 
diagram  below  (figure  11  )  which  assumes  2048  words  per  line. 
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Figure  11  .  DA  Converter 

Thus  each  computer  word  (amplitude)  is  replaced  by  a  constant  voltage  of  32  u  tec.  The 
higher  the  voltage  supplied  to  the  film  recorder,  the  brighter  the  spot  will  be  on  th>  face  of 
the  CRT. 

In  addition,  the  following  line  start  and  end  of  line  pulses  are  supplied  to  the  film  .  e- 
corder:  (See  figure  12). 


Start  of  Line  Pulse 


1 


+8  volts 
0  volts 


End  of  Line  Pulse 


8  volts 
0  volts 
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Figure  12  .  Line  Start  and  Stop  Pulses 


5.  FILM  DEVELOPMENT 

Kodak  2490  Film  may  be  processed  and  loaded  in  the  camera  in  the  presence  of  red  light 
with  wavelengths  longer  titan  5750A°. 

A  special  light  source  has  been  assembled  which  provides  variable  illumination  without 
affecting  flint  fog  level. 

Film  development  is  done  in  a  special  dark  tank  which  has  a  spool  for  70mm  film.  Chemi¬ 
cals  are  poured  into  and  front  vents  in  the  top  at  the  cannister. 

Steps  in  the  development  are  listed  in  sequence. 

(1)  After  film  is  ioaded  on  S|xtol  and  placed  in  cannister  with  cover  secured,  fiil  with 
developer  and  continuously  agitate.  Time  in  developer  is  taken  front  gamma  vs.  time 
curve.  (Figure  13  ). 

After  elapsed  time  the  developer  is  |>oured  out.  If  ciose  time  control  is  desired  it  will 
be  necessary  to  remove  the  cannister  cover  in  red  light  and  pour  Ihe  solution  out 
rapidiy  and  proceed  immediately  to  Step  2. 

(2)  Fill  cannister  with  KODAK  stop-bath ,  pour  out  after  approximately  25  seconds. 
Because  of  the  long  time  required  to  fill  and  empty  the  cannister  with  the  top  on  it  is 
recommended  that  this  be  done  with  top  off  in  red  light.  The  purpose  of  litis  step  is 
to  stop  the  action  of  the  developer.  Instead  of  KODAK  Stop  Bath  the  cannister  may  i>e 
flushed  with  water  for  approximately  10  seconds  or  longer  and  poured  out. 

(3)  Fill  cannister  with  FIXER  Solution.  The  time  for  KODAK  Fixer  or  KODAK  Fixing 
Bath  F-5  is  2-4  minutes.  For  KODAK  Rapid  Fixer  allow  2-3  minutes.  Pour  out.  The 
film  is  now  no  longer  sensitive  to  room  light. 

(4)  Flush  cannister  with  65-95°F  water  for  approximately  30  seconds  and  pour  out. 

(5)  Fill  cannister  with  Hypo-cleaning  Agent  for  1-2  minutes.  Pour  out. 

(6)  Wash  film  in  running  water  for  2-4  minutes. 

(7)  (Optional)  To  avoid  water  spots  and  streaks  on  the  film  the  negative  can  be  washed 
in  photo -flo  solution  before  drying. 
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Figure  13  .  Gamma  vs  Development  Time 


G.  COMPUTER  SOFTWARE 

There  are  two  main  programs  that  are  used  to  process  calibrated  radar  data  and  put  it 
In  a  form  suitable  for  use  by  the  film  recorder  and  the  optical  processor.  The  first,  the 
Optical  Processor  Program,  is  used  on  t  tie  UNTV  AC  1108  computer.  It  accepts  pulses  of 
radar  date,  interpolates  each  pulse  to  produce  2048  or  1024  points,  applies  the  spatial  car¬ 
rier  with  phase,  and  adds  a  bias  level.  This  data  Is  then  written,  a  pulse  at  a  time,  on  a 
tape  which  is  in  a  format  compatible  with  the  SDS  910  computer.  The  second  program  is 
used  by  the  SDS  910  computer.  It  accepts  the  tape  produced  by  the  Optical  Processor  Program, 
and  alternately  reads  from  the  tape  pulses  of  radar  data  (each  pulse  Is  a  record  of  2048  or 
1024  points)  and  ou’puts  them  at  a  fixed  rate  to  the  D  A  converter  and  then  to  the  CRT.  In 
addition,  it  also  pr.  vides  line  start  and  line  stop  pulses  and  camera  clock  frequencies  to  the 
film  recorder. 

In  addition  there  is  a  third  program,  DATA2X,  which  simulates  radar  data  and  outputs  it 
onto  tape  in  the  calibrated  radar  data  format.  This  program  is  used  to  test  the  UNIVAC  1108 
software,  the  SDS  910  software,  and  the  optical  recording  system  when  real  radar  data  Is  not 
readily  available.  The  program  can  simulate  radar  returns  from  a  rotating  target  of  rigidly 
connected  scatterers.  At  the  present  time  it  Is  capable  of  simulating  up  to  twenty  scatterers 
each  having  an  arbitrary  amplitude,  relative  phase,  and  any  position  relative  to  center  of 
rotation.  These  can  be  rotating  at  a  described  rotation  rate,  and  recorded  at  a  given  pulse 
repetition  frequency,  for  any  number  of  pulses.  This  program  also  retains  the  maximum 
amplitude  for  any  number  of  pulses.  This  maximum  amplitude  is  used  as  a  bias  level  for  a 
frame  of  pulses  by  the  Optical  Processor  Program. 

More  detail  concerning  all  these  programs  can  be  found  In  Appendix  III. 
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SECTION  IV 


RESULTS  OF  SUBSYSTEM  TEST 


1.  FILM  RECORDER 

The  film  recorder  coT.sole  receives  a  nalog  signal  information  which  Intensity  modulates 
the  beam  of  a  CRT.  This  display  is  recorded  by  a  continuous  70  mm  servo-controlled  moving 
film  camera. 

The  display  console  receives  from  its  signal  source  three  signals,  a  pulse  to  start  the 
CRT  sweep,  analog  video  information  to  be  displayed  during  the  sweep,  and  a  sweep  stop 
pulse.  The  time  between  sweep  start  pulses  may  be  varied  at  '“ill  between  limits  but  the 
video  display  time  is  fixed  at  approximately  68  msec. 

A  simplified  block  diagram  of  the  display  console  is  illustrated  in  figure  14  . 

a.  Circuit  Operation 

Upon  receiving  the  sweep  start  pulse  the  logic  circuitry  generates  a  variety  of  pulses 
at  different  time  positions  to  operate  the  sweep,  unblanking  and  optical  feedback  circuits. 

These  pulse  waveforms  are  shown  on  the  system  block  dir»'  am  with  the  times  referenced 
against  the  externally  provided  sweep  start  pulse  at  t  .  .0  msec  rectangular  pulse  causes 

a  gate  switch  to  unclamp  a  Miller  integra'or  which  th?n  generates  a  linear  ramp  of  70  msec 
duration.  This  amplified  ramp  drives  the  CRT  horizontal  deflection  coils  providing  the 
horizontal  sweep  trace.  The  logic  circuitry  provides  a  2  msec  pulse  which  provides  switching 
pulses  to  operate  the  optical  feedback  loop.  This  loop  serves  to  maintain  a  constant  operating 
brightness  of  the  CRT  tracer  by  providing  a  negative  feedback  loop  which  has  as  its  reference 
an  input  voltage  set  by  the  "calibrate  intensity"  control.  This  serves  to  maintain  a  constant 
brightness  as  CRT  parameters  such  as  cathode  emission  and  phosphor  efficiency  are  degraded 
with  long  term  aging  effects.  At  the  beginning  of  the  CRT  trace  a  phototransistor  views  the 
first  2  msec  of  the  sweep.  This  2  msec  sample  is  then  held  by  a  "hold"  circuit  for  the  dura¬ 
tion  of  the  sweep.  The  output  of  this  sample  and  hold  circuit  provides  the  feedback  voltage 
which  is  compared  against  the  "calibrate  Intensity"  input  by  their  mutual  summing  into  the 
input  of  the  video  amplifier  as  shown.  The  2  msec  display  of  the  sweep  viewed  by  the  photo - 
transistor  is  turned  on  by  the  unblanking  amplifier  which  is  driven  by  an  output  of  the  logic 
circuitry  shown  in  the  figure.  The  Brightness  control  varies  the  amplitude  of  the  unblanking 
pulse  but  this  has  little  effect  on  the  trace  brightness  because  of  the  intensity  stabilization  of 
the  optical  feedback  and  serves  only  to  establish  loop  operating  parameters. 

The  video  input  is  amplified  and  applied  to  the  CRT  control  grid.  The  dynamic  range 
of  the  CRT  intensity  modulation  is  controlled  by  the  video  gain  control. 

Dynamic  focussing  is  applied  to  the  CRT  to  compensate  for  spot  nonuniformity  across 
the  trace.  Due  to  the  geometry  of  the  tube  the  focal  length  of  the  electrostatic  lens  must  be 
varied  as  the  electron  beam  is  swept  across  the  tube  face  to  maintain  optimum  focus.  Th*s 
focus  correction  is  a  120  volt  parabolic  pulse  which  is  applied  to  the  CRT  in  series  with  the 

2.  2  kv  focus  potential. 

b.  Transfer  Curves 

Curves  are  provided  which  enable  the  operator  to  determine  the  recorded  film  ,:barac- 
terlstics.  Figure  15  ,  percent  transmission  versus  calibrate  Intensity  setiing,  g'vei'  a 
plot  of  'O'  level  operating  points  for  various  control  settings.  This  curve  is  valid  for  2490  film 
processed  for  unity  gamma,  and  a  pulse  density  of  16.  4  line  pairs/mm  on  '.he  film. 
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Figure  15  .  Calibrated  Intensity  Setting 


Figure  16  gives  the  dynamic  swing  of  the  video  Input  signal  in  terms  of  calibrate  intensity 
setting  versus  gain  setting  and  may  be  read  graphically.  It  should  be  pointed  out  that  the 
operating  conditions  are  subject  to  other  effei  ts  such  as  CRT  and  film  grain  noise,  resolution 
limits  and  line  density  on  the  film.  At  present  the  CRT  spot  has  an  astigmatic  focus  with  ap¬ 
proximately  a  2: 1  height  to  width  ratio.  With  a  600  cps  camera  Input  frequency  or  a  line 
density  of  16.4  lines/ mm  there  is  line  to  line  overlap  which  yield  le3S  than  a  5%  Increase  in 
light  transmission  between  pulses  on  the  film  negative.  This  should  effectively  suppress  the 
Raster  effect  and  relieve  the  impression  that  there  is  insufficient  exposure  of  the  film.  The 
line  density  recorded  on  the  film  is  given  by 


T  *  ‘c 

where  x  is  the  period  between  line  start  pulses  in  seconds  and  f  ts  the  camera  input  fre¬ 
quency  in  cps.  c 

c.  Spatial  Resolution 

There  ti  some  question  as  to  the  determination  of  the  system  spatial  resolution  along 
the  sweep.  Measurements  indicate  that  the  CRT  spoi  size  Is  approximately  1.  3  mils  and  is 
within  the  expected  limits.  There  ts  spot  growth  over  the  dynamic  range  of  the  system  and 
thus  the  resolution  limit  imposed  by  the  CRT  spoi  size  is  a  function  of  the  instantaneous  ampli¬ 
tude  of  the  input  video  signal  and  the  *0'  level  setting  of  the  calibrate  intensity  control.  With 
a  modulation  transfer  of  0.  5  the  resolution  on  the  film  may  be  said  to  lie  approximately  55 

lines  mm  due  to  spot  size  alone. 

\ 


Figure  16  .  Dynamic  Swing  of  the  Video  Input  Signal 
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The  sample  time  of  the  D/A  converter  yield  4.  9  X  10  mm/sample  on  the  CRT  face 
or  2.  45  X  10*2  mm/sample  on  the  film  due  to  2:1  reduction  by  camera  op.ics.  li  ts  not 
meaningful  to  speak  of  line  pr.  resolution  as  limited  directly  by  the  spatial  sample  length  but 
to  find  how  many  samples  are  required  per  unit  length  to  yield  the  desir*  i  results  in  the 
optical  processor  output  plane. 

The  combined  f/6.  8  earner*  lens  and  2490  film  with  unity  gamma  has  a  modulation 
transfer  function  of  0.  5  at  approximately  20  llnes.-4'nm  with  the  modulation  transfer  being 
defined  as 

I  max  -  Imin 
Imax  +  Imin 

where  I  is  the  transmitted  light  intensity  through  the  film  negative.  The  system  spatial  re¬ 
solution  is  a  function  of  all  the  above  mentioned  items. 

d.  Sweep  Repeatability  and  Linearity 

It  is  necessary  to  provide  close  control  of  the  beam  as  it  is  swept  across  the  CRT 
face.  In  order  to  minimize  phase  jitter  of  the  special  carrier,  special  care  must  be  taken 
to  reduce  the  electrical  noise  level  in  the  deflection  circuits.  After  evidence  had  shown  that 
such  disturbances  in  the  film  record  could  not  be  tolerated,  special  efforts  were  made  to  re¬ 
duce  circuit  noise  levels  until  there  was  no  longer  any  observable  phase  disturbances  recorded 
on  the  film.  Acceptable  noise  levels  would  produce  jitter  in  the  trace  which  is  on  thi  order 
of  one  spatial  sample  of  the  D/A  converter  or  less.  This  requires  spatial  noise  of  5  X  10*2mm 
or  less.  It  was  difficult  to  translate  this  to  circuit  noise  levels  because  the  voltage  observed 
is  related  to  the  deflection  coil  current  by 


and.  thus  the  duration  of  a  noise  pulse  is  as  important  as  its  magnitude. 

Effective  noise  that  appears  in  the  video  chain  has  been  reduced  below  the  gross  fog 
level  in  the  recorded  film.  Spurious  noise,  and  ringing  transients  that  are  of  short  duration 
arc  effectively  suppressed  both  by  electrical  filtering  and  by  the  effective  low  pass  effect  of 
the  speed  of  the  film  emulsion. 
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SECTION  V 


RESULTS  OF  SIMULATED  TARGET  TFS1S 


1.  BACKGROUND 

In  addition  to  sub- system  and  system  tests  which  were  performed  to  check  out  the  equip¬ 
ment,  a  number  of  moving  targets  were  simulated  for  subsequent  processing  in  an  optical 
correlator  to  demonstrate  system  operatior.  In  these  tests,  the  returns  from  point  scatterers 
were  simulated  as  the  scatterers  rotated  rdative  to  the  radar  line  of  sight.  These  returns 
were  simuiated  once  each  0.01°  rotation,  and  w  itten  onto  film.  The  films  were  then  sent 
to  the  Institute  of  Science  and  Technology  of  the  University  of  Michigan  where  they  were  pro¬ 
cessed  inthlfr  experimental  optical  correlator.  The  processed  film  demonstrated  the  cap¬ 
ability  of  the  optical  processor  (both  the  film  recorder  and  the  optical  correlator)  and  was 
useful  in  pointing  out  abnormalities  in  the  recording  process,  as  well  as  in  demonstrating 
successful  recording  for  processing 


2.  FILM  FRAMES  PRODUCED 


A  listing  of  the  film  frames  provided  to  1ST  is  given  in  Table  2  ,  all  of  which  used  simula¬ 
ted  data.  Further  explanation  of  the  target  configuration  and  orientation  for  which  the 
simulated  data  was  generated  follows  the  table. 

TABLE  2  FILM  FRAMES  FOR  1ST 


FRAME  NO.  TYPE  OF  FRAME 


DATA  TAPE  REQUIRED 


*1 

a  / 


SPATIAL 

CARRIER 

ONLY 


2  SCATTERERS 
> FOR MAT  1 

8 1 

9)  2  SCATTERERS 
KH  FORMAT  I 
1 U 
12 
13 

14^2  SCATTERERS 
15  FORMAT  II 

16J 

17  FChMAT  I 
Id  FORMAT  II 


10  line  pr/mni 
12  line  pr/mni 
14  line  prOum 
16  line  pr/mni 
CASE  A  FIGURE  17 
CASE  B  FIGURE  17 
CASE  C  FIGURE  17 
CASE  D  FIGURE  17 
CASE  A  FIGURE  17 
CASE  B  FIGURE  17 
CASE  C  FIGURE  17 
CASE  D  FIGURE  17 
CASE  A  FIGURE  17 
CASE  B  FIGURE  17 
CASE  C  FIGURE  17 
CASE  D  FIGURE  17 
1 1-ELEMENT  1ST  TARGET 
U-ELEMENT  1ST  TARGET 


Oi<»*  Record  or  sweep 
One  Record  or  sweep 
One  Record  or  sweep 
One  Rerord  or  sweep 
One  frame  511  pulses  5.  11°  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11°  rotation 
One  frame  511  pulses  5.  11°  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11°  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11*  rotation 
One  frame  511  pulses  5.  11*  rotation 
6  frames  continuous  5/500  pulses 
6  frames  continuous  5/500  pulses 


Frames  I  through  4  were  recordings  of  the  spatial  carrier  only  with  the  carrier  being  intro¬ 
duced  in  the  sweep  direction.  A  frame  of  about  50  mm  length  was  recorded  by  reproducing  this 
trace  repetitively.  For  these  frames  several  values  of  spatial  frequency  were  used  as  indicated 
in  the  table. 


Frames  5  through  16  were  recorded  using  a  target  having  two  point  scatterers  separated 
by  20  feet.  The  orientation  of  this  target  was  varied  as  shown  in  Figure  17  in  order  to  get 
several  types  of  recordings. 

Four  target  orientations  are  shown  in  this  figure  and  labeled  Cases  A,  B,  C  and  D.  In  all 
these,  data  was  simulated  as  the  target  routed  through  5. 11*  (or  approximately  0.  1  radlai ). 
The  orienUtions  shown  in  the  figure  are  for  the  center  of  the  recording  interval. 
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Frames  9  through  12  duplicated  frames  5  through  8.  Processing  of  frames  5  through  8 
revealed  Improper  recording  of  data  and  the  necessity  for  further  equipment  adjustments. 

After  adjustments  were  made  to  the  equipment,  frames  9  through  12  were  recorded.  These 
frames  had  data  recorded  more  densely  than  before,  having  a  frame  length  of  approximately 
20  nun  rather  than  the  original  50  mm  with  the  same  number  of  pulses  being  used. 

Frames  13  through  16  were  also  recorded  using  the  simulated  data  for  the  target  orient¬ 
ations  shown  In  Figure  17  ,  except  here  Format  II  Is  used.  In  Format  II  the  Fourier  trans¬ 
form  of  the  range  response  coefficients  is  taken  prior  to  recording.  Processing  of  these 
frames  revealed  that  there  are  excessive  range  sldelobes  appearing  In  the  correlator  output, 
which  preclude  best  processing  In  Format  II.  The  source  of  the  range  sldclobes  has  been 
traced  to  a  variation  In  the  spatial  frequency  with  distance  along  the  trace  (across  the  film). 

The  cause  of  this  problem  has  not  been  firmly  established  at  the  writing  of  this  report.  How¬ 
ever,  as  may  be  seen  In  the  processed  output  Included  later  In  this  section,  there  Is  recognis¬ 
able  Information  In  the  output  for  the  simple  target  configuration,  although  serious  degradation 
In  sldelobes  Is  most  evident. 

Frames  17  and  18  were  recorded  using  an  approximate  model  of  the  Michigan  11-elemei 
target.  The  target  model  used  is  shown  in  Figure  18,  In  which  the  position  of  eleven  corner 
reflectors  from  the  center  of  rotation  of  the  target  Is  given  by  distance  (in  feet)  and  angle  (In 
degrees).  The  simulation  was  performed  assuming  3  rotational  rate  of  1.  22  degrees/seccnd. 
For  the  simulation,  the  target  was  rotated  through  30°  aspect  change.  Frame  17  was  record¬ 
ed  using  Format  I;  frame  18  was  recorded  using  Format  II. 

Processing  of  frame  17  yielded  no  output.  However,  the  cause  of  this  is  In  the  normaliz¬ 
ation  procedure  used.  When  this  target  Is  at  the  aspect  shown  (with  one  cross  ar;r>  normal  to 
the  llne-of-slght),  seven  of  the  corner  reflectors  are  at  the  same  range  and  a  large  return  is 
observed  at  the  radar.  As  the  target  moves,  there  is  constructive  and  destmctlve  interference* 
and  the  amplitude  of  the  return  is  much  lower  than  for  the  above  aspect.  The  ratio  of  the 
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Figure  18.  Michigan  11  -  Element  Target 

maximum  return  to  the  peak  return  at  other  aspects  is  quit*  large.  Since  the  d.  c.  bias  is 
chosen  according  to  the  maximum  return  and  the  fluctuation  about  this  value  corresponds  to 
the  returns  at  other  aspects,  there  is  in  this  frame  insufficient  contrast  in  the  recorded  in¬ 
formation  to  yield  an  output.  A  simple  expedient  to  avoid  this  problem  is  to  normalize  the 
returns  only  over  an  interval  of  time  corresponding  to  the  integration  interval.  Then  only 
those  intervals  containing  large  returns  are  examined  with  their  corresponding  high  bias. 
Adjacent  intervals  containing  lower  return  amplitudes  are  then  processed  with  a  different 
bias  and  outputs  will  result. 

Frame  18,  recorded  in  Format  II,  couid  be  processed  to  yield  a  recognizable  output. 
Although  the  high  range  sidelobes  are  apparent  in  this  figure  also,  the  target  is  clearly  iden- 

f  ifloKlo  °  * 


3.  PROCESSED  OUTPUTS 

Three  processed  outputs  are  shown  in  this  section,  two  for  dumbbell  conflgui  ition  and 
one  for  the  Michigan  Pi-element  target.  These  are  given  in  Figures  19,  7.V  and  21  .  In 
each  of  these  figures,  a  sample  of  the  rerorded  film  is  shown  and  an  enlargement  of  the  pro¬ 
cessed  output  is  shown.  The  output  films  were  all  made  by  the  University  of  Michigan  The 
high  range  sidelobes  due  to  recording  In  Format  n  are  observed  in  the  processed  output  of 
Figures  20  and  21  . 

4.  CONCLUSIONS 

The  results  obtained  with  the  simulated  targets  demonstrate  ccnclusively  the  feasibility 
of  optically  processing  synthetic  spectrum  data.  Although  adjustments  need  to  be  made  to 
the  present  film  recorder  to  optimize  the  recording  of  Format  II  data,  the  equipment  can  now 
be  used  to  process  data  in  both  Formats  I  and  IT. 
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(;t)  Simple  ol  {Recorded  Film  Frame 


■M 


(d)  Processed  Output  (ISP 
Figure  19.  Frame  9,  Format  I 


(a)  Sample  of  Recorded  Film  Frame 


(1))  Processed  Output  (1ST) 


Figure  20.  Frame  14,  Format  II 
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SECTION  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  EQUIPMENT  CAPABILITIES 

The  film  recorder  equipment  built  on  this  con  -act  has  met  all  design  goals  set  forth  at 
the  beginning  of  the  program  and,  in  addition,  it  has  been  operated  satisfactorily  beyond 
the  design  specifications  In  order  to  inmrove  the  quality  of  the  processed  data.  The  equip¬ 
ment  can  write  from  1024  to  2048  spots  per  line.  These  are  recorded  on  a  50  mm  width  of 
film.  It  has  been  used  to  write  up  to  approximately  25  line  pair  per  millimeter  along  the 
length  of  the  film.  The  film  drive  control  frequency  can  be  as  low  as  600  Hz  (7i4  Hz  was 
design  goal)  and  still  achieve  reliable  operation.  The  stability  of  the  trace  is  regulated  so 
that  reliable  recording  of  the  data  can  be  achieved.  The  equipment  can  be  operated  to  record 
approximately  3550  puises  in  a  continuous  strip  (which,  at  a  line  writing  density  of  25  lines/ 
mm,  is  a  length  of  about  5  1/2  inihes). 

2.  LIMITING  FACTORS 

The  SDS-910  computer  unit  used  with  the  film  recorder  has  only  one  magnetic  tape  unit, 
which  limits  the  number  of  pulses  which  can  be  recorded  continuously  on  film.  However, 
the  fit  ill  recorder  has  no  other  limitation  in  recording  a  continuous  record  up  to  the  length 
of  film  on  a  spool,  which  is  60  feet.  By  use  of  more  than  one  input  tape  unit  these  would  be  no 
difficulty  in  recording  radar  data  even  if  the  number  of  puises  oa  a  run  exceeded  3550. 

At  present,  the  quality  of  the  recording  In  Format  n  is  not  sufficiently  good  to  give 
adequate  dynamic  range  for  processing  returns  from  complex  targets.  Similarly,  a  minor 
modification  needs  to  be  made  to  the  software  for  Format  I  such  that  the  maximum  return 
amplitude  is  obtained  for  data  blocks  corresponding  to  the  Integration  interval,  rather  than 
for  a  complete  run.  This  can  be  easily  accomplished.  For  processing  of  radar  data,  both 
of  these  limitations  will  have  to  be  corrected  to  achieve  best  results  in  the  two  formats. 

To  write  data  more  compactly  than  was  originally  the  design  goal,  the  speed  of  the  film 
drive  unit  was  decreased.  This  required  that  the  control  frequency  to  the  film  drive  unit  be 
lower  than  714  Hz.  In  order  to  obtain  a  Irequency  lower  than  this,  an  external  source  was 
used.  To  continue  recording  at  the  lower  s|X?ed  the  external  source  will  be  required.  Mod¬ 
ification  could  be  made  to  the  timing  circuitry  to  obtain  the  lower  control  frequency  if  con¬ 
tinued  operation  is  desired  here.  Present  and  planned  usage  of  the  film  recorder  using  the 
external  source  will  not  hinder  recording  and  processing  of  radar  data. 

3.  FOLLOW-ON  USE  OF  EQUIPMENT 

Whereas  the  basic  feasibility  of  optically  processing  synthetic  spectrum  data  has  been 
demonstrated  on  this  contract  using  simulated  data,  there  are  several  uses  of  the  equipment 
which  have  not  been  attempted  as  yet.  Foremost,  of  course,  is  the  processing  of  data  re¬ 
corded  by  the  radar,  which  was  not  available  during  the  period  of  this  contract.  Also  a 
number  of  simple  geometric  shapes  could  be  used  in  a  simulation  program  to  demonstrate 
not  only  the  capabilities  of  the  optical  processing  but  also  the  capabilities  of  having  a  time 
history  of  the  output  from  which  the  target  configuration  may  be  deduced. 

It  is  recommenued  that  additional  experimentation  be  continued  to  study  the  two-dim- 
ensionai  imaging  of  geometric  configurations  of  Interest,  especially  with  regard  to  the  time 
history  of  the  correlator  output.  In  addition,  the  film  recorder  and  optical  correlator  should 
be  used  to  process  radar  data.  Utilization  of  the  equipment  in  this  manner  will  provide  im¬ 
portant  information  on  the  quality  of  processed  data  available  and  how  it  may  be  used  most 
efficiently.  In  addition,  some  systematic  procedures  for  information  sorting  and  interpreta¬ 
tion  should  be  gained  by  this  approach. 
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APPENDIX  I 


EQUIPMENT  OPERATION  PROCEDURES 


IA  SDS-910  OPERATIONAL  PROCEDURE 


IB  FILM  RECORDER  OPERATIONAL  PROCEDURE 


APPENDIX  I A 


SDS-910  OPERATIONAL  PROCEDURE 


A.  Action  C:  Comment 

Ai.  Turn  computer  power  on  and  connect  cables  (See  Al.  1  for  details). 

Cl.  Clock  pulses  for  camera  drive  will  be  present  anytime  the  computer  power  is  on. 

A2.  With  Magnetic  Tape  Transport  off,  mount  reel  of  magnetic  tape  and  thread  tape. 

Turn  transport  power  on  and  position  tape  at  load  point.  (See  A2. 1  for  details). 

A3.  Set  transport  switches  to 
Unit  Number  -  3 
Density  -  556 

and  set  AUTO  mode. 

A4.  On  computer,  reset  aii  breakpoint  switches. 

A5.  Load  Optical  Processor  Program  into  computer  from  paper  tape  using  "Standard 
Fill  Procedure".  (See  A5.  1  for  details) 

C5.  Program  is  now  in  the  idling  mode  and  ready  for  operation.  The  mode  of  operation 
is  determined  by  the  breakpoint  settings  as  described  in  next  paragraph. 

A6.  0  Breakpoint  Operation  (See  Table  3. ) 

With  ail  breakpoints  in  the  reset  position,  the  program  is  in  the  idling  mode. 

A6.  1  If  the  operator  wants  to  change  the  timing  between  the  Line  Start  pulses,  he  should  - 

A6.  1.  1  Set  and  then  reset  breakpoint  1.  -  -  - 

The  light  on  the  typewriter  will  come  on. 

A6.  1.2  Type  a  four  digit  number  equal  to  the  desired  number  of  milliseconds.  -  - 
The  program  will  now  automatically  return  to  the  idling  mode. 

C6.  1.2  The  number  of  milliseconds  is  initially  set  to  200  ms.  The  new  number  must  be 
equal  to  or  greater  than  this  number. 

A6.  2  With  Breakpoint  2  Set,  the  operation  is  the  same  as  in  the  normal  operating  mode 

except  that  the  magnetic  tape  is  not  read.  Line  Start  and  Line  Stop  pulses  are  genera 
ted  as  in  the  normal  mode.  Whatever  data  is  in  the  record  storage  area  (e.  g.  the 
last  record  read  from  magnetic  tape)  will  be  output  respectively. 

A6.  3  Set  breakpoint  3  when  normal  operation  is  desired.  Normal  operation  will  continue 
as  long  as  breakpoint  3  is  set  or  until  the  end  of  file  mark  on  the  magnetic  tape  is 
reached.  When  the  end  of  file  mark  is  reached,  ali  outputs  will  cease  until  break¬ 
point  3  is  reset  and  set  again. 

When  the  end  of  the  reel  of  magnetic  tape  is  reached  the  typewriter  will  type  a  mes¬ 
sage  indicating  same  and  the  program  will  halt. 
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TABLE  3.  BREAKPOINT  OPERATION 


Set 

Breakpoint 

No.* 


3 

(Normal) 


2 

(Test) 


1 

(Set  and  Reset 
Immediately) 

(Change  Timing) 


Ope  rat  ion 


Program  alternately  reads  a  block  oil  of  mag¬ 
netic  tape  and  then  outputs  it  at  a  fixed  rate  to 
the  CRT.  Line  start  and  line  stop  pulses  are 
output  at  a  preset  rate. 

Aii  outputs  (except  camera  clock  frequency)  are 
inhibited  at  the  end  of  a  frame  until  Breakpoint 
3  is  reset  and  set  again. 

An  end  of  magnetic  tape  reel  signal  will  'nhibit 
all  further  operation. 

The  magnetic  tape  is  not  read.  The  program 
repeatedly  outputs  the  last  block  that  was  read 
from  magnetic  tape.  All  timing  is  exactly  the 
same  as  when  breakpoint  3  is  set. 

A  new  number  to  preset  the  number  of  milli¬ 
seconds  between  line  start  puiscs  is  accepted 
at  this  time. 


*  Breakpoint  3  inhibits  2  and  1.  Breakpoint  2  inhibits  1. 


A2.  1  MAGNETIC  TAPE 


OPERATING  INSTRUCTIONS 


A2. 1.1  GENERAL. 


This  section  contains  operating  instructions  In-  the  Tape  Transport  and  describes  the 
threading  procedu^res?'  e,eC,r‘Cal  meChan'Cal  fro"'  Pa"el  controls  and  tape 


A2.  1.2  OPERATING  CONTROLS.  INDICATORS,  AND  MECHANISMS. 

All  controls,  indicators  and  mechanisms  required  for  operation  of  the  Tape  Transoort 

is  dlTa  r;  3  ed  °^errfatC:r  ,0.  fr0nt  panel  0f  ,he  ,ape  Thes^ontToZri 

ustea  in  i  able  3.  Interface  signals  for  remote  operation  are  shown  in  Table  4 

Note 


When  changing  speed  or  direction  manually,  it  is  required  that  th.' 
STOP  pushbutton  be  depressed  before  making  a  new  selection. 


TABLE  4.  FRONT  PANEL  CONTROLS 


Control  Pushbutton 


Descript  ion 


Function 


OFF 

(TRANSPORT  OFF) 


Momentary  contact  pushbutton 
switch  with  integral  lamp 


When  the  switch  is  operated  and 
pushbutton  is  lit,  tape  transport 
operation  is  inhibited  except  for 
the  tension  arm  retraction  motor 
OFF  button  may  be  operated  only 
if  machine  Is  first  placed  in 
STANDBY  condition. 


STOP 


FAST  FORWARD 


FORWARD 


REVERSE 


Momentary  contact  pushbutton 
switch  with  Integral  lamp 


Locking  pushbutton  switch 
with  integral  lamp 


Locking  pushbutton  switch 
with  integral  lamp 


Locking  pushbutton  switch 
with  integral  lamp 


When  switch  is  operated  and  push¬ 
button  is  lit,  tape  transport  is  en¬ 
ables  and  arm  retraction  motor 
is  inhibited.  Any  previously  opera¬ 
ted  locking  pushbutton  switches 
are  released. 

Provides  manual  selection  of  tape 
motion  In  the  forward  direction  at 
high  speed. 

Provides  manual  selection  of  tape 
motion  in  the  forward  direction  at 
normal  speed. 

Provides  manual  selection  of  tape 
motion  in  the  reverse  direction  at 
normal  speed. 
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TABLE  4.  FRONT  PANEL  CONTROLS  (Cont) 


Control  Pushbutton 


Description 


Function 


FAST  REVERSE 


Locking  pushbutton  switch  Provides  manual  selection  of  tape 

with  integral  lamp  motion  In  the  reverse  direction  at 

high  speed. 


AUTO 


Locking  pushbutton  switch  Transfer  all  control  and  indicator 

with  integral  lamp  functions  to  remote  locations. 


A2.  1.  3  TURNING  POWER  ON. 

Before  tape  threading  can  take  place,  the  circuits  associated  with  automatic  tension 
arm  retraction  must  be  energized. 

Power  is  remotely  controlled  by  computer  power  switch. 

a.  Threading  Tape 

When  threading  tape,  it  is  recommended  that  only  the  end  of  the  tape  is  handled.  Avoid 
touching  the  tape  if  any  other  point. 

The  steps  in  paragraph  3-3  must  be  satisfied  first  before  proceeding  to  load  tape. 


CAUTION 

To  prevent  tape  break  at  end  of  run,  do  NOT  slip  free  end  of  tape 
into  slot  in  take-up  reel  core  or  in  any  way  fasten  the  free  end  to 
the  take-up  reel. 

b.  Load  Point  Marker.  It  is  desired  to  begin  operating  at  the  marker,  press  the 
FORWARD  pushbutton  on  control  unit  until  marker  reaches  the  take-up  reel.  Press  STOP 
and  then  REVERSE  pushbutton.  Tape  will  reverse  direction  and  stop  on  the  LOAD  POINT 
marker. 


c.  Tape  Removal.  When  reflective  markers  are  used,  the  tape  will  stop  on  these 
markers  when  they  reach  the  sensor  and  the  driving  pinch  roller  will  be  released.  The  trans¬ 
port  will  be  inhibited  from  driving  tape  passed  the  marker  once  it  has  been  stopped  by  this 
marker. 


CAUTION 

When  the  tape  has  been  stopped  from  a  fast  mode  of  operation, 
it  is  necessary  to  wait  for  two  to  five  seconds  before  selecting 
another  mode.  This  will  allow  time  for  the  capstans  to  slow 
down  to  normal  speed. 
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To  remove  tape,  select  a  fast  reverse  or  reverse  mode  until  tape  reaches  LOAD  POINT 
marker.  If  nc  marker  is  used,  stop  the  transport  when  approximately  ten  feet  of  tape  re¬ 
mains  on  the  pay-out  reel, 

•  Press  OFF  pushbutton  on  control  unit 

•  Rotate  reels  by  hand  until  tension  arm  rollers  just  pass  the  stationary  rollers 

•  Open  read/write  head  cover,  vacuum  buffer  cover  and  press  arm  retraction  button 

•  When  tension  arms  stop,  wind  the  remaining  tape  onto  the  take-up  reel  while 
turning  the  pay-out  reel  to  avoid  stretching  the  tape, 

•  Close  read -write  head  cover,  vacuum  buffer  cover  and  remove  tape  reel. 

AS.  1  LOADING  OPTICAL  PROCESSOR  PROGRAM 

A5,  1,  1  Standard  Fill  Procedure 

a)  Set  Operating  Switch  to  IDLE. 

b)  Press  Start  Button 

c)  Mount  paper  tape  in  photo-reader.  * 

d)  Place  Operating  Switch  in  RUN  position 

e)  Raise  and  release  Fill  Switch. 

A5.  1.2  Standard  Restart  Procedure 


a)  Set  Operating  Switch  to  IDLE. 

b)  Press  Start  Button 

c)  Move  the  Operating  Switch  to  the  STEP  position  and  then  to  the  RUN  position. 
A5.  1.  3  Standard  Turn-Off  Procedure 

a)  Set  Operating  Switch  to  IDLE 

b)  Press  Start  Button 

c)  Press  POWER  Button. 


*  Tape  feeds  from  left  to  right,  in  the  direction  of  the  arrow  on  the  tape.  The  "TOP"  side 
of  the  tape  should  be  up.  The  tape  should  lie  flat  in  the  reader  under  all  guide  posts.  Knob 
on  photoreader  should  be  turned  from  LOAD  to  RUN  after  mounting  tape. 
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APPENDIX  IB 


FILM  RECORDER  OPERATIONAL  PROCEDURE 

Al.  Turn  power  on  SDS -9 10  computer 

Cl.  See  operational  procedure  for  computer 

Cl.  1  Obtain  start -stop  pulses  from  computer  (set  Break  Point  2). 

A2.  Turn  power  on  for  cathode  ray  tube  assembly  and  film  drive  power  supplies. 

C2.  Allow  a  minimum  of  5  minutes  to  obtain  thermal  stability  in  equipment. 

A3.  Turn  on  power  for  film  drive  irequency  reference. 

C3.  This  signal  was  previously  supplied  by  SDS -910.  However  to  obtain  greater  density 

in  line  spacing  this  modification  was  incorporated  in  the  operating  procedure. 

A4.  Set  video  gain,  calibrate  reference  and  brightness  controls  to  desired  values. 

C4.  Desired  values  are  obtained  from  figures  15  and  16. 

A5.  Set  film  drive  range,  push  film  drive  motor  control  and  synchronize  film  drive  to 

desired  film  speed. 

C4.  Film  drive  speed  now  controlled  by  film  drive  frequency  reference. 

A6.  Remove  start -stop  pulses  from  computer  after  examination  of  video  signal. 

A7.  Set  breakpoint  3  on  computer. 

C7.  1  This  will  start  recording  of  data  on  film. 

C7.2  Should  a  missing  line  stop  pulse  occur,  reset  computer  as  quickly  as  possible. 

A8.  When  frame  has  been  computed,  push  fast  film  advance. 

C8.  Allow  to  run  for  5-10  seconds  to  give  a  trailing  leader  on  the  film. 

A9.  Push  film  drive  motor  control. 

C9.  This  will  stop  film. 

A10.  Remove  camera  assembly  and  develop  film. 
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APPENDIX  D 


FORMAT  ERROR  ANALYSIS 


ERROR  ANALYSIS  FOR  THE 
OPTICAL  RECORDING  AND  PROCESSING 
OF  SYNTHETIC  SPECTRUM  DATA 


bv  T  C  Grlskey 
Au^-ust  5.  1966 


The  work  reported  herein  was  performed  on 
contract  AF30(602)-4287.  ARPA  Order  837. 


Westlnghouse  Defense  and  Space  Center 
Surface  Division 
P.  O.  Box  1897 
Bal.imorc,  Maryland  21203 


ERROR  ANALYSIS  FOR  THE  OPTICAL  RECORDING  AND 
PROCESSING  OF  SYNTHETIC  SPECTRUM  DATA 

L  INTRODUCTION 


The  purpose  of  this  report  is  to  determine  the  errors  associated  with  the  optical  recording 
and  processing  of  synthetic  spectrum  data  The  sidelobe  degt  'dation  due  to  the  inaccuracies 
of  the  system  are  calculated  for  two  data  formats.  The  format  that  yields  the  best  system 
performance  based  on  its  inherent  sidelobe  degradation  is  determined  from  this  analysis. 

A  brief  mathematical  sketch  erf  the  system  is  performed  to  describe  the  system.  The  nota¬ 
tion  of  P  M.  Woodward  will  be  used  throughout  this  analysis.  1  This  notation  is  given  in 
Table  L 


H.  FORMAT  DESCRIPTION 


1.  GENERAL 

A  brief  introduction  to  the  general  format  of  the  recorded  Information  will  be  helpful  in 
understanding  the  developments  which  follow.  Information  obtained  on  one  radar  transmitted 
pulse  is  recorded  on  one  line  across  the  film.  The  number  of  (bta  points  per  pulse  and  the 
number  of  'spots*  per  line  is  selectable  and  is  denoted  by  N.  Information  from  two  succes¬ 
sive  pulses  is  recorded  on  successive  lines  adjacent  to  one  another  along  the  length  of  the 
film.  Variation  of  the  line  density  is  possible  in  the  equipment. 

Each  line  recorded  on  fllru  is  positioned  to  be  perpendicular  to  the  edge  of  the  film  by 
adjusting  the  yoke  of  the  CRT.  Information  from  each  pulse  Is  intensity -modulated  onto  the 
CRT  phosp^xjr  on  a  line  which  is  regulated  quite  precisely.  Every  pulse  is  written  onto  the 
same  11  ic  on  the  CRT  phosphor.  Serration  of  adjacent  lines  on  the  film  is  controlled  by 
having  a  uniform  film  drive  s|X?ed  and  delaying  the  start  of  writing  each  line  of  information. 

Information  can  be  recorded  on  film  using  two  formats.  In  Format  L  the  amplitude 
information  from  a  pulse  is  uted  to  intensity -modulate  the  CRT  and  expose  the  film.  In 
Format  n.  the  Fourier  transform  of  the  amplitude  Information  from  a  pulse  is  used  to  lnten- 
slty-modulate  the  CRT  and  expose  the  film.  This  memorandum  considers  the  errors  for 
these  two  formats: 

2.  FORMAT  I 


In  Format  I.  the  amplitude  of  a  block  of 
be  represented  by: 

0)  T  (x.  y)  rect  rect 

I  (x  '  •  y)  *  rect 

*  A (y )  cos  f  2rfy  ♦  0(y) 


radar  pulses  on  the  negative  transparency  can 
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TABLE  5.  WOODWARD'S  NOTATION 


Definitions: 


u(t)  =  f  u(T)  J(t-T)  dr  --  J  U(l)  el 
-gd  -ao 

ao 

u  *  v  J  u(r)  v(t-  r )  dr 

-CD 

U(f)  --  J  u(t)  e*2**1  dt 
-ao 

ao 

repT  u(t)  =  £  u(t-nT) 

-ao 

ao 

conibF  U(f)  -  £  U(nF)  J(f-nF) 


reel  t  = 


sine  f  - 


1.  |t|<{ 

0.  It  I  >  ' 


sin  *f 


Fourier  Transform  Pairs: 

Waveform  |  Spectrum 


U(f) 

c  U  ♦  c,V 
i  Z 

U(-f) 

u*(f) 

I2"f  U(f) 


Waveform  Spectrum 


uv 

U*V 

rcpT  u 

1 1/T  |  combj  T  U 

combT  u 

|l/T|  rop|/T  U 

U(t) 

u(-f) 

8  (t) 

1 

U(f)  exp{-i23  r f) 


u(t)  cxp(i2?0t) 


exp  (-fft  )  |  exp  (-rf") 


where 


r  (x,  y)  is  the  ratio  of  the  light  amplitude  after  being  modulated  by  the  film  to  the 
light  amplitude  before  being  modilated  by  the  film 

x  is  distance  measured  a  long  the  film 

y  is  distance  measured  across  the  film  (and  is  proportional  to  range) 

Y  is  the  width  of  the  film  frame  being  used  in  both  the  x  and  y  directions 
A  is  determined  by  the  exposure  of  the  background  on  the  film 
A j  is  determined  by  the  exposure  associated  with  the  information  on  the  film 

M  nu«aaer  of  pulses  per  film  frame 
N'  nunuer  of  data  points  per  pulse 
P(x)  =  persistence  of  the  CRT  phosphor 
I(x,  y)  amplitude  distribution  of  a  dot  on  the  Him  due  to  a  dot  on  the  CRT 
Ao  '  I A^  nax  =  'ittxtmum  puise  amplitude 

A(y)  =  ampiilude  information  obtained  from  a  pulse 
fl(y)  =  phase  shift  between  pulses  due  to  doppier 
f  -  carrier  frequency. 

The  function,  ^  reel  ( ^  j  )  rect  (y  ]  ,  is  unity  within  the  data  frame  and  zero  out¬ 
side  it.  The  delta  function,  8  (y  -  )•  represents  sampling  of  the  ctata,  A  ♦  A(y)  cos 

[2rfy  *  0(_-)  ] ,  at  points  separated  by  Y/N.  The  function,  rect  (*£’)  Indicates  tint  these 

sampled  values  exist  on  the  Him  for  a  distance  Y/N.  The  characteristics  of  the  CRT  are 
taken  into  account  by  convolving  the  persistence  of  the  phosphor  P(x)  and  the  distribution  of 
the  Intensity  of  a  dot  I(x,  y)  with  this  data. 

A  possible  configuration  of  the  (simplified)  optical  processor  Is  shown  in  Figure  22. 

The  Fourier  transform  of  the  amplitude  Information  is  performed  in  both  the  x  and  y  direc¬ 
tions.  A  one -dimensional  transformation  is  then  taken  on  one  of  the  sidebands  of  the  infor¬ 
mation.  A  sideband  is  used  here  to  allow  recovery  of  both  amplitude  and  phase  information 
and  to  eliminate  some  of  the  noise  from  the  system. 

The  two  dimensional  Fourier  ‘ransferm  of  the  data  can  be  represented  by 

W  E(v,  u)  Kj  r  (x,  y )  exp  [-J2s(uy  -  vx)  1  dx  dy 

-oo  1  1 

where: 

v  =  x,/XF 
u  =  yj/XF 

X  *  wavelength  of  the  coherent  light  used  in  the  processor 
F  =  focai  length  of  the  transforming  ions 
*i .  y i  *  new  dimensions  in  the  transform  plane. 
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FINE 

RANGE 


To  illustrate  the  operation  of  the  optical  processor  let  us  assume  that  there  are  two 
point  targets,  separated  a  distance  T,  being  detected  by  the  radar  and  that  one  is  rotating 
about  the  other. 


The  amplitude  transmissivity  of  the  film  for  this  particular  case  will  be  given  bv: 

M  N 

(3)  T(x,y)  rect  (y  -  j)  rcct  (^  -  j)  A-Aj  E  E  fp(x)  • 

m*l  n-1  ^ 

. 

1  *x  ‘  ItT’  y*  *  rcct  *^Jao  *  s^'c  £ \  (y- 
cos  (2-fy  ♦  niP)  ♦  sine  (y-  ^)j  cos  (2rfy)j  &(y-^D>] 

This  assumes  that  the  amplitude  information  is  of  the  form  sine  (^-)  and  a  range  cell  is  Z 

units  long  on  the  film.  The  general  appearance  of  this  Information  on  the  film  is  illustrated 
by  Figure  23. 


Carrying  out  the  operation  indicated  by  Equation  (1)  we  have: 


(4) 


E(v.u)  *  I  sine  Yv  exp(-jsYv) 
M  N  s 


H 


sine  Yu  cxp(-JrYu) 
rnY 


)•[■ 


A  8  (u)  8  (v) 

,  I  I  [Pj(v)  Ijfv.u)  oxp(-Jl'irv  n^)  sine  (Ju/(a  8  <u)  8  (v) 
m*l  n=l  ^  ‘  \l 

♦  rcct  ( )  exp  j^-j=(Yu  ♦  2T  u)  ]  •  8  (u-f)  cxpO  mP) 

♦  ^  8  (u*f)  cxp(-J  nip)  J  ♦  rcct  ( exp  (-)’fYu)  • 

|  j  8  (u +f)  ♦  |  8  (u-f)]  |  •  exp(-]2ru  jjf  )^>]J 
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Performing  the  summation  we  have. 


E(v,  u)  =  Kj  r  sine  Yv  exp<-jsYv)  ]  •  [  sine  Yu  exp  (-JrYu)  j  • 

£A  S  (u)  8  (v)  -  Aj  j^Pj(v)  Ij(v,u)slnc  (^j)<^?Ao5(u)  8  (v) 
^yVJ)  exP  pj(M«-l)r^vJ  ♦  rect(^)exp  "'-JfffY^THiJ 


8  (u-f) 


sin  M:(^  -  ^y) 

8in  r  {h  *  H  v) 


1  sin  Mr  (1  ♦  I  V)  r  „  v  ] 

♦  s  8  (v.*f) - -y  -yM  exp  z  (£-  ♦  I  v)  1  > 

sin  s  ^  v)  L  2r  M  J  J 

rect  (  j)  exp(-)rYu)  sJP(*Vv>  exp  r  ^  v  ] 

sin(rf.v)  L  M  J 


sin  v  {^t  ♦  gj  v) 


*  [2  8(u‘f)  5  <u*f)l  I  •  8J?_i£Yu)  exp  [jr(N*l)^ul  ^)1 

J  '  sin  (r-^u)  L  ni  j  j 

y 

The  term  sin(~Yu)  sin(?^u)  eauses  the  spectrum  to  be  repeated  in  the  u  direction  at 

integer  multiples  of  N  Y.  This  is  illustrated  in  Figure  23.  where  the  DC  term  and  the  other 
information  on  Uk  film  appear  as  bright  spots.  An  aperture  can  be  placed  about  the  lines  at 
V'  0r  1!  r  f  k  ‘1,U8tratlon  it  is  placed  at  u  *  -f.  A  cylindrical  lens  is  now  placed 
at  the  other  side  cf  this  aperture  and  a  one  dimensional  Fourier  transform  is  taken  in  the  u 
dlrecMon.  The  amplitude  information  on  the  other  side  of  the  aperture  can  be  approximated 

(6)  E(v.'u)  =  Kj  rect  (Uy)  |^slnc  Yv  exp(-)irYv)  J  •  ^stnc  Yu  exp(-j»Yu)J 

*  [p,(v)  l,(v,u)sinc(£u)(  |  rect  J  exp  £ -Jj(Y*2T)(u*f)J 

sin  Ms  (£-  ♦  n  v)  r  . 

s^Tr  <£.  *  V  v)  exp  *  (EF  *1  v>]  +  rect[Z-^jexp  [-J«rY(u+f)j 

—f  exp  8i5i£Yy}  exp  [j(N.l)rIul  >1 

8  n  *  v  L  JJ  stn(iXu)  L  N  JV 


sin  ^  v 
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"AOA*  PutStS  —x 

(a) 


toe 

2TII  2n7 

(b) 


toe 

2  m  2nt 

(c) 


Figure  23.  Optical  Processing  Using  the  Amplitude  Information  for  Each  Pulse 
Now  taking  the  one  dimensional  transformation  in  the  u  direction,  we  have: 


(7) 


E(v.y  )  =  K3  |slnc(Xy  )  exp(-2sjy  f)J  •  ^sinc(Yv)  exp(-jrYv)  j  •  reel  -  g) 

j^Pj(v)  lj(v)l(y  )  •  reett^-)  •  jsinc  (y  -  T)  j  expH2?y  f) 
sin  M*  (£:  ♦  n  v)  [  Y  “I  f2  •  v  “I 

;  > -v  -  exp  h<M<I)I(^  *rv)J  48fac|z(y  -jM 

sin  “  *  M  v  L  J  L  J 

exp( *j2sy  f)  8i-n  *  y-  exp  f-J(M-l)  ary  vl  E  5(y  -  "I))] 

sin s  v  L  J  n=l  ' ' 


This  is  displayed  on  the  film  shown  in  Figure  23.  It  is  seen  that  the  two  targets  are 
now  separated  It.  range  and  doppler. 
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c. 


FORMAT  D 


The  processing  of  the  data  using  Format  II  is  different  from  that  using  Format  I  in  that 
the  Fourier  transform  of  the  amplitude  information  is  initially  written  on  the  film.  The  am¬ 
plitude  information  on  the  film  for  a  block  of  data  can  be  represented  by: 


(8) 


T  (x, y)  rect  (^  -  |)  reel  (^  -  i) 


M  N 


HI  iw  > 

A  -  A  £  £  I  P(x) 

1  m=l  n=l  V 


1  (x  *  y)  *  rect ($).<| 

“*  [2.(f.Tp)y.n»pj|S(y.f))J 


P 

A  ♦  £ 

°  P=1 


ap(y) 


Where  a  (y)  is  the  Fourier  transform  of  the  amplitude  information  from  the  pth  target, 
Tp  represents*5 the  delay  Information  from  the  pth  target,  and  the  other  parameters~are  as 
defined  in  Equation  (1). 


The  optical  processor  in  this  case  is  shown  in  Figure  24.  Just  one  two-dimensional 
transformation  is  taken  in  this  case. 


Let  us  examine  the  example  illustrated  in  the  analysis  of  Format  L  Assume  that  there 
are  two  point  targets  spaced  a  distance  proportional  to  T  and  that  one  is  moving  with  respect 
to  the  other  such  that  there  is  a  dopplrr  shift  of  0  radians  between  each  pulse.  The  informa¬ 
tion  on  the  film  can  now  be  represented  by: 


Figure  24.  Simplified  Processor  Diagrams 
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This  information  is  illustrated  by  Figure  25.  Now  talcing  the  two  dimensional  Fourier 
transform  indicated  by  Equation  (2)  we  have : 


E(v.u) 


Ki[5 


>]*[ 


sinc(Yv)  exp(-jrYv)  •  sine  (Yu)  exp(-jsYu) 


M  N 

A  8  (u)  S  (v)  -  Aj  Y/  X/  fPjfvHjtv.u)  exp(-j2r^v)  sine  (^r) 
m=l  n-1 


^|Ao  8  (u)  8(v)  +  ^sinc(Yu)  exp(-jsYu)J  *  ^  8  (u-f-T)  exp(jnfl) 
+  \  8  (u+f+T)  exp(-jnfl)  ♦  ^  8  (u+f)  ♦  \  &  (u4)J  |  * 


exp  (*j2i7 


Carrying  out  the  summation,  we  have: 


E(v,u)  Kj  Jj 

j^A  8 


sinc(Yv)  exp(-j”Yv 


sinc(Yu)  exp(-jjrYu) 


A  S(u)  8(v)  -  A 


j  ^Pj(v)  Ij(v,u)  sinc(^)  ^  Aq  8  (u) 


-  -^V— -  exp  r-js(M+l)£v  1  ♦  sine  Y(u-f-T)  exp  f-jsY(n-f-T)l 

sin  (x^  v)  L  M  J  L  J 


sin  I1  (ijv-fl) 


sin  2  (it  j^v-fl) 


exp  -  j77  Y(n+f  *T) 


’-j(N^)  (2-^v-fl)  1  ♦  sine  Y(u*f+T) 


sin  -j  (2uj^  v+0) 


sin  2 


+  sine  Y(u-f) 


sin  (~Yv) 
sin  ("|^v) 


sin  it  Yv 
Y 

sin  17  —  v 
M 


exp  r.!(^H2^v.J 

H’*«  L  J 

exp  £-jir(M  +  l>^  v  J  +  sine  Y(u+f) 


[•WM.^v]))].{«gexp[WN.>,|U]] 


This  information  now  appears  as  bright  spots  separated  in  range  and  doppler.  This  is 
illustrated  in  Figure  25.  To  remove  the  ambiguities,  this  information  is  now  passed  through 
an  aperture  located  in  the  region  ofu-ftou-f  +  T  and  is  imaged  onto  the  film.  The  out¬ 
put  film  is  shown  in  Figure  25. 
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(0)  ( b)  (c) 


Figure  25.  Optical  Processing  Using  the  Fourier  Transform  of  the  Amplitude 

Information  for  Each  Pulse 


111.  ERROR  ANALYSIS 


1.  FORMAT  11 


The  system  Inaccuracies  that  will  degrade  the  processed!  data  will  be  mainly  associated 
uonlim'arlties  of  the  film  drive  and  CRT.  Using  equation  (91  it  Is  possible  to  Introduce  errors 
In  the  positioning  of  the  data  points.  The  data  on  the  film  with  errors  can  be  represented  by: 


(12) 


r(x,y)  - 


rcet  (^  -  rect  (^ 


M 


l 

m«l 


N 


I 

n»l 
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where: 


$  represents  the  error  In  positioning  the  mth  pulse  in  the  x  dire  'i'll  across 
m  the  film 

g  represents  the  error  In  positioning  the  nth  data  point  of  the  mth  pulse  In  the 
m,n  y  direction  along  the  (Urn 

A  represents  the  fractional  variation  I.,  the  amplitude  Information  on  the  film 
m,  n 


(13) 


Taking  the  two-dimensional  Fourier  transform  of  this  film  we  have: 

E(v,u)  -  [sine  Yu  exp(-J*Yv)]  •  [sine  Yu  exp(-)*vu)J  •  £a  8  (u)  8  (v) 

M  N  f 

-A  [  P  (v)  1  (u,v)  E  E  sine  (^*8  n>u  4  8  (u)  8  (v) 

V  m«l  n«l  l 

♦  cos  [2s(f*T)^  *  m0  ]  ♦  cos  (2af^)  jexp  [-J2»v(^  ♦  8m)J 

•  (A  ♦  A  A  )  exp  [-J2*u  (— ♦  S  H  ]  | 

m,n  m,n  m,n  L  N  m.nJ^J 


The  functions  slnc(Yv)  cxp(-JaYv)  and  slnc(Yu)  exp(-)sYu)  act  like  Impulses  being  convolved 
with  the  remainder  of  the  Information.  Hence  a  good  approximation  to  this  expression  Is: 


M  N 

(14)  E(v,  u)  ■  AS  (u)  5  (v)  -  A  fp^v)  I^v.u)  I  Z  sine  (^  ♦  S^u 

k  m*l  n»l 

<  8  (u)  8(v)  .  cos  [2s(f»T)^»  m0  "|  *  cos  (2«f^) 

o.p  [-12=  •  *„,>*] 

e«P 

We  arc  not  concerned  with  the  Information  at  u  ■  0  and  v  •  0  (on  the  optic  axis).  Also, 

v 

the  (unctions  Pj(v),  Ij(v, u)  and  slnc(^  ♦  8fl)u  act  like  wideband  filters  passing  all  the  In¬ 
formation  of  Interest.  Hence,  these  terms  may  be  neglected  such  that  we  have: 
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Since  the  process  of  taking  the  two-dimensional  transformation  la  a  linear  process  we  need 
only  consider  one  of  the  target*  here,  keeping  in  mind  that  another  target  can  be  added  at  any 
time  using  the  principle  of  superposition.  Using  the  target  that  Is  delayed  In  range  by  T, 
we  have 


(16) 


M  N 


"P  [-,J*  (W  •  8mlv  ] 


All  the  information  concerning  this  target  is  located  at  u  ■  a  (f*T).  It  is  only  necessary 
to  consider  one  of  these  sidebands  in  order  to  determine  the  sidelobe  degradation  associated 
with  the  sy-.tem  Inaccuracies.  Using  this,  Equation  (16)  may  be  rewritten  as: 


(17) 


E(v,u)  - 


M 


K  I 
m*  1 


N 


l 

n«l 


(A 

m,  n 


♦  A 

m,  n 


A 

m,  n 


)  exp 


J^l2r  [(l*T-u) 


nY 

N 


.Yv  - 

-  ( IT  -  W  m  - 


We  will  be  mainly  concerned  with  the  envelope  of  this  function  in  die  sidelobe  region. 
Since  it  is  practical  to  measure  the  sideiobes  with  respect  to  the  peak  at  the  main  lobe,  this 
function  will  be  normalized  with  respect  to  its  magnitude  at  u  -  f*T  and  v  .  -0M.  We  can 
approximate  this  by:  Y 


(18) 


m 


M  N 

I  I 

1  n-1 


(A 


in,  n 


♦  A 


m 


nA . 


M 

l 

m*l 


N 

I 

n  i 


m,  n 


It  is  assumed  that  the  probability  distribution  of  the  random  errors  A  8  and 

m,  n  m 

fim,n  C#n  **  rcpre8cn,cd  P(Am,  „>•  p(  p(  „)  each  having  a  mean  of  zero. 

This  means  that  the  film  drive  has  a  velocity  that  varies  about  the  desired  operating  velocity 
In  a  random  manner,  the  CRT  has  an  intensity  variation  that  varies  about  the  desired  inten¬ 
sity  in  a  random  manner,  and  the  errors  In  the  sweep  on  the  CRT  are  mainly  due  to  the 
sweep  drive  and  not  the  triggers  that  start  the  sweep,  thereby  causing  a  random  error  In  the 
displacement  of  the  data  points  about  their  desired  positions.  Hence,  the  probability  of  a 

particular  value  of  A  ,8  and  8  is  given  by 
m, n  m,  m, n  7 


(,9>  P(Am  n’  8m  8m  n]  “  P(Am  P<8m>  P<  8  > 

m,  n  m,  m,n  m,  n  m,n 

The  standard  deviation  of  the  envelope  |  E  |  can  be  represented  by: 


(20) 
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where: 


on  |e|  2 .  |E|!wvp|vp".,»>dV,«i 

(22)  ]i[  ■  ///  |e|  P(4„  „)  P(»J 

In  order  to  calculate  the  above  expression  it  is  necessary  to  examine  the  intensity  of  the  In¬ 
formation  contained  in  equation  (18).  This  is  given  by: 

/  N 

M  N 

£  £  (A  ♦  A  A  )  r  \ 

.  „  .  **,  m,  n  m.n  m, n  r  „v  v«  1 

(23)  |E|2.  ro-‘  gi  T, - 

r  r  »  ' 


I  Z  A 

m-1  n-1 

M  N 


y  y  (A  A  ) 

L  L  m.n  m.n  m,  n 
n«l  n«l  *  ’ 


M  N 

I  I  Amn 

mil  n  1  ’ 


m.n  r 

- sln<  2i  (i 


((♦T-u)^-£f -D-  8  v-  8 
N  M  mm 


vj) 


Expanding  the  cosine  and  sine  in  this  expression  we  have: 
( M  N 


y 

I  <A 

< 

< 

♦ 

C 

lEl  2  - 

i-* 

m-l 

n-1 

m,  n 

m, n  m, n 

jcos  2s  r(UT-u)^-(~0)ml 

|E|  - 

M 

N 

L 

L 

A 

m  n 

[  L  J 

["-2s(  8  v  *  8  u)l  ♦  sin  2r  r(f*T-u)^-£j-0)m|  sin[2::(  8  v*  8  u)l 
L  m  m.nj  L  NM  J  [_  mm,  n  'J 


M  N 

£  I  (A  «A  A  ) 
m  i  n  i  m’n  m-n  m.n 

M  N 
I  I  A 
m-1  n*l  ’ 


—  |sin  2»  ^(f*T-u)^-(^-0)m")cos  £-2r(8 


v*  8 

m  m, 


.CO.  2.  [(r.T.u)$.(£.(0n,]  .In  [■!■( 

It  is  assumed  that  these  errors  are  small  and  using  a  small  angle  approximation  we 


M  N 

I  Z  (A  ♦  *  A  ) 


/  A  in  vu  n  i 

.  0  _  .  .  '  ro,n  m.n  m.n'  r  „  „  T 

lEl  ■  M  H -  '  co.2.[(f.T-o)=f.(|iv.(r)mJ 

l  l  A  L 


m=l  n=l  ro,n 
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*2«  V*8m,nU)  tln2t  [('♦T-u)£-<£-flra]jj 


I  Z  (Aran*AronAmJ 

m*l  n-1  m'n  m'n  ra'n 


M  N 

Z  Z  A 


|  sin  2 1  |(f*T-u)|~(|J-<l)n>  ] 


l 


m  l  n*l 


m,n 


-  2,<8mV*  8m 

m  ra 


[«-T-“>7T-<£  -«n.]J] 


Expanding  this  expression  we  have: 

1 


(26) 


El2  . 


5M  N 

Jl  n^l  (Am.n'Am,n  Am,n)  cos  [Wt-u)^ 


-  2.( 


-2r( 
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/  M  N  N 

(nM  n^lAm-n) 

{M  N 

J,  /,  j 

(M  N 

y  y  (A  A  )  [2*(  8l_v*  8  u)l  sin  [2s(f*T-u)!~ 

m.|  n*l  m>n  m'n  ro'n  L  ^  ra’n  J  L  N 
^  -B)m  J  } 

‘ill  (A™^*^.nAm.n),,n[2'(',T-")F'2,(S!-H  J' 

•  2(  j  ,  j,  <*„,„•  V»  *„.»>  *"  [2XI-T-u)=?-i.(S-«n.]  J 

{J,  [•2'<V*»m^]“,[2'<,*T-“'^-2'<S-|,>“]} 

•{J,  J,  <An.yAm^Am,»[-2,<»n.''*»n./l  ”*  [2^.T-U)f  -24'-»mjJj 


I 


Using  equation  (21)  and  carrying  out  the  indicated  integrations  keeping  in  mind  that  the  random 
parameters  have  a  zero  mean,  we  have: 


(27) 
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Rearranging  these  terms  we  have: 
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X  wavelength  of  the  coherent  light  used  in  the  processor 
F  *  focal  length  of  the  processor 
Xj,  y j  3  new  dimensions  in  the  image  plane 

<7g  =  standard  deviation  of  the  envelope 

2 

This  equation  is  similar  to  that  developed  by  Ruze  who  showed  that  the  variance  of  the 
envelope  is  related  to  the  variance  of  the  amplitude  r.id  phase  errors  for  a  phased  array  antenna. 

The  statistical  behavior  of  the  sidelobes,  due  to  amplitude  and  phase  errors,  for  this 
normalized  envelope  of  equation  (18)  is  described  by  the  modified  Rayleigh  distribution. 
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The  probability  density  function  of  the  modified  Rayleigh  distribution  is  given  by: 


(31)  p(r)dr  »  *4  exp  )  !o  (  )  dr 

where 

p(r)dr  =  probability  that  the  sidelobes  of  the  envelope  lies  r  and  r  +  dr 
*  mean-square  value  of  r 

lQ  «  modified  Dcssel  function  of  the  first  kind 
Good  values  to  use  for  v  and  u  in  equation  (30)  for  a ■  are  those  where  the  peak  of  the  infor¬ 
mation  is  contained.  This  occurs  at  v  =  L  u  *  f  ♦  T,  where  0  is  the  phase  shift  from  pulse 

y 

to  pulse  due  to  doppier,  q  is  the  distance  between  pulses  on  the  film,  f  is  the  spatial  carrier 

at  which  the  information  is  written,  and  T  is  the  range  delay  of  the  object  Combining  all  of 
this  information  we  have: 
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A  more  detailed  explanation  is  required  to  define  the  standard  deviation  of  the  error 
a  .  This  error  contains  velocity  errors  due  to  the  velocity  inaccuracies  of  the  film  dr've 

and  the  CRT.  The  yoke  of  the  tube  is  set  such  that  the  lines  on  the  film  are  at  a  right  angle  to 
the  edge  of  the  film.  The  distance  between  the  beginning  and  ending  of  the  pulse  along  the 
direction  of  motion  of  the  film  is  given  by: 

(33)  x  vyt  -  |  vgt  sin  0 

where 

Vj  =  velocity  of  film 

Vg  =  sweep  velocity  of  the  CRT 

8  *  angle  at  which  the  line  is  skewed  on  the  CRT 

The  scan  velocity  of  the  CRT  is  divided  by  two  since  it  undergoes  a  two  to  one  reduction 
in  size  when  it  is  imaged  on  the  film.  In  equation  (33)  the  angle  0  is  such  that  x  =  0  when 
there  are  no  errors  in  the  velocity  controls  of  the  film  drive  or  the  sweep  drive  of  the  CRT. 
The  sweep  velocity  of  the  CRT  will  be  much  greater  than  the  film  drive  velocity.  Hence,  the 
function  sin  0  can  be  approximated  by:  sin  0  =  2  v^/v 
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where 


Vj  =  mean  velocity  of  the  film 
vg  -  mean  scan  rate  of  the  CRT. 

Substituting  this  into  equation  (33)  we  have: 

(34)  x-(vf-v.^)t 

Therefore  an  error  in  vf  of  8  vf  and  an  error  in  vg  of  8  vg  yields  a  position  error  of  8  .  which 
can  be  represented  by : 


If  the  time  is  set  equal  to  the  time  between  traces  we  have: 

(36)  8  =  8f  -  8  ? 

X  y  v8 

where 

8  f  =  8vft  *  displacement  inaccuracy  of  the  film  drive  from  pulse  to  pulse 
8  y  =  8  V  *  crror  ln  1,10  r*P°atibllity  of  the  CRT  from  pulse  to  pulse 
«x  *  error  in  spacing  of  the  lines 
Calculating  the  var Lance  of  equation  (36)  we  have: 


where 


am  *  8tandard  deviation  of  the  displacement  between  pulses  on  the  film 

%i.  n  *  standard  deviation  of  the  rcpeatlbUtty  of  the  traces  on  the  CRT 

<t  f  -  standard  deviation  of  the  positioning  inaccuracy  of  the  film  drive 

The  angular  accuracy  of  the  film  drive  is  usually  given  In  angular  measurements.  Using 
this  fact  and  combining  the  above  equation  with  our  original  expression  for  the  variance  of  the 
jldelobes  ol  the  envelope  ln  equation  (32)  we  have: 


(37) 
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where 


m  n 

a 

m,  n 

f 

0 


Y 

M 


v 

s 


represents  the  amplitude  of  the  contribution  from  the  m'th 
pulse  of  the  processor  and  the  n'th  data  point 

standard  deviation  of  the  normalized  fractional  amplitude  on 
the  film 

spatial  frequency  of  the  data  in  cycles  per  millimeter 

phase  shift  from  pulse  to  pulse  due  to  doppler  measured  in 
radians 

distance  between  pulses  on  film  in  millimeters 
mean  velocity  of  the  film  drive 
mean  scan  rate  for  the  CRT 


Y  =  pulse  length  of  the  film,  50  millimeters 

<rR  *  standard  deviation  on  the  ropeatibility  of  the  Information 

from  pulse  to  pulse  on  the  film,  (Y <r„  *  <r  ) 

R  m,  n 


R  -  radius  of  the  capstan  in  the  camera 

<r  standard  dc  'lation  of  the  positional  error  of  the  film  drive, 

p  <%  V 


The  standard  deviation  of  the  rcpeatlbility  of  the  information  from  pulse  to  pulse  is 

dependent  on  the  repeatibility  of  the  CRT  and  the  errors  associated  with  the  film  sliding  back 
and  forth  in  the  y-dlrertlon. 


In  order  to  calculate  the  maximum  standard  deviation  <r,  of  the  sidolobos  of  the  envelope 

it  is  necessary  to  use  the  worst  possible  parameters  in  the  above  equation.  Typical  Inaccuracies 
for  this  device  are  already  known.  However,  it  is  not  known  how  these  "typical"  Inaccuracies 
are  related  to  the  standard  deviation  of  the  inaccuracy  of  the  device.  For  the  purposes  of 
this  calculation  it  Is  assumed  that  the  "typical"  inaccuracies  arc  equal  to  the  standard  devia* 
tion  of  the  inaccuracies  of  the  system.  This  means  that  the  Inaccuracies  are  less  than  the 
"typical"  inaccuracies  for  about  68. 31  of  the  time. 

Typical  inaccuracies  of  the  system  arc: 

(a)  repeatibility  of  the  CRT  of  1  part  in  10,  000 

(b)  instantaneous  positional  perturbation  of  the  film  drive  is  *1  second  of  arc 

(c)  intensity  deviation  on  the  CRT  is  5  percent 

(d)  the  film  will  shift  by  1  part  in  5  x  10®  in  the  y-dlrectlon  per  pulse. 

We  can  now  specify  some  of  the  other  parameters  associated  with  equation  (37): 

a  2  =  0.05 
m,  n 

f  +  T  =  a  maximum  of  20  cycles  per  millimeter 
0  =  a  maximum  of  n  radians 

Y  =50  millimeters 
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<rR  =  10- 

R  *  a  maximum  of  1/2  inch  or  12.7  mm 

(j-  *  4.  85  x  10* 6  radians 

P 

Some  additional  explanation  is  required  to  specify  the  other  parameters  in  this  equation. 
The  average  sweep  velocity  of  the  CRT  is  equal  to  vg  =  10G/32  N  x  10‘®  in  millimeters  per 

second.  This  was  obtained  from  the  fact  that  100  millimeters  of  the  tube  surface  will  be  used, 
there  will  be  N  data  points  displayed,  and  it  takes  32  microseconds  to  write  each  data  point  on 

-6 

the  film.  The  average  film  drive  velocity  is  vf  =  50/83  NM  x  10  .  This  was  obtained  as 

follows.  The  pulses  to  be  processed  occupy  50  millimeters  of  film,  it  takes  83  milliseconds 

from  the  start  of  one  pulse  to  the  start  o'  a  second  pulse  at  a  thousand  data  points  per  pulse, 

and  there  arc  M  pulses  in  the  50  millin<eters  of  film.  This  means  that  the  ratio  of  V.  to  V 

0  193  1  8 

is  '  ,  where  M  is  the  number  of  pulses  along  50  millimeters  of  film. 


Now  wc  can  investigate  the  factor  involving  amplitude  information  in  the  above  equation. 
It  is  assumed  that  the  amplitude  information  Am  n  can  be  separated  into  two  independent 

functions  a  dependent  only  on  m  and  a  dependent  only  oi  n.  Hence  A  >=  a  a  ,  and  the 
m  n  m,nmn 

above  equation  becomes : 


Now  we  can  define  the  factorl 


(ify  i 

'm=l  '  m*l 


This  factor  is  equivalent  to 


the  gain  of  an  antenna  array  normalized  with  respect  to  the  gain  of  a  uniformly  illuminated 
array.  The  optical  processor  will  probably  use  an  amplitude  taper  to  maintain  a  low  sidelobe 
level  in  the  doppler  dimension.  It  has  been  shown  by  D.  Sleplan3  that,  to  concentrate  this 
energy  into  a  region  of  the  film  such  that  the  sidelobe  energy  is  minimized,  it  is  necessary 
to  produce  an  amplitude  taper  proportional  to  a  prolate  spheroidal  wave  function  of  zero 
order.  This  same  criteria  was  also  used  for  a  discrete  20  element  array  and  it  is  felt  by 
the  author  that  the  results  are  approximately  these  that  were  obtained  by  Slepian.  Using  this 
work  the  variation  of  oc  with  respect  to  sidelobe  level  is  given  by  Figure  26-  Now,  using  this 
in  the  above  equation  and  the  fact  that  there  will  probably  be  no  taper  in  the  y -direction  (other 
than  that  already  included  in  the  radar  waveform)  we  have: 
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„  2  0.445  1.5  x  10'11  x  M 

^E  “cTKlff  *  oc“N - 

N  will  probably  take  on  the  value  of  2000.  Using  this  fact  and  simplifying  the  above 
equatt  an  we  have : 

„1  2.22  xlO*4  7.5xl0‘15xM 

E  r  OTM  *  a 

Since  M  will  range  between  500  and  1000  pulces,  it  is  seen  that  500  is  the  worst  case  number 
to  be  used.  This  yields: 

a  2  °- 444  x  10*6 
E  < x 

Now  it  is  possible  to  calculate  the  degradation  of  the  sidelobes  using  the  Rayleigh  dis¬ 
tribution  that  describes  the  behavior  of  the  sidelobes,  see  equation  (31).  Using  reference  5 
it  was  found  that: 

(Y- a)  =  3.75 

providing  the  sidclobe  level  of  the  envelope  stays  below  r  for  99. 90c  of  the  time  using  a  nominal 
sidolobe  level  of  a.  Solving  for  the  ratio  of  r/a  untng  the  above  two  equations  we  hive: 

r  i/n  x  10'3  +1 

a  a/ar 

Using  this  expression  it  is  now  possible  to  calculate  the  degradation  of  the  sidclobe  level  for 
various  nominal  sidelobe  levels.  Figure  27displays  this  degradation  as  a  function  of  nominal 
sidciobo  levels. 


Figure  27.  Sidelobes  Will  Do  below  The  Degraded  Sidelobe  Level  for  99. 9^  Of  The  Time 
For  A  System  Designed  To  Yield  A  Given  Nomii<al  Sidelobe  Level 
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2.  FORMAT  I 


Making  tic  same  assumption*  that  were  made  using  format  I,  and  Introducing  the  ap¬ 
propriate  positional  errors  in  equation  (3)  we  have 


M  N 


j)re«($  -  })[*  -  A,  Z  Z  <\,*Am  nAm) 

m»l  n*l 

(«*>  *'  <«  -  w  ■  \ 

+  m, n  ) 

|  [a0  •  .U>c !  (  f  -  *  -  T)J 


cos  (2r  f  ^  •m|l) 


sine  J  -  y )  cos  (2  *  f 


where 


»V)  »  (y  .  "Y  _  jW 
N'  ^  N  m,n'  |yfj 


5m  represents  the  error  in  positioning  the  m'th  pulse  in  x -direction  on  the 


film 


5  represents  the  error  in  positioning  the  n’th  data  point  of  the  m*th  pulses 
in  the  y-dlrcction  along  the  film 

Am  n  represents  the  fractional  variation  of  the  amplitude  Information  on  the  film 
Taking  the  two  dimensional  Fourier  transform  of  the  film,  we  have- 
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At  this  point  it  is  best  to  consider  only  one  of  these  targets,  keeping  in  mind  that  the  other 
target  may  be  added  using  the  principal  of  superposition.  The  one  sidelobe  is  now  passed 
through  another  aperture  located  at  U  =  f.  The  Information  at  the  eecond  aperture  can  be 
represented  by: 
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Taking  the  Fourier  transform  in  the  u  direction  we  have 

E  (v,  y’)  «  K  Jjilnc  xy*  exp  J2rfy'j  •  Maine  Yv  exp  (-J»Yv)l  • 

M  N  ^  . 
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The  terms  I(y’),  sine  Yv  and  sine  Xy’  act  like  impulses  which  have  no  great  effect  on  the 
system.  The  functions,  Pj(v),  lj(v)  and  rect  -  j)  act  like  wideband  filters  which  slightly 
attenuate  the  Information.  Using  these  facts  this  equation  can  be  simplified  to: 
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However,  the  functions  l(y’),  wine  Yv  and  sine  Xy*  cause  the  discrete  sampled  values 
to  run  together.  This  means  that  the  variation  in  the  rect|  — ^ - —  1  can  be  neglected 


vftr)  ™“ 


2.  hY' 

*  *  TT 


and  this  function  can  be  simplified  to  recti  -y^ I.  Rewriting  this  equation,  we  tave: 
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(40) 
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This  summation  In  it  represents  the  envelope  along  the  range  dimension.  This  envelope  has 
no  sidelobes  and  can  be  Ignored.  The  normalised  envelope  of  interest  is: 

I  A  H 


Currying  out  the  same  analysis  that  was  performed  on  equation  (18)  we  have: 
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We  are  Interested  In  this  equation  at  v  ■ 
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The  parameters  In  the  above  equation  arc  defined  in  equation  (37).  This  reduces  to: 


•2  =  VM  L?*05  +  M2  *  >-5  x  10'11  ♦  3. 


66  x  10’ 


Using  this  the  ratio  of  degraded  sidelobe  level  to  nomimU  sidelobe  level  is  given  by 

r  _  2. 65  x  10~2  ,  n 
a  uja~ 
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Calculated  values  of  degraded  side  lobe  level  vs.  nominal  level  are  shown  in  Figure  27, 
which  also  contains  the  curve  for  Format  II. 

IV.  CONCLUSIONS 

The  signal  power  in  the  random  error  component  sets  a  limit  on  attainable  side  lobe 
suppression,  with  given  system  parameters.  An  Inspection  of  the  curves  for  Formats  I  and  0 
(shown  in  Figure  27)  indicates  that  with  one -dimensional  integration  (Format  I)  the  suppression 
level  is  limited  at  about  30  dt>.  With  two-dimensional  Integration,  greater  suppression  of 
random  errors  is  possible,  with  45.5  db  capability  being  Indicated  by  the  curve.  It  is  antici¬ 
pated  that  the  poorer  sidelobes  in  Format  I  will  be  concentrated  in  a  ridge  (at  constant  range) 
across  the  output  plane.  This  is  because  the  error  power  before  integration  is  concentrated 
in  a  swath  across  the  output  plane  about  one  radar  range  cell  wide.  Format  n,  on  the  other 
hand,  has  ‘he  error  power  distributed  more  or  less  uniformly  over  the  input  plane  and  should 
not  have  ridge-like  doppler  sidelobes.  The  structure  with  Format  I  has  been  observed  with 
processed  dat  from  digital  computer  runs.  The  primary  error  sources  in  this  case  were  in 
radar  range  lobe  adjustments  and  in  phase  alignment  of  pulses  to  provide  a  coherent  doppler 
reference. 
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APPENDIX  m 


SOFTWARE 


1.  OPTICAL  PROCESSOR  PROGRAM 

This  is  the  main  program  In  preparing  radar  data  so  that  it  may  be  used  with  the  SDS  910 
and  the  film  recorder  to  produce  photographic  film  which  is  optically  processed.  The  input 
to  this  program  is  a  magnetic  tape  which  contains  radar  data  in  the  form  of  84  amplitude  and 
84  phase  quantities  per  pulse.  The  output  is  another  magnetic  tape  with  records  of  2048  or 
1024  words.  Each  record  represents  a  radar  pulse  that  has  been  interpolated  to  produce 
more  sample  points,  aligned  (alignment  is  not  necessary  with  simulated  radar  data  proAiced 
by  DATA2X),  applied  to  a  spatial  carrier  with  phase,  and  then  added  to  a  bias  level.  This 
output  tape  is  in  a  format  compatible  with  the  SDS  910. 

a.  Logic 

This  program  begins  by  rewinding  FORTRAN  tape  units  11  and  12.  It  then  sets  up 
the  sin  X/X  table  which  Is  ueed  to  Interpolate  the  radar  data.  In  this  table  X  varies  from  0 
to  85.  5  by  Increments  of  »/10.  The  program  now  reads  all  the  data  cards  (cards  1  thru  6). 

It  also  calls  subroutines  LSFIT4.  This  subroutine  sets  up  a  look  up  table  to  compensate  for 
the  overall  film  recorder  nonilnearities.  Then  a  table  Is  set  up  for  the  spatial  carrier  and 
all  the  variables  needed  in  calling  any  of  the  alignment  subroutines  are  defined.  Following 
this  is  the  logic  for  Input  tape  control.  Here,  according  to  the  values  of  the  variables  on  the 
dkta  cards,  the  input  tape  can  be  advanced  a  number  of  pulses,  backspaced,  rewound,  or 
given  number  of  pulses  can  be  skipped  between  each  processed  pulse.  After  a  pulse  is  read 
in  from  tape,  one  of  the  three  alignment  subroutines  (AL1GN1,  AL1GN2,  AL1GN3)  is  called 
depending  on  the  value  of  NAL1GN,  an  Input  variable.  The  align  subroutines  convert  the  84 
amplitudes  and  84  phases  to  real  and  imaginary  parts.  These  are  then  interpolated  using 
sin  X/X  to  give  ten  points  between  every  two  points  of  the  original  84.  Th's  results  in  831 
real,  and  831  imaginary  values  per  pulse.  Then  the  pulses  are  aligned.  The  method  used 
depends  on  which  align  subroutine  is  called.  Simulated  radar  data,  produced  by  DATA2X, 
is  not  aligned. 

The  data  now  is  in  aligned  and  partially  Interpolated  form.  Linear  interpolation  be¬ 
tween  adjacent  points  is  now  used  to  form  the  necessary  1024  or  2048  points  per  radar  pulse. 
Now  the  pulse  number  and  the  maximum  amplitude  for  the  present  frame  (a  frame  is  a  set  of 
pulses  which  are  recorded  adjacent  on  film)  is  printed.  The  maximum  amplitude  for  a  frame 
was  found  by  subroutine  AMPMAX  as  the  data  was  generated.  II  A  (n)  represents  the  n,h 
amplitude  sample  point  of  a  pulse  and  0(n)  represents  the  n,h  phase  sample  point  (there  are 
1024  or  2048  sample  points  per  pulse)  and  AMPMAX  is  the  maximum  amplitude  of  a  given 
frame,  then  need  to  form  the  following  function: 

AMOD(n)  -  AMPMAX  ♦  A(n)  cos  (2»fcn  ♦  0(n) 

where  n  varies  from  1  to  1024  or  2048,  fc  is  the  spatial  carrier  freijiency.  Since  the  data  is 
not  now  in  terms  of  amplitude  and  phase,  but  in  the  form  of  real  and  imaginary  values,  much 
computation  time  can  be  saved  by  not  converting  to  amplitude  and  phases  but  by  rewritting 
(1)  as; 

AMOD(n)  *  AMPMAX  ♦  A(n)  [cos  (2sfcn)  cos  0(n)  -  sin  (2ifcn)  sin  0(n)J 
or 

AMOD(n)  =  AMPMAX  ♦  cos  (2*fcn)  ^A(n)  cos  0  (n)  J  -  sin  (2»fcn)  [A(n)  sin  0(n)J 
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where  the  terms  in  brackets  A(n)  cos  0(n)  and  A  (n)  sin  0(n)  are  the  real  and  imaginary  values. 
The  ^lantities  cos  (2?fcn)  and  sin  (2xfcn)  are  the  spatial  carrier  terms  which  are  given  in  a 
look-up  table  formed  In  the  beginning  of  this  program,  and  which  do  not  change  from  pulse  to 
pulse  in  a  given  frame. 


Jh*  values  of  AM0D<n)  are  now  fixed,  and  scaled  so  that  they  are  within  the  values 
-255  to  *256,  i.  e.  9  bits.  If  an  end  of  frame  is  now  reached  an  end  of  file  is  written  on  the 
output  tape.  Here  subroutine  PACK  is  called.  This  subroutine  checks  the  input  data  and 
makes  certain  each  value  is  within  -255  and  *256.  It  also  packs  the  data  on  the  output  magnetic 
tape  in  a  format  compatible  with  the  SDS  910.  Now  the  total  time  of  execution  of  this  program 
is  written  out,  and  end  of  file  is  written  on  the  output  tape,  and  the  program  pauses  before 
exit  so  that  tapes  may  be  removed.  See  page  97  for  flow  chart. 

b.  Subroutines  for  Optical  Processor  Program 
(1)  ALIGNI  -  Prominent  Point  Align  Subroutine 

The  prominent  Point  Align  Subroutine  aligns  the  radar  pulse  returns  by  making 
the  phase  at  the  prominent  point  zero.  All  other  phase  values  are  altered  by  the  amount  of 
this  prominent  point  phase. 

Logic 

The  routine  determines  a  number  which  is  the  last  desired  channel  to  be  processed  The 
routine  multiplies  the  amplitude  channel  value  by  the  cosine  of  the  corresponding  phase  chan¬ 
nel,  and  also  by  the  sine  of  the  phase  channel.  These  two  values  arc  the  real  and  the  imaginary 
portions  of  the  complex  return.  The  above  is  repeated  with  each  desired  channel  The  re¬ 
sults  of  the  above  is  a  conversion  from  amplitude  and  phase  to  real  imaginary  arrays. 

The  routine  takes  these  two  arrays  and  fills-ln  between  the  sampling  points.  The  inter¬ 
polating  is  with  SIN  X  X  values  calculated  in  the  main  program.  For  each  position,  the 
routine  multiplies  each  sample  value  by  the  correct  value  of  SINXA:  and  sums  all  these  pro¬ 
ducts  to  give  the  correct  weight  at  that  position.  The  routine  goes  through  the  above  with 
each  position  to  give  new  real  and  Imaginary  arrays.  These  arrays  are  the  arrays  with  nine 
vatues  between  each  of  the  old  values.  For  each  desired  position  in  the  above  arrays,  the 
routine  squares  the  real  and  the  imaginary  values;  and  adds  these  together.  This  is  the  squar® 
of  the  amplitude  of  the  return  at  this  position.  Each  squared  amplitude  over  the  desired  re¬ 
gion  is  calculated.  The  maximum  value  and  the  position  of  this  value  is  found.  The  phase  at 
this  amplitude  is  subtracted  from  the  phase  of  each  position.  The  real  and  imaginary  portions 
with  the  adjusted  phases  are  determined.  This  routine  outputs  these  arrays  and  the  position 
of  the  maximum  point  in  the  arrays.  For  simulated  radar  data  produced  by  DATA2X  align¬ 
ment  is  not  necessary.  ’  * 

The  calling  for  this  routine  is 


CALL  ALIGNI  (ICH,  DATA,  R,  F,  SXDX,  NFA,  N0C,  MID,  NDEL,  INITL) 
where  the  inputs  are: 

ICH  -  number  of  amplitude  and  phase  channels  in  the  pulse  of  data. 

DATA  -  name  of  data  array  which  is  half  amplitudes  and  half  phases. 

SXDX  *  ,he  sir.  XOC  array  used  for  interpolation. 


NFA 

N0C 

MID 


number  of  the  first  desired  channel  of  amplitude  data  to  be  used. 

number  of  amplitude  channels  to  be  used. 

number  of  ’he  middle  amplitude  channel  in  the  pulse. 
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NDEL  -  number  of  positions  In  the  interpolated  amplitude  array  to  be  searched  each  side 
of  the  middle  position  for  the  maximum  amplitude. 

INITL  -  control  constant  which  is  zero  for  the  first  pulse  and  non-zero  for  the  rest  of 
the  pulses. 


and  the  outputs  are: 


-  The  complex  radar  return  is  converted  to  real  and  imaginary  arrays,  then  these 
arrays  are  aligned.  R  is  this  aligned  real  array  for  the  pulse. 

-  The  Imaginary  aligned  array  for  the  pulse. 


MID 


-  The  position  number  in  the  interpolated  array  of  the  maximum  amplitude. 
(2)  ALIGN2  -  Cross-correlation  Align  Subroutine 


ALIGN2  aligns  the  radar  pulse  returns  by  cross -correlating  each  return  with  a 
reference  return.  The  maximum  correlation  is  determined.  The  phase  difference  of  the  two 
returns  at  this  maximum  correlation  is  subtracted  irom  each  phase  value  In  the  return  to  be 
aligned.  Also,  the  delay  between  the  two  returns  at  the  maximum  correlation  Is  subtracted 
from  the  return  to  be  aligned. 


Logic 

The  program  defines  the  midpoint  of  the  arrays  which  are  outputs  from  this  routine. 
These  arrays  are  the  real  and  imaginary  returns  sampled  at  every  0.  2  nanoseconds.  The 
routine  defines  the  number  of  the  last  desired  channel  to  be  processed,  and  then  starts  pro¬ 
cessing  with  i he  first  desired  channel.  The  cosine  of  the  phase  channel  is  multiplied  by  the 
value  in  the  corresponding  amplitude  channel,  and  also  the  sine  of  the  phase  channel  is  multi¬ 
plied  by  the  value  «n  the  corresponding  amplitude  channel.  These  two  products  are  the  real 
and  the  im  binary  portions  of  the  complex  radar  return.  Each  amplitude  -  phase  channel  is 
converted  to  a  real  and  an  imaginary  value  in  the  above  way. 

The  routine  fllls-tn  the  above  arrays  using  sin  X/X  Interpolation.  This  interpolation 
converts  the  above  arrays  which  are  sampled  at  every  2.0  nanoseconds  to  arrays  which  are 
sampled  at  every  0.2  nanoseconds.  This  fill-in  is  accomplished  for  a  point  by  summing  the 
products  of  each  2.  0  nanosecond  sample  and  a  sin  X/X  value.  The  sin  X/X  value  used  with 
a  sample  represents  the  distance  of  the  cample  from  the  point  to  be  filled  In.  The  sin  X/X 
array  has  a  zero  crossover  at  every  tenth  position,  thus,  the  2.  0  nanosecond  samples  which 
arc  at  every  tenth  position  on  the  new  position  scale  are  uncharged  by  this  Interpolation. 

After  this  interpolation,  the  routine  tests  a  quantity,  which  was  set  equal  to  zero  in  the 
main  program,  for  being  zero.  For  the  first  pulse  this  quantity  is  zero,  thus  the  routine 
multiplies  each  element  in  the  interpolation  arrays  oi  the  first  pulse  return  by  the  card  Input 
B.  This  modified  return  becomes  the  reference  with  which  the  second  return  is  cross- 
correlated.  After  this  first  reference  is  established,  the  routine  defines  the  quantity  which 
had  been  zero  to  be  equal  to  77,  and  then  the  routine  returns  to  the  main  program.  For  all 
pulse  returns  after  the  first  to  be  used  in  a  2D  map,  the  test  on  the  above  quantity  allows  the 
logic  flow  that  gives  the  cross -correlation  of  the  pulse  return  with  the  reference. 

When  two  pulse  returns  differ  only  in  their  relative  range  and  overall  phase,  cross- 
correlation  may  be  used  to  find  both  (fianttties.  If  one  function  is  shifted  along  the  time  axis 
while  the  other  function  is  fixed,  a  set  of  correlation  values  is  determined.  If  two  functions 
have  the  same  wave  form  the  maximum  correlation  will  occur  when  the  two  functions  have 
no  relative  delay.  When  the  two  functions  are  complex,  the  complex  conjugate  of  one  function 
must  be  used  in  the  cross-correlation  integral.  In  this  routine  the  complex  return  to  be 
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aligned  Is  cross-co“related  with  the  complex  reference  return.  The  routine  defines  the  num¬ 
ber  of  positions  on  the  0.  2  nanosecond  sampling  scale  that  the  return  to  be  aligned  Is  to  be 
shifted  relative  to  the  reference  return.  The  shift  that  gives  the  maximum  correlation  re¬ 
presents  the  relative  delay  between  the  two  returns.  The  routine  detei  mines  the  cross¬ 
correlation  for  each  shifted  position  and  compares  the  correlation  values  to  find  the  maximum. 
The  shift  that  gave  the  maximum  correlation,  and  the  real  and  imaginary  values  of  the  cor¬ 
relation  sumn  ation  at  the  maximum  correlation  are  used  to  align  the  return.  Each  real  and 
Imaginary  valie  of  the  return  is  shifted  by  the  relative  delay  determined  in  the  above  cross- 
correlation.  The  real  and  imaginary  values  of  the  above  cross -correlation  summation  at  the 
maximum  correlation  are  converted  to  a  phase  which  is  the  overall  phase  difference  between 
the  returns  at  the  maximum  ccrrelatlon.  This  phase  difference  Is  subtracted  from  each  phase 
of  the  return  to  be  aligned.  Ti  e  result  of  all  the  above  is  to  subtract  out  the  range  and  gross 
phase  changes  between  returns  The  routine  defines  the  real  and  the  imaginary  arrays  re¬ 
presenting  this  aligned  return.  The  routine  also  forms  a  new  reference  return  as:  one  minus 
B  times  the  old  reference,  pl-^a,  B  times  the  just  aligned  return. 

The  routine  goes  back  to  the  main  program  with  this  aligned  pulse  return,  and  leaving  an 
updated  reference.  No  alignment  is  necessary  for  simulated  data  produced  by  DATA2X.  The 
calling  sequence  for  this  routine  is 

CALL  AL1GN2  (1CH,  DATA,  R,  F,  SXDX,  NFA,  NOC,  LLMIN,  LLMAX,  1N1TL, 

MID,  B) 


where  the  Inputs  are: 

ICH  -  number  of  amplitude  and  phase  channels  In  the  pulse  of  data 

HATA  -  name  of  data  ray  which  is  half  amplitudes  and  half  phases 

-  the  sin  X/X  array  used  for  interpolation 

-  Number  of  the  first  desired  channel  of  an.  >! i t u  Je  data  to  be  used 

-  nMnt>er  of  amplitude  channels  to  tx?  used 


SXDX 

NFA 

NOC 


LLMIN  -  number  of  positions  below  center  to  be  shifted  for  correlation  values 
LLMAX  -  numlier  of  positions  atwvc  center 

1NITL  -  control  constant  which  is  zero  for  the  first  pulse  and  non-zero  for  the  rest  of 
pulses  making  up  a  2D  map 

MID  -  number  of  the  middle  amplitude  channel  in  the  radar  pulse 
B  *  smoothing  constant  used  on  reference  pulse 
and  the  outputs  are: 


-  The  complex  radar  return  is  converted  to  real  and  imaginary  arrays,  then  these 
arrays  are  aligned.  R  is  this  aligned  real  array  for  the  pulse. 

-  the  imaginary  aligned  array  for  the  pulse 

-  the  position  number  in  the  interpolation  array  of  the  maximum  amplitude. 
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(3)  ALIGNS  -  Smooth  Track  Align  Subroutine 

ALIGNS  allgas  the  radar  returns  by  calculating  from  orbital  parameters  th* 
range  and  phase  that  the  return  should  have.  The  first  return  is  used  as  a  range  reference 
and  all  other  returns  are  corrected  relative  to  this  first  return.  The  r»n-llnear  portion  of 
this  calculated  equation  is  used  to  find  the  phase  to  be  subtracted  from  each  phase  In  the 
return. 

Logic 

The  program  defines  the  midpoint  of  the  arrays  which  are  outputs  from  this  routine. 
These  arrays  are  the  real  and  Imaginary  returns  sampled  at  every  0.  2  nanoseconds.  The 
routine  defines  the  number  of  the  last  desired  channel  to  be  processed  and  then  starts  oroces- 
sing  with  the  first  desired  channel.  The  cosine  of  the  phase  channel  is  multiplied  by  the 
\alue  in  the  corresponding  amplitude  channel,  and  also  the  sine  of  the  phase  channel  is  multi¬ 
plied  by  the  value  In  the  corresponding  amplitude  channel.  These  two  products  are  the  real 
and  the  imaginary  portions  of  the  complex  radar  return.  Each  amplitude-phase  channel  is 
converted  to  a  real  and  an  Imaginary  value  in  the  above  way. 

The  routine  fills  -in  the  above  arrays  using  sin  (X)/X  interpolations.  This  interpolation 
converts  the  above  arrays  which  are  sampled  at  every  2.  0  nanoseconds  to  arrays  which  are 
sampled  at  every  0.  2  nanoseconds.  This  fill-in  is  accomplished  for  a  point  by  summing  the 
products  of  each  2.  0  nanosecond  sample  and  a  sin  (X)/X  value.  The  sin  X/X  value  used  with 
a  sample  represents  the  distance  of  the  sample  from  the  point  to  be  filled  in.  The  sin 
(X)  X  array  has  a  zero  crossover  at  every  tenth  position,  thus,  the  2.  0  nanosecond  samples 
which  are  at  every  tenth  position  on  the  new  position  scale  are  unchanged  by  this  interpolation. 

After  this  interpolation,  the  routine  tests  a  quantity,  which  was  set  equal  to  zero  in  the 
main  program,  for  being  zero.  For  the  first  pulse  return  this  quantity  is  zero,  thus  the 
routine  follows  the  logic  which  establishes  the  reference  range  quantity  to  be  used  with  all  the 
remaining  pulse  returns.  This  range  quantify  ’DELTA"  is  the  recorded  "flnerange"  minus 
the  calculated  time  dependent  range.  Thl3  calculated  time  dependent  range  has  the  following 
equation: 

R(t)  *  Ajt  ♦  Ajt  ♦  Ajt^  ♦  A^t*  ♦  A,.!5  *  Agt6 

Where  the  A's  arc  the  inputs  from  the  main  program,  and  t  is  the  relative  time  M:  >  pulse  was 
predicted  to  reach  the  Jarget.  The  first  term  (A.t)  is  subtracted  from  the  above  range  equa¬ 
tion  leaving  only  the  non-linear  terms.  This  resultant  non-linear  value  is  converted  into 
radians  which  represent  the  phase  the  return  should  have.  The  routine  subtracts  this  phase 
from  each  value  of  phase  in  the  0.  2  nanosecond  sample,  and  returns  to  the  main  program  with 
the  above  reference  and  with  the  real  and  imaginary  phase  corrected  radar  return. 

For  the  second  and  all  remaining  radar  returns  to  make-up  a  2D  map,  the  routine  follows 
the  above  logic  until  after  the  test  on  the  above  quantity  which  had  been  zero  for  the  first  radar 
return.  After  the  first  radar  return,  this  quantity  was  set  equal  to  +77,  thus,  the  test  results 
in  the  routine  following  the  logic  which  uses  the  range  reference  quantity  DELTA  with  the 
radar  return  to  be  aligned.  The  routine  determines  the  new  time  t,  and  uses  it  in  the  time 
dependent  range  equation  in  the  same  way  as  for  the  first  pulse.  The  routine  now  determines 
the  range  shift  (time  difference)  between  the  first  pulse  and  the  one  to  be  aligned.  The  dif¬ 
ference  between  the  pulse  to  be  aligned  the  first  in  nanoseconds  is  equal  to  the  time  dependent 
range  for  the  pulse  plus  DELTA,  minus  the  fine  range  of  the  pulse,  all  divided  by  ten  times 
the  speed  of  light  in  meters  per  nanosecond. 

DR  =  (R  (t)  *  DELTA  -  FRANGE)/10C 

The  number  of  positions  the  0.  2  nanosecond  sampled  pulse  must  be  shifted  to  align  it  with  the 
first  is  ten  times  DR.  The  phase  of  the  pulse  is  determined  In  the  same  way  as  it  was  found 
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10*Dff  return>  The  real  and  the  lm*S*nary  samples  for  the  pulse  are  shifted 

-a,  hRh«  .  *!  k  al  *"  the  pu,8es  ln  range*  The  calculated  phase  is  subtracted  from 

viues  Thr^uH3  "““‘T  a,igned  radar  re,urn  converted  to  real  and  Imanna-y 

CALL  AUGN3  (1CK,  DATA,  R,  F,  SXDX,  NFA,  NOC,  fNITL,  MID,  FRUP,  FRLO, 
NL,  Al,  A2,  A3,  A4,  A5,  A6.  OMC,  TTS,  NPULSE,  VEL) 

where  the  inputs  are: 

—  -  number  of  amplitude  and  phase  channels  in  the  pulse  of  data 
DATA  -  name  of  data  array  which  Is  half  amplitudes  and  half  phases 
— D—  *  ,he  sin  X  X  array  used  for  Interpolation 

NFA_  .  number  of  the  first  desired  channel  of  amplitude  data  to  be  used 
NOC_  -  number  of  amplitude  channels  to  be  used 

—  ’  Cp°2lf  con•s,an,  which  is  2er°  *or  the  first  pulse  and  non-aero  for  res.  of  the 

MID  -  number  of  the  middle  amplitude  channel  In  the  radar  pulse 

FRUP  -  the  16  most  significant  bits  of  the  40  bit  ’’fine  range”  number 

FRL0  •  the  24  lower  order  bits  of  "fine  range” 

NL_  -  number  of  pulses 

Al  -  modified  first  range  derivative 

A2  .  modified  second  range  derivative 

A3  -  modified  third  range  derivative 

A4  -  modified  fourth  range  derivative 

A5  -  modified  fifth  range  derivative 

A6  -  modified  sixth  range  derivative 

°MC  -  center  frequency  of  radar 

TTS  -  time  between  processed  pulses 

NPULSE  -  number  of  pulse 
VEL  -  speed  of  light 

and  tne  ot  puts  are: 

—  ‘  *  li'*  comP,ex  radar  return  is  converted  to  real  and  imaginary  arrays  then  these 

arrays  are  aligned.  R  Is  this  aligned  real  array  for  the  pufse 
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F_  -  the  Imaginary  aligned  array  for  the  pulse 

MID  -  the  position  number  In  the  interpolated  array  of  the  maximum  amplitude. 

(4)  LSFIT4  -  Least  Squares  Fit  Subroutine 

This  subroutine  accepts  N  points  defined  by  X(I),  Y(I).  It  then  fits  a  fourth  order 
pjlynominal  to  these  points  by  the  "method  of  least  squares".  The  output  of  this  subroutine 
are  the  five  coefficients  that  define  the  polynomtnal: 

Y(x)  -  A  ♦  Bx  ♦  Cx2  ♦  Dx3  ♦  Ex4  (2) 

This  function  y(x)  will  give  a  better  approximation  to  the  original  points  X(I),  Y(I)  as  N  in¬ 
creases:  l.  e.  the  subroutine 's  given  more  points  in  a  given  interval. 

It  is  with  this  function  y(x)  that  the  compensation  look  up  table  is  formed.  The 
inputs  to  the  table  are  X  values,  the  outputs  are  Y  values.  Therefore,  it  is  necessary  to  choose 
the  X(l),  Y(l)  so  that  the  amplitude  values  after  going  through  the  compensation  look  up  table, 
when  applied  to  the  grid  of  the  CRT  in  the  film  recorder,  result  in  a  film  whose  amplitude 
transmission  is  linearly  proportional  to  the  uncompensated  amplitude  values. 

The  calling  sequence  for  this  routine  is: 

CALL  LSFIT4  (X,  Y,  N,  A.  B,  C,  D,  E) 

where  he  inputs  are; 

X,  Y  -  These  are  two  arrays  defined  »n  the  main  program.  They  are  read  in  from  input 
cards  3  on,  depending  on  N,  i.  e.  if  N  is  greater  than  4  .norc  than  one  card  is  re¬ 
quired. 

N  -  This  Is  the  numt>cr  of  points  (defined  by  X.  Y  sets)  to  which  a  curve  is  fitted, 

the  outputs  arc; 

A,  B,  C,  D,  E  *  These  are  the  coefficients  of  the  fourth  order  polynonnnal  as  defined  by 
Equation  (2). 

(5)  S1MEQ  -  Sul) routine  for  Solution  of  Simultaneous  Equations 

This  subroutine  Is  used  by  LSFIT4  In  solving  a  set  of  simultaneous  equations 
needed  in  determining  the  coefficients  A.  B,C,  D,  E.  by  the  method  of  least  squares. 

The  calling  sequence  for  this  routine  is: 

CALLSIMEQ(A,  NROWA,  NCOLB,  DET,  MAXKOW) 

where  the  inputs  arc; 

A  -  a  matrix  array  containing  the  B  vectors 

NROWA  -  number  of  rows  in  matrix  A 

NCOLB  -  number  of  B  vectors 

MAX  RON  -  maximum  row  dimension  of  A  array. 
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the  outputs  are; 


DET  -  value  of  determinate  of  A  matrix 

(The  solutions  to  a  set  of  equations  are  placed  In  first  NCOLB  columns  of  matrix  A). 

(6)  PACK  -  Subroutine  Which  Writes  Tape  for  SDS  910 

The  UNIVAC  1108  uses  a  36  bit  word.  The  SDS  requires  the  1024  or  2048  word 
arrays  to  be  written  on  tape  such  that  each  word  consists  of  12  bits.  These  words  arc  adja¬ 
cent  on  the  tape  with  no  gaps  except  a  record  gap  every  1024  or  2048  words.  The  first  3  bits 
of  the  12  bit  word  are  not  used.  The  next  9  bits  are  used  to  represent  amplitude  values  rang¬ 
ing  from  -255  to  256.  Positive  numbers  are  represented  in  the  usual  way.  A  positive  number 
becomes  negative  by  taking  a  two’s  compliment  of  the  9  bits. 

This  subroutine  accepts  the  36  bit  UNIVAC  1108  word  and  packs  it  onto  magnetic 
tape  as  described. 

The  calling  sequence  for  this  routine  is; 

CALL  PACK  (IW1,  L) 
the  inputs  arc; 

1W1  -  arrays  containing  values  to  be  written  on  magnetic  tape  for  SDS  910. 

L  -  size  of  IWI  array.  After  each  L  numbers  a  record  gap  is  put  on  tape. 

(The  output  is  a  magnetic  tape). 


4  Optical  Processor  Program  Data  Cards 


Card  1.  Format  8F10.  2 


Card  2.  Format  1415 


Each  number  Is  the  maximum  amplitude  lor  a  given 
frame.  This  maxlnvim  for  each  frame  Is  obtained  from 
AMPMAX  subroutine,  not  part  of  this  program.  AMX(l) 
Is  the  maximum  for  the  first  frame,  AMX(2)  the  second 
frame,  etc. 


NIPP 


This  is  the  number  of  interpolations  per  pulse  and  must 
be  either  1024  or  2048. 


NPPF 


Total  number  of  points  X,  Y  to  be  read  in  on  the  next 
set  of  data  cards  (number  of  cards  depends  on  N). 

Number  of  pulses  per  frame.  This  Is  fixed  for  a  given 
run. 


Card  3.  Format  8F10.  3 
X(l).  Y(l)  to 

xM:  mr 


Card  4.  Format  5E14. 8 


These  arc  points  which  are  used  by  the  program  to  form 
a  look  up  table  to  compensate  for  the  nonttnearttles  of 
the  film  recorder.  The  explanation  of  LSFIT4  defines 
how  these  points  are  detei  mined.  If  N  is  greater  than 
4.  more  than  one  card  is  required. 


Th's  Is  the  smoothing  constant  for  the  cross-correla¬ 
tion  alignment  subroutine  (AL1GN2)  where  each  pulse 
Is  aligned  to  give  maximum  correlation  wl'.h 


where  Pn  is  the  immediately  preceding  aligned  pulse 
and  Pn  is  the  weighted  average  of  all  aligned  preceding 
pulses. 


This  is  the  center  frequency  of  the  radar  in  MHz. 


FC  -  This  is  a  quantity  that  determines  the  spatial  carrier. 

For  10  line  pair  per  millimeter  FC  *  .293. 

Card  S.  Format  1415 

ITAP  -  This  is  the  data  tape  rewind  control.  If  ITAP  /  0.  the 

data  tape  on  Fortran  IV  tape  unit  11  is  rewound  If 
ITAP  -  0  it  is  not  rewound. 

LF  -  This  is  the  number  of  pulses  skipped  on  the  data  tape 

before  the  first  one  used. 
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MODUL 


-  If  this  pulse  skipping  control  equals  N,  every  nth  pulse 
is  skipped.  If  MODUL  «  I,  every  pulse  is  pr^~essed. 


NALIGN 


NFC 

NL 

NOC 

Card  6.  Format  6F12. 6 


-  This  determines  which  align  subroutine  Is  to  be  used. 
Simulated  radar  data  produced  by  DATA2X  is  not 
aligned.  If  NALIGN  •  1,  ALIGN!  is  used.  If 
NALIGN  -  2,  ALIGN2  Is  used.  If  NALIGN*3  ALICN3 
is  used. 

-  This  is  the  number  of  the  first  channel  of  the  84  to  be 
processed.  Normally  NFC  *  1. 

•  This  Is  the  total  number  of  pulses  to  be  processed  by 
given  run. 

■  This  is  the  number  of  channels,  starting  at  NFC. 
actually  used  for  each  pulse. 


RH01,  RM02 . RHOG 


These  six  numbers  are  the  coefficients  for  a 
power  series  representing  range  versus  time 
used  in  ALIGN3. 


2.  DATA2X  -  RADAR  DATA  SIMULATION  PROGRAM 

This  program  produces  a  data  tape  in  the  calibrated  radar  data  format.  This  tape  can 
be  used  to  test  the  optical  recording  system  and  software.  The  program  simulates  the  radar 
returns  from  a  rigid  rotating  system  of  point  scatters.  At  the  present  time  it  is  possible  to 
simulate  up  to  twenty  point  scattcrers.  Each  point  scatterer  can  have  a  arbitrary  amplitude 
position  and  relative  phase.  The  rigid  system  of  point  scattcrers  can  then  be  rotated  at  any 
desired  rotation  rate  and  recorded  at  a  given  pulse  repetition  frequency.  There  arc  two  modes 
of  operation  of  this  program.  In  the  first  mode  the  return  of  each  point  scatterer  of  a  given 
amplitude  is  of  tin?  form  sin  x2_c2  •  whcrc  c  is  chosen  to  give  a  desired  peak  to 

side  lobe  ratio.  This  mode  is  designated  as  Format  I.  In  the  second  mode,  the  return  from 
each  point  scatterer  is  the  Fourier  transform  in  the  range  dimension  of  sin  f "5  2 //  2  2 

This  mode  produces  Format  II  data.  \x  *C  Al  x  -C 


a.  Logic 

The  program  begins  by  calling  sub  routine  SX0XS.  This  subroutine  sets  up  a  sin  x/x 
table  which  is  used  by  subroutine  AMPMAX  to  find  the  maximum  amplitude  of  a  pulse.  Now 
rata  cards  1  and  2  arc  read  in.  Here  the  variable  M FORMAT  Is  tested.  If  M FORMAT  •!, 
Format  I  data  Is  produced I.  Thlsj?  done  by  calling  subroutine  RZERO.  This  subroutine  sets 
up  a  a  >lc  of  sin  _  <3-  which  is  used  in  producing  Format  I  data.  If  MFORMAT=2, 

subroutine  BAMPjs  cajled^  This  produces  a  table  which  Is  the  Fourier  Transform  of  the 
tuncllon  sln^x2.cyv/x2.c2  in  the  range  dimension.  This  Is  then  used  to  produce  Format  II 

data.  For  this  Format  NSORT  is  set  equal  to  2.  this  determines  how  subroutines  AMPMAX 
locates  the  maximum  amplitude  of  a  pulse.  Now  data  card  3  is  read  in  and  various  parameters 
are  defined.  The  remaining  data  cards  arc  now  read.  The  number  of  cards  depends  on  the 
number  of  point  scattcrers  and  on  the  number  of  different  target  configurations. 

Each  pulse  of  data  consists  of  84  amplitude  and  84  phase  channels.  Therefore  It  is 
necessary  to  calculate  the  radar  response  of  the  point  scatterer  configuration  at  84  points. 

This  is  done  by  assuming  that  the  radar  return  from  a  point  scatterer  is  given  by  the  table 
computed  by  subroutine  RZEROM  for  Format  I.  or  by  the  table  from  subroutine  BAMP  for 
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To.'nul  n.  In  range,  the  axis  of  rotation  of  the  point  scatters  is  assumed  to  be  at  the  center 
of  the  84  range  channels.  The  program  now  computes  the  relative  phase  and  range  of  each 
point  scatterer  The  function  that  was  calculated  by  either  RZEROM  or  BAMP  is  now  multi¬ 
plied  by  the  amplitude  of  each  point  scatterer  and  the  contribution  from  all  the  point  scatterers 
is  summed  at  every  channel.  This  is  done  for  both  the  real  and  imaginary  components.  Now 
the  data  is  converted  to  amplitude  and  phase  at  each  channel  and  written  onto  magnetic  tape. 
Subroutine  AMPMAX  is  called  to  find  the  maximum  amplitude  of  this  pulse.  The  target  is  now 
rotated  through  an  angle,  which  depends  on  the  pulse  repetition  frequency  and  the  rotation 
rate,  and  the  same  computation  Is  begun  for  a  new  pulse. 

After  ali  the  puises  for  a  given  target  are  completed,  the  maximum  amplitude  of  these 
puises  is  found  from  the  maximum  of  each  pulse.  The  program  now  reads  In  the  value  of  IEND 
on  the  last  card  for  a  given  target.  If  LEND  ■  0  a  new  set  of  target  configuration  data  cards 
is  read  in.  If  IEND  /  0  an  END  OF  FILE  is  placed  on  the  tape  and  EXIT  is  called.  See  page  98 
for  flow  chart. 

b.  Subroutines  for  DATA2X  Program 

(1)  RZFROM  -  Subroutine  for  Format  1  response  Pattern. 

This  routine  sets  up  a  table  of  the  function  sin  The  table  is 

the  array  RZ.  RZ(I)  is  the  value  of  the  above  function  at  X  >0.  Then  X  is  incremented  by 
the  amount  DTAU  until  there  are  855  values  of  RZ.  The  peak  amplitude  of  the  above  function 
is  (ec  -  c*c)/2c  and  the  first  zero  crossing  is  at  ^ +  ^,2  .  The  callir^  sequence  is: 

CALL  RZEROM  (K,  DTAU.  C,  OMB,  TAU,  RZ,  IPT) 


where  the  inputs  are: 


K  -  number  of  sample  points  is  K  ♦  1. 

DTAU  -  increment  of  X. 


C  -  determines  the  peak  to  side  lobe  ratio  of  the  response  pattern. 

OMB  -  This  is  the  scale  factor  of  the  X  values.  With  this  scale  factor  the 

new  variable  U  »  (OMB)*X/2.  Now  the  first  zero  crossing  occurs 

**  X  -  (2*  fT?  )/OMB. 

IPT  -  Print  Control.  If  IPT  ■  0  there  is  no  print  out.  If  IPT  /  0  then 

the  TAU  and  RZ  arrays  are  printed  out. 

and  the  outputs  are: 


TAU  •  an  array  containing  the  K  ♦  I  X  values. 

RZ  -  an  array  containing  the  K  ♦  1  values  of  the  Format  I  response 

pattern. 

(2)  BAMP  -  Subroutine  for  Format  II  Response  Pattern. 


This  routine  uses  the  modified  Bessel  function  of  order  zero  to  calculate  a  table 
B  which  is  an  array  containing  the  Fourier  Transform  of  the  table  formed  by  subroutine 
RZEROM.  This  table  is  then  used  by  DATA2X  to  produce  Format  II  data.  The  calling  se¬ 
quence  is 


CALL  BAMP  (DDD,  LLL,  B,  IPT) 
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where  the  Inputs  are: 

-  the  side  lobe  suppression  constant. 

-  number  of  samples  In  the  B  array. 

-  print  control,  if  IPT  *  0  there  ts  no  print  out.  If  IPT  4  0  the  B 
array  Is  printed  out. 


-  an  array  containing  the  LLL  values  of  the  Format  II  response 
pattern. 

(3)  AMPMAX  -  Subroutine  For  Finding  The  Maximum  Amplitude  of  a  Pulse 

Since  radar  amplitudes  for  a  pulse  are  sampled  at  only  84  points,  the  maximum 
of  these  may  not  be  the  true  maximum.  Interpolation  of  all  84  points  to  find  the  maximum  Is 
time  consuming.  This  routine,  depending  on  the  variable  NSORT,  Is  used  to  accurately  find 
the  maximum  amplitude  of  a  pulse  by  two  different  rapid  methods. 

If  NSORT  *  1  the  following  method  Is  used.  First  the  maximum  of  the  84  ampli¬ 
tudes  Is  found.  Then  another  search  Is  performed  to  find  the  relative  maxima,  these  are 
amplitudes  that  are  less  than  the  maximum  of  the  84,  and  greater  than  this  maximum  minus 
DA.  About  the  maximum  of  the  84  amplitudes,  and  about  each  relative  maxima  a  sin  X/X 
Interpolation  is  performed.  Tne  number  of  Interpolated  points  between  each  sample  Is  NIPC. 
This  Interpolation  Is  done  In  a  region  -  IDX  to  IDX  about  each  location.  From  these  Inter¬ 
polations  a  maxima  Is  found.  This  Is  the  maximum  amplitude  for  this  pulse.  Its  value  Is 
stored  In  the  array  AMXT,  and  a  return  Is  made  to  the  calling  program. 

If  NSORT  »  2  the  region  between  channels  NFT  and  NLT  Is  Interpolated  by  sin 
X/X.  The  number  of  Interpolations  between  channel  Is  given  by  NIPC.  The  maximum  of  these 
Interpolated  values  Is  then  the  maximum  for  this  pulse.  It  Is  stored  In  the  array  AMXT,  and 
a  return  Is  made  to  the  calling  program. 

The  calling  sequence  Is: 

CALL  AMPMAX  (NIPC,  NCH,  DA,  IDX,  NPPF,  AMP,  SXDXS,  NFT,  NLT, 
NSORT,  AMXT,  INP) 


DDD 

LLL 

IPT 


The  output  is: 

B 


where  the  Inputs  arc: 

NIPC  -  number  of  Interpolations  per  channel. 

NCH  -  number  of  channels 

PA  -  used  In  finding  relative  maxima  when  NSORT  ■  1 

1PX  -  used  In  interpolating  about  relative  maxima  when  NSORT  =1 

NPPF  -  number  of  pulses  per  frame 

AMP  -  this  ts  the  amplitude  array  that  Is  searched  for  a  maxima 

SXDXS  -  Slnx/X  array  used  In  Interpolation 

NFT  -  number  of  first  channel  for  Interpolation  when  NSORT  =  2 
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NLT  -  number  of  last  channel  for  interpolation  when  NSORT  *  2 

1NP  -  this  variable  must  be  one  the  first  time  AlriPMAX  is  called,  and 

incremented  by  one  every  following  call. 

the  output  is: 

AMXT  -  an  array  containing  the  maximum  of  each  pulse.  The  number  of 
values  in  the  array  is  equal  to  the  number  of  timer  AMPMAX  is 
called. 

(4)  SX0XS  -  Subroutine  for  producing  sinX/X  array. 

This  routine  sets  up  the  sinX/X  table  SXDXS  SXDXS(l)  =  1,  for  einX  X  at 
X  =  0.  Then  X  is  incremented  by  r/NIPC  until  there  are  (NCH)*(NIPC)  ♦  10  values  In  the 
SXDXS  array. 

The  calling  sequence  is 

CALL  SX0XS  (SXDXS,  N1PC,  NCH). 
where  the  inputs  are: 

N1PC  -  jf/NIPC  is  the  increment  for  X  used  in  calculating  the  array  of 

slnX'X. 

NCA  -  There  are  NCH* Nl PC  ♦  10  values  In  the  SXDXS  array, 

the  output  is: 

SXDXS  -  an  array  of  sinx/X  values. 
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i  DATA2X  Program  Data  Cards 
Card  1.  Format  110,  F10.  3,  2110 


NIPC 

DA 

IPX 

NSORT 


Card  2. 


-  number  of  interpolations  per  channel,  used  by  AMPMAX 

-  used  by  AMPMAX  to  find  relative  maximum  when  NSORT  *  1 

used  by  AMPMAX  in  interpolating  about  relative  maximum  when 
NSORT  =  1 

-  determines  which  method  is  used  by  AMPMAX  to  find  maximum  ampli¬ 
tude  of  a  pulse. 

Format  615,  4F10.  3 


MFORMT  - 

K 

LLL 

IPT 

NFT 

NLT 

DTAU 

C 


If  M FORMAT  »  1,  format  I  data  is  produced,  if  MFORMT  *  2  format  fl 
data  is  produced. 

there  are  K*1  values  in  the  RZ  array  produced  by  subroutine  RZEROM. 

there  are  LLL  values  in  the  B  array  produced  by  subroutine  BAMP. 

print  control,  if  IPT  4  0  the  RZ  and  the  B  array  produced  by  RZEROM 
and  BAMP  will  be  printed  out,  if  IPT  -  0  they  will  not  be  printed. 

This  is  used  by  subroutine  AMPMAX.  and  is  the  first  channel  to  be 
interpolated  when  NSORT  »  2 

This  is  used  by  subroutine  AMPMAX,  and  is  the  last  channel  to  be 
interpolated  when  NSORT  *  2. 

this  is  the  increment  for  the  X  variable  at  which  sin  Jx*  -  c2//x2  -  C2 
is  calculated  by  subroutine  RZEROM.  1  * 

determines  the  peak  to  sidelobe  ratio  of  the  response  pattern  produced 
by  the  RZEROM  routine. 


OMB 

DDD 

Card  3 

NPT 

M 

AK 

ROT 

PRF 

DELX 


-  scale  factor  for  the  variable  x  in  the  RZEROM  routine 

-  sidelobe  suppression  constant  for  the  BAMP  subroutine. 
Format  215,  5F10.  5 

-  number  of  point  scatterers 

"  *°r  *ach  target  2M+1  pulses  are  generated. 

-  increment  for  the  range,  for  84  channels,  AK  =  1. 

-  rotation  rate  of  point  scatterers  in  degrees  per  second. 

-  pulse  repetition  frequency  in  Hertz. 

-  actual  range  Increment  for  84  channels  AK  *  1 


X ST  ART  -  this  is  the  first  range  value  relative  to  center  of  rotation.  If  the  center 
of  rotation  of  the  target  is  to  be  in  the  center  of  the  84  channels  then 
XSTART  »  -42  (i.  e. ,  will  have  -42  to  +41  channel  locations,  the  center 
of  rotation  at  cham>el  0. ) 

Card  4.  Format  8F10.  3 

R(l),  R(2) . R(NPT)  -  This  is  an  array  of  the  ranges  of  each  point  scatterer  from 

the  center  of  rotation.  If  NPT  Is  greater  than  8  more  than 
one  card  is  required. 


Card  5.  Format  8F10. 3 

A(l),  A(2) . A(NPT)  -  This  Is  an  array  of  the  amplitudes  of  each  point  scatterer. 

If  NPT  Is  greater  than  8  more  than  one  caro  if.  required. 


Card  6.  Format  8F10.  3 


PHIS(l),  PHIS(2), . . . .  -  This  is  array  of  the  phase  change,  in  degrees,  on  reflec 
PlflS(NPfy~  lion  from  each  point  scatterer.  If  NPT  is  greater  than  8 

more  than  one  card  is  required. 


Card  7.  Format  8F10.  3 

TETS(l),  TETS(2), ...  -  This  is  the  array  of  the  angles  of  the  lines  drawn  from  each 
TETS^NPT)  point  scatterer  to  center  of  rotation  of  the  target  when  it 

is  at  the  M  ♦  1  pulse  (there  arc  2M+1  pulses  for  each  target). 
Each  angle  is  at  the  center  of  the  rotation  interval.  There 
are  M  pulses  to  each  side  of  TETS  for  each  point  scatter. 
The  figure  below  shows  the  target  configuration: 


LOCATION  OF  All) 
^Oat  PULSE  2M  ♦  I 


(The  total  angle  through  which  the  target  rotates  for 
2M+1  pulses  is  2a  *  (2M+1)  (ROT)  (1/PRF)  degrees) 
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Card  8.  Format  13 

IEND  -  If  IEND  1  0  read  another  set  of  cards  3  through  7  for  a  new  target  and 
a  new  set  of  pulses.  If  IEND  >  0  program  writes  an  END  FILE  II  on 
tape,  pauses  to  remove  tape  then,  exits. 

3.  SDS  010  PROGRAM 

The  SDS  910  program  accepts  pulses  of  rada  data  from  a  magnetic  tape  produced  on  the 
UNIVAC  1108  by  the  Optical  Processor  Program.  On  this  tape  are  records  of  either  2048  or 
1024  words  each.  Each  word  represents  a  puise  width  which  is,  in  normal  operation,  a 
single  sweep  on  the  CR  T.  The  SDS  910  reads  these  records  into  storage  and  then  outputs 
them  at  a  given  rale  to  che  D/A  converter  and  the  film  recorder.  In  addition,  line  start  pulses, 
line  stop  pulses,  and  camera  clock  frequencies  are  generated. 

The  SDS  910  can  accept  Information  from  magnetic  tape,  paper  tape,  or  an  on-line  type¬ 
writer.  Also,  control  of  a  program  can  be  initiated  by  means  of  breakpoint  switches  i  through 

4.  The  SDS  910  Program  is  read  in  from  paper  tape  by  the  standard  fill  procedure  (see 
sec*,  on  on  operation  of  the  film  recorder  for  more  detail  on  the  program  and  operation  of  the 
SDS  910).  Ail  breakpoint  switches  should  be  reset.  The  program  is  now  In  the  idling  mode. 

It  Is  necessary  to  provide  time  between  Line  Start  pulses.  This  is  done  by  setting  and  then 
resetting  breakpoint  1,  and  then  typing  in  a  four  digit  number  equal  to  the  desired  number  of 
milliseconds.  If  this  is  not  done  the  time  between  Line  Start  will  be  set  to  200  milliseconds 
(1.  e. ,  5  pulses  per  second).  At  the  present  time  the  number  of  milliseconds  between  line 
start,  pulses  must  be  greater  than  200  and  less  than  400. 

For  normal  operation  break -point  3  Is  now  set.  The  program  now  alternately  reads  a 
bio^k  off  of  magnetic  tape  and  outputs  it  at  the  already  ftxed  rate  to  the  CRT.  Line  start 
and  line  stop  pulses  are  also  generated.  At  the  end  of  each  frame  of  pulses  all  outputs  are 
Inhibited,  except  the  camera  clock  frequency.  If  breakpoint  3  is  reset  and  2  Is  set  further 
magnetic  tape  is  not  read.  The  program  repeatedly  outputs  the  last  block  that  was  read  from 
a  magnetic  tape.  All  timing  is  exactly  the  same  as  when  breakpoint  3  was  set. 


Flow  Chart  for  Optical  Processor  Pro^im. 


Flow  Chart  for  DATA2X  Program 
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A  film  recorder  for  high  resolution  radar  echoes  has  been  designed  and  fab¬ 
ricated.  It  employs  a  modulated  cathode  ray  tube  trace  which  is  imaged  onto  a 
moving  film.  Individual  pulse  echoed  are  written  across  one  dimension  of  the  film, 
corresponding  to  the  direction  of  the  CRT  sweep.  Successive  echoes  are  accumu¬ 
lated  side-by-side  in  the  orthogonal  direction.  The  data  is  recorded  in  a  format 
suitable  for  subsequent  processing  in  a  c<  herent  optical  system.  The  purpose  is 
to  obtain  doppler  resolution  by  integration  over  a  sequence  of  pulses  This  is  to 
be  carried  out  simultaneously  for  each  range  resolution  element,  giving  an  output 
with  two  dimensions  of  resolution. 

The  specifications,  design,  fabrication,  ard  check-out  of  the  film  recording 
equipment  are  described  In  this  report.  Simulated  data  has  been  employed  to 
verify  the  overall  performance  of  the  device.  A  review  is  given  of  the  dot.  frames 
recorded  and  the  test  results  obtained  from  them.  The  system  can  record  data  in 
two  different  formats,  to  allow  some  flexibility  in  implementing  the  coherent 
imaging  processor. 
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